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SYNERGETICS AND QUANTUM PHYSICS OF THE ALIVE 


RESONANT ACTIVATION OF METASTABLE NON- 
EQUILIBRIUM SYSTEMS 


V. А. СНІМАКОУ, S. Р. SIT'KO 


Scientific Research Center «Vidhuk», Kiev, Ukraine 


Abstract. In this article we consider different approaches to the study of peculiarities of escape 
of a non-equilibrium system out of. metastable states in a variety of physically different pótential 
configurations which may be of relevance to the fundamental concepts of quantum physics of 
the alive (Sit'ko, 1993). We consider also the dynamical behaviour of some noise-free systems 
driven by an external periodic input and show how the change of the phase state of different 
systems under the influence of periodic driving force may occur. It is shown that à coherent 
action on a system of both periodic and random forces may substantially increase the activation 
rates. For resonantly driven bistable system the novel features (an existence of metastable limit 
cycle and appearance of caustics on the phase plane) are found. We have carried out a series 
of MC simulations of the rate activation processes in the considered non-equilibrium systems. 


Such simulations yielded also a remarkable phenomenon — the resonance-like response to 
periodic input in the phase space of the system within a small interval of the external frequency 
values. 


I. INTRODUCTION 


A study of resonantly activated non-equilibrium systems may have a wide- 
ranging significance for detailed understanding of the onset and characterization of 
the physical aspects of microwave radiation action upon biological systems in 
particular and, more generally, may concern the fundamental questions of new 
exciting field of research of the modern non-linear science — Physics of the Alive. Its 
successful applications to the treatment and control of various somatic and 
psychoneurological diseases in humans were discussed recently by Sit’ko and 
Mkrtchian (1994). 

For the historical background of this new field of research connected with a 
phenomenon of the resonant response of living systems to the low-intensity 
electromagnetic radiation in mm-band range (EHF EM radiation), we refer the reader 
to the recent discussion by Sitk'o and Mkrtchian (1994) (see also, Adey and 
Lawrence (1984) for a revicw). The scientific background is well represented in a 
variety of well-known papers (Devyatkov et al., 1973; Webb, 1984; Frohlich, 1988; 
Sit'ko, 1993). Last decade has produced new insights and developments in the field, 
both in experiment and theory (Belyaev, 1992; Serikov, 1993; Chaly, 1993). 

The main observations which constitute the heart of the concept of Physics of the 
Alive are as follows: 

1. An extremely low level of the EHF EM radiation power density which is able 
to produce effective action upon living matter. 

2. A sharp-resonant nature of the interaction between EM radiation and a system 
subjected to this radiation. There is a line character of spectra of the EHF EM 
radiation action upon living organisms. 

A very important feature of living structures is the way of their functioning as an 
open system under permanent exchange of matter and energy with ambient medium. 
From the thermodynamic point of view they are highly non-equilibrium systems 
operating far from the equilibrium point. The corresponding approaches for the 
description of such systems are widely known (e.g., Glansdorff and Prigogine, 1971). 
But from the point of view of Physics of the Alive with its unusual peculiarities, 
further theoretical study as well as conceptual developments of the ideas concerning 
the non-linear interaction between biological systems and EHF EM radiation ARI 
needed. Therefore, we will try to consider here some mechanisms of statistical escape 
in non-equilibrium systems which might be challenging to this new ficld of research. 

Non-equilibrium driven systems in general may have several coexisting attractors 
taken AS the limit points, or limit cycles. The simplest limit set of the points in the 
state space of a dynamical system in two dimensions is a limit point, while a circular 
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limit set is called a periodic limit set, or a limit cycle. The attractive point has a two- 
dimensional set called its basin of attraction. The state зразе is decomposed into a set 
of basins. Their boundaries are called separatrices. 

Coexisting attractors become metastable in the presence of noise. In this case, 
the fluctuational transition probability between system attractors (or the rate of 
escape from a metastable state), is one of the most useful characteristic quantity 
describing its dynamics. However, because non-equilibrium systems lack detailed 
balance, one can obtain (even using substantial approximations) closed-form 
expressions for the probability density only when a system is considered as one- 
dimensional Markovian (Hanggi, Marchesoni, and Sodano, 1984; Hanggi and Jung, 
1988), or it may be reduced to the effective one-dimensional Markovian system 
(Risken, 1984; Graham and Tel, 1985; 1986). The main problem here is to obtain the 
Stationary probability density — stationary solution of a Fokker-Planck equation 
(FPE), which is always known only for equilibrium systems, but not for non- 
equilibrium ones. 

There are many well established approaches to the problem of escape from a 
potential well and, in general, to the problem of diffusion in multistable systems 
(Ludwig, 1975; Ben-Jacob et al., 1982; Caldeira and Leggett, 1983; Freidlin and 
Ventzel’, 1984; Graham and Tel’, 1986; Luciani and Verga, 1987; Leiber, 
Marchesoni, and Risken, 1987; Forster and Michailov, 1988; Fox, 1988; Grabert, 
1988; Grabert and Linkwitz, 1988; Klosek-Dygas, Matkowsky, and Schuss, 1988; 
McKane, 1989; Dykman, 1990). The escape rate (this notion was introduced by 
Arrhenius in 1889), at which the system moves from one metastable state to the 
other, was first treated in terms of the escape process that is governed by Brownian 
motion dynamics influenced by forces which in turn are connected with the 
temperature and friction in the pioneering paper of Kramers (1940). Kramers’ paper 
during the past two decades still remains the cornerstone source for many 
generalizations — rate theory with memory friction, the turnover between weak and 
strong friction (Hanggi, Talkner, and Borkovec, 1990), ets.— and developments in 
the non-equilibrium rate theory, both in classical (Risken, 1984; Luciani and Verga, 
1987; Weiss et al., 1987) and quantum mechanics (Caldeira and Leggett, 1983; 
Larkin and Ovchinnikov, 1984; Pollak, 1986; Grabert and Linkwitz, 1988). The 
sketch of these and some other developments will be given below. 

The situation becomes much more complex when we start to deal with resonantly 
activated processes because a system then is driven not only by thermal noise (or any 
other type of noise), but by a periodic force, too. The archetypical situation here is a 
non-equilibrium bistable stochastic system in the presence of noise and periodic input. 
It is well known that different non-linear systems driven by a periodic input may have 
several stable states that correspond to forced vibrations differing in their periods, 
amplitudes, and phases. For example, a frequency of oscillation of a non-linear 
oscillator subjected to a strong periodic field depends on their amplitude and such a 
system displays a bistability (Dykman and Smelyauskiy, 1990). 

Recent emphasis in this area is put around dynamical aspects of influence of both 
a coherent and oscillatory driving force on the escape from metastable wells (resonant 
activation of a Brownian particle out of a potential wells) (Carmeli and Nitzan, 1985; 
Chow and Ambegaokar, 1988; Mannella and Grigolini, 1989; Linkwitz and Grabert, 
1991; Gammaitoni et al., 1993). 

Closely related to the problem of resonant activation is the problem of so-called 
stochastic resonance (SR) when a transition probability between different system 
attractors increases exponentially with increasing noise intensity. An importance of 
SR for application to different systems in physics, chemistry and biology (neurons, 
membrane channels, etc.) is clearly attested by the fact that literature concerning the 
SR is growing very rapidly (McNamara and Wiesenfeld, 1989; Debuath, Zhou, and 
Moss, 1989; Presilla, Marchesoni, and Gammaitoni, 1989; Gang, G. Nicolis, and C. 
Nicolis, 1990; Moss, 1991) (see also the special issue of J. Stat. Phys., 1993 (Eds. 
Moss, Bulsara, and Shiesinger) that provide an overview of the field). The latest 
developments in this field concerning the possibility of existence of the high frequency 
SR (HFSR) (Dykman et al, 1993, 1994; Gammaitoni et al, 1993), as well as 
interesting results concerning the first passage time problem in time-dependent fields 
(Fletcher, Havlin, and Weiss, 1988) and cooperative effects in a coherently: activated 
system of coupled elements where unusually large amplification of the periodic 
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modulation is observed (Shiino, 1987; Jung et al., 1992), are also of special interest 
for the problems discussed here. 

Recently very encouraging results were obtained by Petracchi et al. (1994) in a 
study of periodic modulation of the membrane potential of K* channels which clearly 
indicate that a problem of SR might be important for the transmission of weak 
electromagnetic signals at the cellular level of organization of living systems. A 
sinusoidally modulated ionic channel in this study was considered in terms of 
multistable system coherently driven by internal noise in addition to the external pe- 
riodic force (cf. Chinarov et al., 1992). 

It is only in recent years that detailed investigation of the role of dissipation in 
non-linear dynamical systems has begun (see, for example, a discussion by Wolynes 
(1987) and references therein). Dissipation in condensed phase causes a process more 
likely to be adiabatic. Frequency dependence of the dissipation may lead to a 
profound influence on results (this is especially important in biomolecules where 
friction may be involved). 

In connection with the above mentioned necessity to look for the new 
developments of the mechanisms of interaction between biological systems and EHF 
EM radiation itself, it seems worth mentioning here that there exist a number of 
recent results in which it was observed the revelation of macroscopic quantum nature 
of the Alive. If living systems are able to respond to single quanta of flux in a 
magnetic field (its value is h/2e = 2.07 - 10735 Wb), they also have the 500 MHz/mV 
frequency/voltage Josephson interconversion and the Josephson effect may have 
relevance to them (Frohlich, 1988). In was also recently declared by Smith and Best 
(1990) that «Microwave effects in yeast cells strongly suggest that the Josephson 
relation (500 MHz/mV) should be applied to all experiments with living systems». 

It seems natural to treat some actual problems of resonant response of living 
systems to the EHF EM radiation in mm-band range, which still remain under 
discussion up-to-date, as well as mechanisms of appearance of characteristic eigen- 
frequencies in living matter (Andreyev, Bely, and Sit'ko, 1982), as closely connected 
to the study of resonant activation processes in metastable non-equilibrium systems. 

The aim of this article is to consider possible approaches to study the peculiari- 
ties of escape in a variety of physically different metastable potential configurations, 
which may have relevance to the fundamental concepts of «Quantum Physics of the 
Alive» recently discussed by Sit'ko (1993). We underline here that understanding of 
mechanisms of functional and stuctural organization of living matter is possible if one 
assume that there exist a highly non-linear interaction between the latter and EHF 
EM radiation. | 

The outline of the paper is as follows. In Sec. Il we consider the dynamical 
behaviour of some noise-free systems driven by an external periodic input and show 
how the change of the phase state of different systems under the influence of periodic 
driving force may occur. In Sec. III we discuss the problem of diffusion in multistable 
systems and. show that coherent action on a system of both periodic and random 
forces may substantially increase the activation rates. In Sec. IV, using the optimal 
path method, the fluctuational transitions in bistable systems are studied. For 
resonantly driven bistable system the novel features (an existence of metastable limit 
cycle and appearance of caustics on the phase plane) are found. A caustic is the 
boundary of the domain of the state space accessible to the fluctuational trajectories 
emanating from a given stable stationary state. Section V is devoted to the summary 
and discussion. 


П. DETERMINISTIC FLOW OF SOME SYSTEMS DRIVEN BY PERIODIC INPUT 


Firsr, we consider the model of a Brownian particle moving in a tilted non- 
confining potential and subjected to periodic force. In this case the system to be 
described is an ideal Josephson junction. With the forcing frequency closed to the 
oscillation eigenfrequency, a resonant microwave enhancement of activation was first 
observed experimentally (in the GHz-band region) by Devoret et al. (1984; 1987). 
The equation of motion of a_single macroscopic variable in this problem (phase 
difference across the junction Ф) is (Esteve, Devoret, and Martinis, 1986) 
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where U(@) is the tilted cosine potential, I,(t) is the noise current associated with the 
Josephson noise of the resistance, In(t) = Insin(Qt+ Фо) is microwave current, 
Фо = h/2e is the flux quantum, and К and C the are the junction resistance and 
capacitance, respectively. 


Fox further analysis, rewrite Eq. (2.1) using normalized units and omitting the 
noise current, in the following deterministic form 


Ф, Ф ОФ). 

E +7 "of =- aD + i,sin(Qt t Фо), (2.2) 
where Ф is dimensionless phase difference across the junction, 7 is dimensionless 
damping parameter (17! = 2е/. R? C/A), І, is critical current shunted by R and С and 
driven by dc bias Го, and і, = 1,/ I. (Barone and Paterne, 1982). 


Eq. (2.2) describes deterministic dissipative system which possesses steady-state 
solutions corresponding to its attractors. 


The potential U(®) which has local minima at the points 
Prax = insin (іо) + 2zn, is the so-called washboard potential (Fig. 1) 


U(®) = ij — i,,cos(P), (2.3) 
where ip = J,/J, is the normalized bias. 


Figure Е. Washboard potential described by Eq. (2.3). 


Consider the deterministic flow of this system. A dynamical system may be 
modelled on a state space by a vector assigned to each point. The curved arrows 
associated with each trajectory represent the flow of the dynamics — the simultaneous 
movement of all initial states along the system trajectory. 

We first discuss a deterministic flow of a system without periodic input (i, = 0). 
In this case for the values of a bias current i; satisfying the inequality ij > i, where 
і is some critical current value, the phase of a system comsists of regions of attracting 
domains of many coexisting focus-like attractors. These domains (with corresponding 
family of the phase trajectories for each attractor) are shown in Fig. 2a. For the values 
of i; larger than the current values which characterize coexisting, focus-like attractors, 
there exist in the phase space both domain of attraction of the running state and 
focus-like attractors with their own basins of attraction (Fig. 2b). For 
io > itr, (it > іс) the system has only monostable running state with the whole 
phase space being only the basin of attraction for this state (Fig. 2c). 

In order to illustrate what is involved when the periodic input is switched on 
(in # 0), let us examine the dynamical behaviour of a system for some certain values 
of the driving frequency © keeping unchanged other parameters. Figure 3 illustrates a 
remarkable phenomenon — an appearance of resonance-like response of the 
considered system to such input. For small values of ©, the phase space is practically 
the same as in the case without periodic input (Fig. 2a) — one can see there only 
coexisting focus-like attractors (Fig. 3a). For larger values of the driving frequency, 
the phase space is transformed in such a way, that besides the focus-like attractors 
there also appears monostable running state (Fig. 3b). Note that such transformation 
of the phase space takes place only within a small interval of the frequency values 
(Q, < Q « Q5) and, moreover, the bias current ip is constant. In contrary to the case 
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without periodic modulation (Fig. 2c), the structure of the phase space recovers again 
for the frequency greater than a certain value (Q > Q5) (Fig. 3c). And only in the 
sense that for some values of driving frequency a dynamical dehaviour of a system is 
characterized by essentially different phase portrain, we mean that it has resonance- 
like response. For the considered system, a value of the quality factor of oscillation Q, 
which is defined as a ratio of the cyclic frequency to the damping parameter, is less 
than 3.0. 


Figure 2. Determonistic flow of the system (2.2) 
without periodic input. a) Domains of coexisting 
focus-like attractors, Q = 0.1; b) running state 
attractor and focus-like attractors with their own 


Figure 3. Resonance-like response of the system 
(2.2) to periodic input. a) Small values of ©; 
b) phase space with larger values of the driving 
frequency (0) < Q < Q5); c) structure of the 


phase space for (Q1 > Q2). See also а discussion 
in the text. 


basins of attraction, © = 0.6; c) monostable 
running slate, Q = 1.2. Parameter y = 0.35. 


Next, we consider other archetypical model — a Brownian particle moving in a 
bistable potential. In the case when this system is driven by periodic input, its 
equation of motion represents the classical dissipating non-linear oscillator — the 
Duffing oscillator with a linear friction and is of the form (Chinarov, Dykman, and 
Smelyanskiy, 1993) 


924 94. 9004) 

o? +y dp 30 + 50с08(07) + ED, (2.4) 
where y is dimensionless damping parameter, potential U(q) has local minima at the 
points qu = £ 1 

4 


с 
О ру (2.5) 


and &(f) is a zero mean Gaussian random process with correlation function defined by 
equation 
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where D is the noise intensity. 

Figure 4 depicts a topology of the phase space of this system represented by two 
stable focuses which correspond to local minima qmim of the potential U(g) in the 
absence of noise and periodic input (a corresponding separatrix is also shown there). 
When the periodic force is switched on, a topology of the phase space changes: for low 
values of the driving frequency there appears a limit cycle which is shown in Fig. 4a, 
only one stable focus exists for some larger values of the driving frequency (see 
Fig. 4b), and for higher values of the driving frequency, there appear two stable 
focuses again as shown in Fig. 4c. A topology of phase plane in this case is very 
similar to that of the modulation of a free noise one-system moves practically along 
the deterministic trajectories existing in the absence of periodic forcing. The Q value 
for this system is varied in the range 0.7—10.0. 

In Figure 5 a topological structure of the phase plase is shown for different values 
of the damping parameter y and for small value of the driving frequency (о = 0.1). 
For low values of the damping parameter y (quality factor Q = 10 in this case), a 
strange attractor may appear (Fig. 5a). For larger values of y(Q = 0.7—4.0), complex 
oscillation patterns appear in the phase space. They are shown in Fig. 5b and Fig. 5c. 


верата Ст | y 
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Figure 4. Topology of the phase space of system 
(2.2) in the absence of noice and periodic input. 
a) Limit cycle for low values of the driving 
frequency (œw = 0.5); b) one stable focus for 
larger values of the frequency (0 = 0.05); c) two 
stable focuses for higher values of the driving 
frequency (w = 20). Parameter y = 1.0. 
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Figure 5. Topological structure of the phase plane 
for different values of the damping parameter y 
and for o = 0.1. a) Strange attractor for low 
values of the damping parameter y = 0.1 (quality 
factor О = 10); b) — co complex oscillation 
patterns for larger values of y = 0.25 (b), 
y = 0.75 (©). (О = 0.7 4.0, so = 0.75). 
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In Figure 6 quite unusual types of attractors describing the complex dynamical 
behaviour of the considered system on the phase plane are shown for different values 
of the damping parameter y and for the larger value of the driving frequency 
(о = 0.596). The structure of the attractor with small value of the damping parameter 
(y = 0.1), shown in Fig. 6a, changes substantially with increasing y, and for the value 
of about 1.5 and more, the oscillations are localized around only one of the system 
stable states (Fig. 6d). 


Figure 6. Attraciors describing the complex dynamical behaviour of the considered system on the phase plane. 
a) Structure of the attractor with small value of the damping parameter (у = 0.1); b) у = 0.25; с) y = 1.0; 
d) oscillations are localized around only one of the system stable states (y = 1.5). In all cases о = 0.596, 
50 = 0.75. 


In Figure 7 some attractors are shown for different values of the damping 
parameter у and the driving frequency w. 

It may be clearly scen from these figures that the dynamics of the deterministic 
systems driven by periodic force strongly depends on the quality factor Q and might 
be well controlled by the driving frequency. 
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Figure 7. Attractors for different values of the damping parameter y and the driving frequency о. 
а) w = 0.05, y = 1.0; b) w = 0.5; с) о = 0.05; d) о = 0.59; е) o = 0.59, y = 0.4. 


ПІ. DIFFUSION IN METASTABLE SYSTEMS 


First we discuss briefly a problem of diffusion in metastable systems where many 
approaches were suggested to evaluate the escape rates. The flux-over-population 
method allows one to calculate first the steady-state current and then to obtain a sta- 
tionary non-equilibrium probability density inside the domain of attraction (Hanggi, 
Talkner, and Borkovec, 1990). For the construction of the current with its density 
р(х; t), one must have, at least, in principle, the explicit solution of corresponding 
master equation which governs the time dependence of density (both in Markovian 
and non-Markovian cases). For the bistable system with archetypical potential U(g) 
described by Eq. (2.5), the total probability flux j over the barrier is given by the 
product of escape rate from the left well to the right one and the population nw of the 
left well, and the rate of escape is given by the following expression 


Ki rw = і піз (3.1) 


which for the case of a Brownian particle with mass М moving in the quartic potential 
U(q) reads as (Hanggi, Talkner, and Borkovec, 1990) 
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where q+ is the value of coordinate q which is larger than its value at the right well, 
y is the damping coefficient, kg is the Boltzman constant, and T is the absolute 
temperature. 

Another well known method is the mean-first-passage-time (MFPT) approach 
where MFPT is just the average time needed for the random walker, starting from 
some point inside the domain of attraction, to reach the separatrix for the first time 
(Matkowsky and Schuss, 1981; Fox, 1988; Fletcher, Havlin and Weiss, 1988; Hanggi, 
Talkner and Borkovec, 1990). This becomes possible if the probability for a recrossing 
of the separatrix decreases rapidly as a function of the distance from it, i. е., the 
backscattering of a trajectory from the regions beyond the separatrix is negligible 
(Kramers, 1940; Talkner and Braun, 1988). The escape rate is usually determined by 
the inverse mean first-passage time. In general, this is not an easy task to calculate 
this quantity, but fortunately, for systems driven by the white noise, there exists 
explicit quadrature formula for one-dimensional processes. Recently some theoretical 
results were obtained also for higher dimensions (non-Markovian processes with 
memory) (Hanggi, Talkner and Borkovec, 1990). 

The problem of diffusion in bistable systems driven by colored noise, where the 
main difficulty is connected with finite noise correlation time, was considered in a 
number of recent papers (Matkowsky and Schuss, 1981; Fox, 1987, 1988; Jung and 
Hanggi, 1988; Trisonis and Grigolini, 1988; Leiber, Marchesoni and Risken, 1987; 
Chinarov et al., 1991), for review, see also Hanggi, Talkner and Borkovec (1990) and 
references therein. To such non-Markovian systems, the formalism of the FPE is not 
applicable, because of the lack of a master equation which describes the evolution of 
the stationary probability density. Colored noise in this case affects not only the 
dynamics of a system, but also the form of the stationary probability density. Existing 
aproximate approaches work best for low and extremely high values of the noise 
correlation time t (noise color). For the former, t serves as a small parameter for 
evaluation of corrections to the white-noise limit (Fox, 1987; Sancho et al., 1982), 
while for the latter case, the dynamics is considered in an approximation replacing the 
diffusion description of one-dimensional system variable x, which is much slower than 
the noise variable ё, by a deterministic motion of coordinate x in some effective 
potential (Trisonis and Grigolini, 1988). Unfortunately, up to now there exists a 
number of contradictory results concerning the escape rates out of metastable states of 
systems driven by colored noise (Jung and Hanggi, 1988). Therefore, precise 
computational methods are strongly needed here. 

As to the computational problems, both in equilibrium and non-equilibrium rate 
calculations, different methods were developed. In one of them, in order to calculate 
the rates, evolution equation for the probability density is solved numerically using 
the eigenvalue expansion which for the onc-dimensional case has the form (Risken, 
1984) 


р(х; t) = bs P(x) exp (—4,), (3.3) 


where the sets {@(x)} and {A;} are the cigenfunctions and eigenvalues of the master 
operator in the master equation which governs the time dependence of the probability 
density. In practice, the basis set is truncated and the problem is reduced to the 
evaluation of the eigenvalue spectrum of some non-Hermitian matrix using the 
matrix-continued fraction methods. This approach was used to calculate the rates of 
escape out of a single metastable state, in a double well, and in a case of periodic 
potential (Risken, 1984). 

The other widcly used approach is based on the numerical integration of the 
stochastic equations of motion «Sancho-et al., 1982; Fox, 1988). In this case the 
average lifetime of the stochastic trajectory in a well is the inverse of the rate. The 
drawback of this approach is ihe complex separation of time scales, which can be 
overcome by using the so-called reactive flux method, where the rate is calculated 
using classical transition-state theory (Eyring, 1935) and some variational procedure 
(Montgomery, Chandler and Berne, 1979). 

In a general case of non-thermal systems, c. g., those driven by a time-depen- 
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dent field, a dissipation would be expected to influence not only prefactors, but also 
activation energies for the rate of the escape from metastable states (Wolynes, 1987; 
Chinarov, Dykman, and Smelyanskiy, 1993). The influence of the environment 
results in the damping of the dynamic system, what makes the corresponding 
Stochastic process to be a non-Markovian in nature. 

Several papers have appeared in recent years dealing with a problem of the 
dissipation influence on the probability of escape from a metastable state. Some of 
them are concentrated on the dissipation influence on escape processes in classical 
systems (Mel'nikov and Meshkov, 1986; Wolynes, 1987; McKane, 1989) (see also 
recently developed theory by Chinarov, Dykman, and Smelyanskiy (1993) for the 
calculation of the dissipative corrections to the arguments of the exponentials for 
escape rates of thermally non-equilibrium classical systems). Other papers are 
concentrated on the influence of dissipation on quantum tunneling (Hanggi, Weiss, 
and Risenborough, 1986; Pollak, 1986; Legget et al., 1987; Hanggi, 1987; Grabert 
and Ingold, 1988). This problem is actual for many physical systems, in particular, 
for underdamped Josephson junctions with coexisting fixed- and oscillatory-current 
regimes (Barone and Paterne, 1982; see also Sec. ID. 

A very useful quantity needed to evaluate the resonant activation of a Brownian 
particle moving in a tilted non-confining potential (cf. Eqs. (2.1) — (2.2)) is a so-called 
enhancement factor R (sometimes it is called an activation factor) introduced by 
Devoret et al., (1987) 


Tis, ©) 
~ TO) " (3.4) 


where Pm, Q) is a rate of escape from metastable well (f. e., a potential described by 
Eq. (2.3)), and Г(0) is a usual Kramers’ value (Kramers, 1940) in the absence of 
oscillatory driving force of amplitude i, and frequency Ф. When the microwave 
frequency is varied, there exist a pronounced peak for the enhancement factor R near 
the frequency of oscillation at the bottom of the well (Devoret et al., 1984; 1987). It 
was shown recently by Linkwitz and Grabert (1991) that the structure of this 
resonance (its asymmetric and quite broad form) may be stipulated by the 
anharmonicity of the potential. Some discrepancy still exists for the high frequency 
fallout — exponential: or power-law tail of the resonance in a case of a cubic 
metastable well (Larkin and Ovchinnikov, 1984; Mannella and Grigolini, 1989; 
Linkwitz and Grabert 1991). There is also a strong effect of a time-delayed friction on 
the rate of escape out of a metastable well (Grabert and Linkwitz, 1988). 
Corresponding non-Markovian damping kernel describes the influence of the 
transmission line shunting a current-biased Josephson junction on its dynamics 
(Esteve, Devoret, and Martinis, 1986). 

Recent numerous studies of SR (some references are given in the Introduction) 
also bare a relation to the problem of resonant activation (enhancement of the rate of 
activation processes by a resonant driving force). Stochastic resonance itself is a 
revelation of interesting non-linear phenomenon associated with noise induced 
switching in multistable systems which are driven also by а periodic forces 
(McNamara and Wiesenfeld, 1989; Presilla, Marchesoni, and Gammaitoni, 1989; 
Moss, 1991; Moss, Bulsara, and Shiesinger, 1993). One of the remarkable features of 
SR is an increase in the signal-to-noise ratio (SNR) with increased input noise (in 
contrast,.for systems having linear dynamics, an increase in noise must lead to a 
decrease in the output SNR). SR is not a resonance in its usual sense (increasing of a 
response when a driving frequency. is tuned to systems frequency), rather the SNR 
takes its maximum value.when the noise intensity is tuned near a certain value. 

The effect of a periodic force on the dynamics of a noisy multistable.system in 
the case of SR is stipulated, as a rule, by the additive external noise. Quite unusual 
situation may arise if the noise and periodic forces are considered multiplicatively 
(Dykman et al., 1992), In contrary to conventional low-frequency SR, in the last case 
SR occurs for asymmetric bistable potential with wells of different depths. Recently a 
new type of resonance was found in a case of overdamped bistable quartic oscillator 
driven by colored noise (Gammaitoni et al., 1993). Resonant crossing process here is 
characterized by frequencies much larger than the SR frequencies. 

Let us discuss now a problem of diffusion from the point of view of simulation of 
the rate activation processes. Below we list several non-trivial results of the Monte 
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Carlo (MC) simulations of the residence times (inverse of the rates) statistics which 
were obtained by us recently for the metastable systems considered earlier (Sec. ID. 

In Figure 8 the dependence of mean residence time (the time which the system 
spends in one of the wells before it reaches the potential barrier top) on frequency is 
shown for the potential described by Eq. (2.5). The driving periodic force amplitude 
multiplied by the amplitude of the oscillation in a well (define this quantity as Zp) is 
nearly equal to the value of the barrier height in this case. Two MC realizations 
averaged over 200 stochastic trajectories are practically identical. 

Figure 9 depicts a frequency dependence of the activation rate in a case when the 
quantity J; is much less than the barrier height. A case when the noise is the only 
driving force gives very low rates (see Fig. 9a). The periodic forcing with small 
amplitude taken alone can not switch the system between its stable states (Fig. 9b). 
Only a coherent influence of both noise and periodic forces produces significant 
enhancement of the escape rate (Fig. 9c). The value of a quality factor Q in this case 
is about 10. 


«T»^ К" k 


w 


Figure 8. Dependence of mean residence time on Figure 9. Frequency dependence of the activation rate 

frequency for the potential described by Eq. in a case when the quantity IF is much less than the 

(2.5). Averaging is done over 100 MC realizations. barrier height. See text for the explanation. a) Noise is 
the only driving force; b) periodic forcing with small 
amplitude; c) coherent influence of both noise and pe- 
riodic forces. Quality factor Q — 10. 


In Figure 10 a pronounced resonance peak (mean rate versus frequency) is 
shown fon the potential described by Eq. (2.5). The background rate level in the 
absence of periodic forcing is also shown herc. A quality factor Q here approximately 
equals 60. The resonance frequency corresponds to the real value of observed 
experimental values ~ 6.3 GHz for small values of the temperature (of the order of 
4—5 К) which gives the very low noise intensity (Devoret et al., 1984; 1987). 


R 


Figure 10. Dependence of the mean rate on the 
driving frequency for the potential described by Eq. 
(2.5). The background rate level in the absence of 
periodic forcing is also shown. Quality factor 
Q = 60. 
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IV. PRINCIPLE OF MINIMUM NOISE ENERGY (OPTIMAL PATH METHOD) 
AND FLUCTUATIONAL TRANSITIONS IN BISTABLE SYSTEMS 


For systems driven by a weak Gaussian noise, the analysis of the escape rates is 
based on the method of the optimal path (Ludwig, 1975; Ben-Jacob et al., 1982; 
Freidlin and Ventzel’, 1984; Luciani and Verga, 1987; Forster and Michailov, 1988; 
Dykman, 1988; Klosek-Dygas et al., 1988; McKane, 1989; Dykman and Smelyanskiy, 
1990). This method based on the Feynman's idea of the interrelation between the 
paths of the system and the driving noise (Feynman and Hibbs, 1975), makes it 
possible to reduce, to logarithmic accuracy, the calculation of escape probability of 
noise-driven system to the solution of a boundary-value problem for certain set of 
differential equations. An explicit form of these equations depends on the form of the 
initial equation of motion. It is substantial here that if the range of attraction to a 
stable state under consideration is limited by a separatrix with a saddle point on it, 
the optimal path along which a system escapes ends with an overwhelming probability 
at this point. 

А convenient approach to the calculation of the logarithm of statistical 
distribution (as well as the logarithm of escape probability) of a system without 
detailed balance is based on numerical analysis of the Hamilton-Jacobi equations, 
obtained by applying a WKB-approximation to FPE for the probability density 
distribution (Ludwig, 1975; Dykman, 1990; Chinarov, Dykman and Smelyanskiy, 
1993; Chinarov, Smelyanskiy, 1993). This approach is analogous to the non- 
equilibrium thermodynamic potential (Graham and Tel', 1985). It is worthy to note 
that such non-equilibrium thermodynamic potential playing the role of a Lyapunov 
functional, can be also written for the spatially extended pattern forming systems 
even in the case when a gradient form for the equations of motion does not exist 
(Graham and Tel’, 1991). 

The optimal path method may be also well suited for the calculation of escape 
probability in a case of resonantly activated processes in different physical (Dykman 
and Smelyanskiy, 1990; Chinarov, Dykman and Smelyanskiy, 1993; Chinarov and 
Smelyanskiy, 1993) as well as biological systems (Chinarov et al., 1992; Chinarov, 
1994; Chinarov, 1995). 

To study the noise influence on the system dynamics, let us write corresponding 
equations of motion (Eqs. (2. 1), (2. 4) in a more general form 


dq/dt = K,(q, p; t) + (q, pi 0 
dp/dt = K,(q, p; t) + &x(q. р; t), (4.1) 


where &(q, р; D is a Gaussian zero mean random field with a correlation function 
defined by the following equation 

g (4, p; 0 5 (4, p; (t— т))) = 2D6(t)B,(q, р). (4.2) 
Here functions K;(x; f) are some non-linear functions on their arguments x; which 
may depend on time 7. Below we consider only the simple case when the functions K; 
do not depend on time (full analysis with К; = КАТ) was done by Dykman and 
Smelyanskiy, 1990; Chinarov, Dykman 'and Smelyanskiy, 1993). The symmetry 
properties of the diffusion coefficients Ва, р) іп Eq. (4.2) follow the assumption that 
noise itself does not break the symmetry of the system with respect to time inversion, 
and therefore we assume that 


Bj (а, р) = By (4, T^ P)(- ту, (і, / = 1,2), 
Big Boy +. Bl, = 0, (Віз = Ва). (4.3) 


In a case of low intensity noise limit (D << 1), the system with overwhelming 
probability approaches one or another stable state, while the transitions between 
steady states practically do not happen. At larger times т (т >> и! in a vicinity of 
one of the stable states the quasi-stationary distribution is formed. For the times 
t >>t, under the large (and, respectively, with small probability) «outburst» of the 
noise, the system starts its motion in the vicinity of this stable state and approaches 
the point of observation X,(1). The problem we deal with consists of in the calculation 
of such a distribution. For the case of white noise in Eqs. (4.1), (4.2), the simplest 
approach to consider the transition probability density o((x;j]; f) for the stochastic 
process X; consists of in searching for the solution of corresponding FPE equation. 
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Let us assume that the distribution p(q, p; 0) is quasi-stationary in the range 
„<<1<< Wl, | (4.4) 


(т, is a relaxation time, т. = у | for the system described by Eq. (2.4), i. e., prior to 
the escape, the system gets thermalized inside the potential well and its intrawell 
quasistationary distribution p(q, p) becomes practically Gibbsian. Therefore, to the 
lowest order in the noise intensity D, the quasi-stationary solution for this equation 
can be found in the form р о exp (—R/D), where К = R(X;) satisfies the following 
non-linear equation 


H(R,, Ry; q, p) = 0, (4.5) 
where 
SITER 2 4.6) 
Н=5 С + 2В,,К,6, + В.К] + K R + KpRp, ( 
- OR _ IR (4.7) 


The function R(q, р) сап be associated with a mechanical action of а two- 
dimensional auxiliary system with the coordinates g, p, and Eq. (4.5) is then a 
Hamilton-Jacobi equation. The related Hamiltonian equations for the «coordinates» 4, 
p and «momenta» R,, R, describe the extreme paths of the system. The probability of 
the escape from a stable state f is given, to logarithmic accuracy, by the following 
expression (Dykman and Smelyanskiy, 1990; Chinarov, Dykman and Smelyanskiy, 
1993) 


W = const : exp[—R/D] ‚ R= Rds; Ps); (Ra; ру) = 0) 3 (4.8) 


where gs, Prs are the values of the coordinate X and momentum Р of the initial 
system at the stable state f and at the saddle point s. Function R satisfies the 
conditions: R(q, р) = 0, 04/041; = OR/óply-— 0; the point д, p; (on the phase 
plane) is the global minimum of R. 

. Further calculation of activation energies R is based on numerical analysis of the 
Hamiltonian equations for an auxiliary system (4.5) — (4.8) ог corresponding 
equations for the «coordinates» g, p and momenta R,, R, of this system 


дз 0H/AR,, p= 0H/OR,, К, = -дніду, Ву = -днідр, (4.9) 


where dots over different variables denote corresponding derivatives with respect to 
the time. 

To find the activation energy, it is mecessary to solve the boundary-value prob- 
lem for Eqs. (4.9) with q > q, p> реа! to апа д > чь P > pp al t> — (qr, pp and 
ds» p, are the values of g, р right at the metastable focus and the saddle point, 
respectively). Hence, the problem of calculation of the activation energy comes to 
finding the value of the parameter for which the extreme path arrives to a saddle 
point to a given accuracy. 

In Figure Па а typical example of fluctuational trajectory (the rotating frame is 
used here) which moves from the vicinity of state f (focus here) to the saddle point s 
is shown for the case of non-linear Duffing oscillator (Chinarov, Dykman and 
Smelyanskiy, 1993). This trajectory represents the most probable suitable path which 
corresponds to the optimal fluctuation of the noise. The equations for this trajectory 
are the Hamiltonian equations (Eq. (4.9)). 

Figure 116 depicts a fluctuational transition between different attractors of the 
oscillator along the optimal path. For resonantly driven bistable systems, the novel 
feature connected with a change of the topological structure of the phase — an exis- 
tence of metastable limit cycle — is found. From Figure 12 (see Fig. 12a and Fig. 
12b), it may be easily seen that metastable limit cycle appears exactly at the same 
place of the phase plane for fluctuational trajectories originated from both attractors of 
the system. 

Eluctuational trajectories can be calculated for different initial values of variables 
q, P, Ryp Rp surrounding the starting point of this trajectory — (е. g., focus-like 
steady state in Fig. Па). And only in this case we can obtain a bundle of trajectories 
which give, taking together, another remarkable qualitative feature — an appearance 
of topological pattern on the phase plane consisting of a family of caustics (Fig. 12b). 
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For systems described by FPE, the caustics associated with a saddle point have been 
first found by Chinarov, Dykman, and Smelyanskiy (1993). The caustics organize the 
system topology in such a way that due to the large fluctuation the system moves 
along the optimal path, say, from the vicinity of the first stable state, and there is no 
penetration of the diffusion process to certain domains of the phase plane surrounding 
the second stable state. Note that any trajectory becomes self-crossing after the first 
reflection from the caustic. 


Figure 11. Typical example of fructuational trajectory described by Eq. (4.9). a) Optimal path for the: nom- 
linear Duffing oscillator (f — focus, s — saddle point); b) fluctuational transition between different attractors 
of the oscillator along the optimal path. 


Metastable LC 


Figure 12. Structure of the phase space. for the auxiliary system described by Eqs. (4.5) —(4.9:. 
a) Metastable limit cycle; b) family of caustics (initial system is described by Eq (2.4)... 


Figure 13 depicts a topology of the phase plane near the vicinity of a saddle 
point. 


0.004 
Р | 
caustic 


Figure 13. Topology of the phase 
plane near the vicinity of a saddle 
point. 
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In Figure 14 an example of fluctuational trajectory is shown for the case of a 
Brownian particle which moves in a tilted non-confining potential described by Eq. 
(2.3). This particle is influenced here by a periodic force. Note that on the contrary 
to the deterministic flow on the phase plane (or corresponding stochastic trajectory 
when the noise is switched on), the auxiliary system described by Eq. (4.9) moves in 
the opposite direction comparatively to the initial system described by Eq. (2.4). 


e 


Figure 14. Fluctuational trajectory for the case of a Brownian particle noving in a tilted non-confining 
potential described by Eq. (2.3). 


V. SUMMARY AND DISCUSSION 


In this paper, we have considered various aspects of simultaneous influence of 
both the random and oscillatory driving forces on the non-equilibrium dynamics of 
multistable stochastic systems and have summarized the modern theoretical results 
concerning the resonant activation of the escape from metastable wells which may be 
closely connected to the problem of resonant response of living systems to the EHF 
EM radiation in mm-band range. In particular, the motion of a Brownian particle in a 
tilted non-confining potential as well as the dynamics of the classical dissipating non- 
linear oscillator driven by periodic force were considered. 

We have reviewed also some recent theoretical approaches suggested to evaluate 
the escape rates in non-equilibrium systems and discussed a problem of diffusion in 
bistable systems driven by colored noise and the influence of dissipation on the 
activation rates. The probability of the escape from metastable states was calculated 
using the instanton technique and the novel features of the dynamical behaviour of 
resonantly driven bistable systems were obtained. These features — a metastable 
limit cycle which. is characterized by quite unusual fine structure and a family of 
caustics for fluctuational trajectories which divide the phase plane on domains with 
different types of dynamic flow — show that the symmetry of the phase portrait of 
such systems gets essentially broken. The non-equilibrium thermodynamic potentials 
which can be constructed using this technique, are very convenient for study of 
instabilities which may lead to the pattern formation in spatially extended systems. 

We have carried out a series of MC simulations of the rate activation processes in 
the considered non-equilibrium systems. They showed that the chaotic behaviour was 
taking place for certain values of the driving frequency and system may posses a new 
stochastic regime (an asymmetrical strange attractor has appeared) as well as other 
regimes having a double spiral phase portrait. Such simulations yielded also a 
remarkable phenomenon — the resonance-like response to periodic input in the phase 
space of the considered system within a small interval of the external frequency 
values. 

It should be noted that the topological nature of system attractors fully describes 
its evolution in such a way that it may be clearly seen from the obtained phase 
portraits how in some cases trajectory curves repeat eight and the stable limit cycle of 
different form appears in the vicinity of the separatrix loop for some values of 
parameters, how a birth of a pair of asymmetrical limit cycles occurs, or strange 
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attractor appears. In general, such phase portraits may be considered as dynamical 
«fingerprints» of the system behaviour. 

It is worthy to underline here, that the further study of the dynamical behaviour 
of non-equilibrium classical and especially quantum systems subjected to periodic 
and/or noise inputs, as well as the search for the conditions which give rise to the 
appearance of the sequence of bifurcations of differing regimes leading to the regular 
or chaotic behaviour, is necessary in order to get a deeper insight into the 
mechanisms of self-organization of living systems. In this connection, we only remark 
that it has been recently proposed that the healthy living systems are chaotic in 
nature, and diseases make their behaviour regular, and even periodic (Goldberger, 
Rigney, and West, 1990). 

The main conclusion which may be made here is that the dynamics of non- 
equilibrium systems driven by periodic force might be well controlled by the driving 
frequency. Hence, a study of resonantly activated non-equilibrium systems may be 
very important for understanding of the physical mechanisms of microwave radiation 
action upon biological systems. 


РЕЗОНАНСНАЯ АКТИВАЦИЯ МЕТАСТАБИЛЬНЫХ НЕРАВНОВЕСНЫХ СИСТЕМ 
В. А. ЧИНАРОВ, С. П. СИТЬКО 


В статье рассматриваются различные подходы к изучению особенностей распада неравновесных систем 
из метастабильного состояния для физически различных потенциальных конфигураций, которые могут 
иметь отношение к фундаментальным концепциям физики живого (Ситько, 1993). Мы рассмотрели 
также динамическое поведение некоторых систем без шума, подверженных внешним периодическим 
воздействиям, и показали, как могут происходить изменения фазового состояния различных систем под 
влиянием периодических сигналов. Показано, что когерентные воздействия на систему периодических 
и непериодических сигналов могут существенно повлиять на скорость активационных процессов в би- 
стабильных системах. Для таких систем были найдены новые особенности (существование метаста- 
бильных предельных циклов и существование каустик на фазовой поверхности). Было проведено много 
вычислительных экспериментов с помощью метода Монте-Карло, моделирующих активационные про- 
цессы в рассматриваемых неравновесных системах. Найден эффект резонансного отклика на периоди- 
ческие воздействия в фазовом пространстве систем в малом частотном интервале. 


РЕЗОНАНСНА АКТИВАЦІЯ МЕТАСТАБІЛЬНИХ СИСТЕМ, ЩО ЗНАХОДЯТЬСЯ В СТАНІ 
НЕРІВНОВАГИ 


В. А. ЧИНАРОВ, С. П. СІТЬКО 


У статті розглядаються різні підходи до вивчення особливості виведення нерівноважних систем з мета- 
стабільного стану у випадку різних фізичних конфігурацій потенціалу, які відповідають фундаменталь- 
ним принципам квантової фізики живого (Сітько, 1993). Окрім цього, розглядається динамічна по- 
ведінка деяких систем, що не мають шумів та керовані зовнішнім періодичним сигналом. У цьому ви- 
падку показуться, як може відбуватися зміна фазового стану під впливом періодичної керуючої сили. 
Показано, як когерентний вплив на систему як періодичною, так і випадковою силами може суттво 
підвищити швидкість активації. Для резонансної керованих бістабільних систем отримані нові характе- 
ристики існування границі метастабільного стану і появи каустик у фазовому просторі. Досліджується 
серія МС-моделей ступеня активації процесів у розглянутих системах, які існують у стані нерівноваги. 
Моделі такого гатунку дають знаменний феномен -- резонансно-подібний відгук до періодичного зов- 
нішнього сигналу у фазовому просторі системи у невеликому інтервалі зовнішніх частотних значень. 
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FLUCTUATION EFFECTS IN STATIONARY OPEN SYSTEMS: SELFORGANIZATION AND 
FRACTAL NATURE OF NEW PHASE NUCLEI 


A. V. CHALYI 


Abstract. The general notion of stationary states for open systems is proposed with fluctuation 
Hamiltonians taken into account. Conditions of limit-cycle and focus solutions for the 
fluctuation models of self-organization are examined. Results of the approach to the problem of 
the new phase formation based on the ideas of fractal geometry and finite-size scaling are 
discussed. 


Одна из наиболее фундаментальньх проблем современного естествознания — 
проблема самоорганизации открьтьх систем на нетермодинамической ветви со- 
стояний, включая процессы появления новых фазовых признаков — продолжает 
вызывать к себе повышенный интерес. К ее решению привлекаются достаточно 
мощные теоретические методы нелинейной механики, кинетических моделей, 
физики фазовых переходов и другие [1, 2]. 

Известный синергетический принцип «порядок — через флуктуации» пред- 
полагает обязательный учет вблизи критических (бифуркационных) точек флук- 
туационных эффектов, вызванных взаимодействием флуктуаций характерных 
параметров порядка (управляющих мод) исследуемой системы. В связи с этим 
представляется актуальным создание и усовершенствование так называемых 
флуктуационных моделей процессов самоорганизации и зарождения новой фазы, 
которые в своей основе опираются на достижения современной теории‘ фазовых 
переходов (теории масштабной инвариантности, ренормализационно-группового 
подхода, метода коллективных переменных) [3—6]. 

Продолжая исследования [7, 8], в данной работе: предполагается обсудить 
следующие важные вопросы: 

1) обоснование возможности использования флуктуационных гамильтониа- 
нов теории фазовых переходов для открытых систем, которые являются стацио- 
нарными в обобщенном смысле, 

2) получение условий возникновения периодических решений типа фокуса и 
предельного цикла для флуктуационных моделей процессов самоорганизации, 

3) привлечение идей фракталов и скейлинга в ограниченных системах к про- 
блеме структурообразования (в частности, для более глубокого понимания осо- 
бенностей процесса зарождения новой фазы). 


ОБОБЩЕНИЕ ПОНЯТИЯ СТАЦИОНАРНОСТИ ОТКРЫТОЙ СИСТЕМЫ Й 
ПО ПРИГОЖИНУ НА СЛУЧАЙ ОБРАЗОВАНИЯ УПОРЯДОЧЕННЫХ 
СТРУКТУР И САМООРГАНИЗАЦИИ 


Рассмотрим открытую систему стационарного типа. Как известно, понятие 
стационарности по Пригожину предполагает, что полная энтропия 5 „олн Открытой 
системы неизменна в силу того, что увеличение энтропии AS;, связанное c необ- 


ратимыми диссипативными процессами внутри системы, полностью компенсиру- 
ется взаимодействием системы с окружением, которое сопровождается оттоком 
энтропии AS, через поверхность, ограничивающую объем системы. Таким обра- 


зом, имеет место соотношение 
ASnonu = AS; + AS, = 0, d) 
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‚ поскольку AS; > 0, a AS, = — AS; < 0. Здесь 

AS; = ffo dV dt» 0 (2) 
Й 


-- положительное производство энтропии внутри системы, a 
AS, = - f f Js dQ dt < 0 (3) 
Q 


— приток отрицательной энтропии (негэнтропии) из окружающей внешней Cpe- 
ды. В (3) использованы стандартные для неравновесной термодинамики обозначе- 
ния, а именно: J, = 425,/40 dt < 0 — поток энтропии через поверхность ©, Ha- 
правленный внутрь системы, а dQ — векторный элемент поверхности, совпадаю- 
щий по направлению с направлением внешней нормали. 

Другими словами, стационарность достигается за счет такого притока нег- 
энтропии извне (AN, > 0), который в точности равен уменьшению негэнтропии 
внутри системы (AN; < 0), так что 


ANnony = AN; + AN, = 0. (ta) 


Сказанное выше подтверждаст известное утверждение: открытая стационарная 
система питается отрицательной знтрописй, или негэнтропией. 

Именно для подобных стационар- 
ных открытых систем выполняется тео- 
рема Пригожина о минимуме произ- 
водства энтропии, которая доказывает- 5 
ся для небольших отклонений системы 
от равновесия. Термин «небольшие 
отклонения» предполагает справедли- 
вость линейной термодинамики необ- 
ратимых процессов, т. е. такого при- 
ближения, для которого верны линей- 
ный закон J; = Х І. Ху между потоками 
J; и термодинамическими силами X, и 
задание изменения энтропии в виде 
квадратичной формы термодинамиче- 
ских сил AS = Х gig X;X,, где Як — KO- 
зффициенть, вычисленные в положе- 
нии равновесия. 

Обобщим понятие открытой систе- 
мы стационарного типа по Пригожину 
на случай, когда в произвольной от- 
крытой системе возможно создание 
упорядоченных (пространственных, 
временных или пространственно-вре- 
менных) структур. Для наглядности обратимся к рис. 1, на котором изображены 
различные фазовые состояния однокомпонентной системы классическая жид- 
кость — газ в координатах давление Р — плотность р. Всю фазовую плоскость 
можно разделить на следующие области: 

1 — область стабильной системы жидкость — пар в близкой (флуктуацион- 
ной) окрестности критической точки К, в которой температурная и плотностная 
переменные т = (T — Т.)/Т. ««1 и Др = (p —p.)/p. << 1 малы настолько, что 
флуктуационные эффекты могут быть описаны на основе «нулевого» приближе- 
ния современных теорий фазовых переходов 2-го рода и критических явлений 
(простого скейлинга) [3]; 

2 — область широкой окрестности критической точки, гдет = 1, Др = 1 ине- 
обходим учет неасимптотических и асимметричных поправок к масштабным за- 
конам «простого скейлинга» [9]; 

3 — регулярная область изменения физических свойств вдали от критической 
точки, где T >> 1, Ao >> 1 и не следует учитывать флуктуации параметра no- 
рядка; 

4—8 — области метастабильных состояний, классификация которых предло- 
жена в [10, 111, 

4 — нефлуктуационная область вблизи линии фазовых переходов 1-го рода 
(кривой сосуществования фаз, или бинодали); 
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5 — переходная (кроссоверная) область от нефлуктуационной к флуктуаци- 
онным областям 6—8; 

6 — область критических флуктуаций, связанных с близостью метастабиль- 
ного состояния к критической точке, где необходим учет флуктуационных эф- 
фектов подобных тем, которые имеют место в области [; 

7 — область гомофазных флуктуаций, связанных с близостью к границе 
устойчивости вдали от критической точки (так называемая псевдокритическая 
область); 

8 — переходная (кроссоверная) область от флуктуационного критического к 
флуктуационному псевдокритическому поведению; 

9 — область неустойчивых в термодинамическом смысле состояний, где про- 
исходят процессы спинодального распада и образования упорядоченных структур 
(самоорганизации). 

Именно эта, последняя, область и является предметом данных исследований. 
При образовании упорядоченных структур в области 9 изменение энтропии внут- 
ри системы AS; состоит в общем случае из двух вкладов: 1) положительного изме- 
нения энтропии AS(Ó за счет диссипативных процессов и 2) отрицательного изме- 


нения энтропии ASÍ9 за счет образования упорядоченных структур, т. е. 
AS; = AS + 4509. (4) 
При этом стационарное состояние может быть обеспечено двумя способами: 


a) если AS(Ó > 45), то AS; > 0 и для выполнения условия AS,,4, = 0 необходимо 
соблюдение неравенства 


AS, = — IAS(Ó — AS(9 | < 0, (5 
сходного с условием (1) стационарности открытой системы по Пригожину; 


б) если ASM < {AS |, то AS; < 0 и в этом случае необходим приток энтропии 
извне (AS, > 0) для обеспечения условия стационарности (А5поли = 0), т. e. 


AS, = IAS(9 | — AS > 0. (6) 


Неравенство (6) назовем условием стацонарности открытых систем в обоб- 
щенном смысле. Заметим, что для произвольной открытой системы может ока- 
заться неєправедливыми условие малости отклонений параметров от положения 
равновесия (малости градиентов, определяющих величину термодинамических 
сил) и тем самым условия применимости линейной термодинамики необратимых 
процессов и, соответственно, теоремы Пригожина о минимуме производства энт- 
ропии в стационарном состоянии. 

И, наконец, необходимо еще раз подчеркнуть, что при рассмотрении про- 
цесса образования упорядоченных пространственных или временных структур 
вблизи критических точек, т. е. при бифуркациях Тьюринга или Хопфа, необхо- 
дим последовательный учет флуктуационных вкладов А5. Поэтому 


ASnom = AS + 459 + AS, + А5. (7) 
Отсюда следует, что при выполнении условий стационарности (5) или (6) опреде- 


ляющим для энтропии становится слагаемое Д5;. Соответственно, для флуктуа- 
ционной части свободной энергии такой системы имеем: AF, = — Т Аби. 


УСЛОВИЯ ВОЗНИКНОВЕНИЯ ПРЕДЕЛЬНОГО ЦИКЛА ИЛИ ФОКУСА ДЛЯ 
ФЛУКТУАЦИОННЫХ МОДЕЛЕЙ ПРОЦЕССОВ САМООРГАНИЗАЦИИ 


В теории фазовых переходов и критических явлений флуктуационную часть 
свободной энергии АХ; часто называют флуктуационным гамильтонианом H [р] 
(например, гамильтонианом Ландау — Гинзбурга [3, 4]). Типичная форма тако- 
ro флуктуационного гамильтониана для двух взаимодействующих параметров по- , 
рядка имеет следующий вид: 
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AF; [p1 9217 Н (pi, 21 = (1/2) f lap? + ар? + (1/2) bipi + 
+ (1/2) 6563 + с. (Уфі)? + c2(Vp2)? 1 dV + Ні Фі» P2); (8) 


где Hin ЇФі, 92] — гамильтониан взаимодействия. В зависящей от времени теории 
фазовых переходов задание гамильтониана (8) полностью определяет систему 
уравнений движения (кинетических уравнений) для параметров порядка P, и P2 
на основании известного соотношения 


pir, t) = - ГӘН Ір дра» t). (9) 


Здесь флуктуационный гамильтониан H [y;] рассматривается как функционал OT 
параметров порядка q;(r, t), а Г; — коэффициенты Онзагера. 

Ранее в [6] были предложены несколько флуктуационных моделей процессов 
самоорганизации, т. е. моделей, использующих флуктуационные гамильтонианы 
типа (8). Для этих моделей был изучен тип особых точек и сформулированы 
условия образования пространственных и временных структур (бифуркаций Тью- 
ринга и Хопфа). Основные результаты этих исследований, касающиеся наиболее 
интересных случаев периодических решений в виде фокусов и предельного цикла 
можно сформулировать следующим образом: 

Флуктуационные модели с двумя взаимодействующими параметрами поряд- 
ка характеризуются системой уравнений движения (9) общего вида 


Pi = (рі, P2) + РАД, 
р. = O(P1, 95) + Р,Ар›, (10) 


me Фи Ө — нелинейные функции параметров порядка. 
Осциллирующие решения для системы уравнений (10) появляются, если 
определитель 
< 
зії вв $4 aD 


где 

Е = D, inm WV. E; = Dyk? — Ө,, Ез = V * Q,'. 
При дополнительных условиях E, + E; >< 0 реализуются решения в виде ста- 
бильного и нестабильного фокуса. 


Для периодических решений в форме предельного цикла необходимо и доста- 
точно, чтобы выполнялись следующие условия: 


а) E, + Е, = 0, Е, + Е, #0, Е, – Е.Е, < 0; 


6) Е + Е, = 0, р, = р, = 0, Е, * Е, #0, р, я D є 0, 
— IUsl — ПО ОТ 7 < Dik? < — Ия + ПО + о, 
= E Ty = aed па | 1 9|. 
0 Ө,’ ©, , 01 0 1 , 02 ©,’ e, , 
B) Е + £,=0, Е, # Е, є 0, lUo l = 1001, Е, < 0, 
LU, — [-Ез |? < БА? < IU, + [-Е |; 


Е + Е. =0, E, # E; #0, Di D, [Dk — 1Uol «0; 


m Е, = Е, = 0, Е, # Е, # 0, D, = 0, Р, = 0, 
Dok? l Uol = LU, « 0; 


е) E, = Е, = 0, Е, <0. (12) 


Не останавливаясь на детальном анализе полученных условий (12), отметим 
самый главный результат: периодические решения системы (10) появляются в 
области 9 на рис. І. Так, условие Е; < 0 означает, что произведение обратных 
восприимчивостей есть величина отрицательная, т. е. по одному из параметров 
порядка исследуемая система находится под спинодалью. Перевод системы под 
спинодаль может достигаться также по механизму восходящей диффузии, т. е. 
при отрицательном значении одного из коэффициентов диффузии D; или D». 
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Другими словами, термодинамически неустойчивая в зтой области система стано- 
вится структурно-устойчивой благодаря появлению пространственньх (времен- 
ных или пространственно-временных) структур. Действительно, если в пределе 
К, © > 0 восприимчивость системы по одному из параметров порядка была отри- 
цательной, то при наличии упорядоченных структур она становится неотрица- 
тельной, а именно: x(k) = x(0) + D? > 0, х(о) = у(0) + Dw? > 0, где D, и D— 
коэффициенты, характсризующие эффекты пространственной (временной) дис- 
персии, т. е. пространственную (временную) диффузию. 

Численное моделирование, проведенное в [12], подтверждает полученные 
результаты. Естественно, что сделанные выводы касаются лишь достаточно близ- 
кой окрестности особых точек (стационарных состояний фу), поскольку для TIO- 
лучения характеристического уравнения использовалось линейное приближение 
по возмущениям др; =P; — фо. Дальнейшая временная эволюция полученных 
упорядоченных структур требует дополнительных исследований с выходом за 
пределы такого линейного приближения. 


ИДЕИ ФРАКТАЛЬНОЙ ГЕОМЕТРИИ И СКЕЙЛИНГА ДЛЯ ОГРАНИЧЕННЫХ 
СИСТЕМ В ПРОБЛЕМЕ ОБРАЗОВАНИЯ ЗАРОДЫШЕЙ НОВОЙ ФАЗЫ 


Флуктуационные эффекты играют важную роль в процессах зарождения но- 
вой фазы. Их следует принимать во внимание как при гомогенном или гетероген- 
ном зародышеобразовании в метастабильной области, особенно в близкой окрест- 
ности границы устойчивости (области 6—8 вблизи спинодали), так и при фазовом 
расслоении термодинамически неустойчивой системы по механизму спинодально- 
го распада в области 9. 

На ранних стадиях процесса образования новой фазы зародыши не имеют 
сферической формы, а являются фрактальными объектами. Как показано в [10], 
флуктуация числа частиц в зародыше новой фазы непосредственно определяется 
фрактальной размерностью такого кластера. Действительно, вблизи критических 
(бифуркационных) точек восприимчивость у становится сильно флуктуирующей 
величиной, вследствие чего средне-квадратичная флуктуация числа частиц 


«AN?» V? ~ (N2/V kg Ty)? ~ V? p(kg Ty)? (13) 


обнаруживаєт достаточно сложную зависимость OT линейного размера L кластера. 
Действительно, вдали от критических (бифуркационных). точек «AN?» V? ~ [42, 
где 4 — пространственная размерность системы. Учет флуктуационных эффектов 
меняет эту зависимость радикально. В силу гипотезы скейлинга для простран- 
ственно-ограниченных систем [13, 14] (зародыши новой фазы являются именно 
такими) для восприимчивости может быть записана следующая формула: 


y= Dy (AzLV", ВЫВ), (14) 


гдс у (х, у) — масштабная функция, А и В — неуниверсальные константы, 
D. у, у — критические индексы, а A — внешнее поле, сопряженное параметру по- 
рядка. Заметим, что именно через величину А может быть исследован интересный 
вопрос о воздействии внешнего электромагнитного излучения на процесс зароды- 
шеобразования (к этому вопросу предполагается вернуться в дальнейшем). 

На основании формул (13) и (14) в окрестности критических точек имеем 


«ANS M2 ~ ДД ~ pP (15) 

или 
«AN? 217 ~ EPA, (ATL, ры”), (15а) 
где р — фрактальная размерность. При получении (15) использованы скейлинго- 


вые и гиперскейлинговые (содержащие пространственную размерность) равенства 
между критическими индексами 


у= (2-”у, а+2В +у=2, (и = 2-а, (16) 
а также известная формула Мандельброта, связывающая фрактальную и про- 
странственные размерности: D = d — p/v. 


Таким образом, из формул (15) и (15а) следует, uro масштабное преобразо- 
вание средне-квадратичной флуктуации числа частиц в зародыше новой фазы ха- 
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рактеризуєтся критическим индексом, равньм в точности фрактальной размерно- 
сти кластера. Типичнье значения фрактальной размерности кластеров вблизи 
границы устойчивости D = 2,5 [10]. Выводы о фрактальной природе кластеров — 
зародышей новой фазы подтверждаются рядом исследований, основанных на ком- 
пьютерном моделировании (см. [10, 15 p. 

Исследование динамики роста фрактальных кластеров для систем со скаляр- 
ным параметром порядка и центрально-симметричным межмолекулярным потен- 
циалом межмолекулярного взаимодействия показывает, что зародыши новой фа- 
зы достаточно скоро принимают сферическую форму. Их дальнейшая судьба за- 
висит от того, достигли ли они критического размера Lyp или нет. Если L < Lyp, TO 
флуктуационно возникающие докритические зародыши новой фазы рассасывают- 
ся. Если же L > Lyp, то такие закритические зародыши начинают неудержимо på- 
сти, поглощая в себя все окружающее вещество. 

Важным примером, иллюстрирующим подобный механизм зародышеобразо- 
вания, является, по-видимому, онкогенез. Используя универсальные представле- 
ния современной физики фазовых переходов, можно сформулировать необходи- 
мые условия, выполнение которых способно воспрепятствовать процессу неудер- 
жимого роста патологических новообразований, а именно: 1) увеличение 
значения критического размера L4, зародыша новой фазы, перерастающего B 
злокачественную опухоль, 2) создание более высокого нуклеационного барьера 
Rwin(Lep), который нужно преодолеть для последующего роста закритических 
зародышей. 

Покажем действенность выполнения этих условий на хорошо изученной 
двухфазной системе «жидкость — пар». Для этого случая известно выражение 
для минимальной энергии, необходимой для создания сферического зародыша 
радиуса L [16]: 


Ки L) =. (4/ Зулі?дидрк + 4n0L?, ат 


где ôu =u(p, Т) - н(ркоь Т) — степень пересыщения (перегрева метастабильной 
фазы пара или жидкости, др, = [Pxc(T) — p.]/p. — значение параметра порядка 
на кривой сосуществования (KC), pxc(T) и р. — плотности вдоль КС и в критиче- 
ской точке, с — коэффициент поверхностного натяжения. Первое и второе слага- 
емые в (17) определяют соответственно объемный и поверхностный вклады в ми- 
нимальную работу образования зародыша новой фазы. На основании (17) легко 
рассчитать критический размер зародыша и минимальную работу, необходимую 
для создания дееспособного зародыша критического размера, т. е. нуклеационный 
барьер: 

Lyp = За/дидрко» 

Roin (Exp) = Ібло"/ З(ди)? (OP xc)?- (18) 


Как следует из (18), вдали от критической точки и границь устойчивости (B 
области 4) возрастание величин Lyp и Rypin(Lxp) связано C а) увеличением с и 
б) уменьшением ди и др (на самой KC Ly, > ® Rmin > © ). По мере приближе- 
ния к критической точке (в области 6--8) имеет место такая ситуация: а) крити- 
ческий размер зародыша растет как радиус корреляции флуктуаций &: [р ~ È, 
поскольку с ~ т(2— 1)”, ди ~ 18, др тб, Від +1) =у + 28 = 2 – а = dv. Отсюда 
следует, что с приближением к критической точке размер зародыша растет по 
скейлинговому закону Г. ~ t " в пространственно-неограниченной системе или 
стремится к характерному линейному размеру L пространственно-ограниченной 
системы; б) нуклеационной барьер не имеет особенности в критической точке или 
на границе устойчивости, поскольку Rmin 7 [2 ~ 109—1)» E? ~ 19.. Однако, как по- 
казано в ряде работ (см. [10, 15 p, минимальная работа по образованию зароды- 
ша новой фазы Киї 1/Gi, где Сі — число Гинзбурга. Для систем с дальнодей- 
ствующим потенциалом, для которых справедливо приближение самосогласован- 
ного поля число Гинзбурга мало (Сі << 1) и нуклеационный барьер оказывается 
достаточно высоким (Rmin > > 1). 

Таким образом, процесс образования новой фазы энергетически невыгоден 
как на кривой сосуществования фаз, так и даже при глубоких вторжениях в ме- 
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тастабильную область для систем с малым числом. Ги нзбурга (жидкие кристаллы, 
полимерные соединения с лальнодействующим межмолекулярным взаимодейст- 
вием и др.). Очевидно, что развитие принципиальной идси изоморфизма фазовых 
переходов в физических и медико-биологических системах является актуальной 


проблемой, требующей дальнейших теоретических и экспериментальных иссле 
ваний. 


ФЛУКТУАЦІЙНІ ЕФЕКТИ В СТАЦІОНАРНИХ ВІДКРИТИХ СИСТЕМАХ: 
САМООРГАНІЗАЦІЯ І ФРАКТАЛЬНА ПРИРОДА ЗАРОДКІВ НОВОЇ ФАЗИ 


O. В. ЧАЛИЙ 


Введено узагальнене поняття стаціонарності відкритих систем, для опису яких можуть бути викори- 
стані флуктуаційні гамільтоніани теорії фазових переходів. Розглянуті умови отримання розв'язків у 
вигляді граничного циклу і фокусу для флуктуаційних моделей процесів самоорганізації. Обговорені 
результати використання ідей фракталів і скейлінга для обмежених систем до проблеми утворення 
зародків нової фази. 


ФЛУКТУАЦИОННЫЕ ЭФФЕКТЫ В СТАЦИОНАРНЫХ ОТКРЫТЫХ СИСТЕМАХ: 
САМООРГАНИЗАЦИЯ И ФРАКТАЛЬНАЯ ПРИРОДА ЗАРОДЫШЕЙ НОВОЙ ФАЗЫ 


А. В. ЧАЛЫЙ 


Введено обобщенное понятие стационарных открытых систем, для описания которых могут быть 
использованы флуктуационные гамильтонианы теории фазовых переходов. Рассмотрены условия 
получения решений в виде предельного цикла и фокуса для флуктуационных моделей процессов 
самоорганизации. Обсуждены результаты применения идей фракталов и скейлинга в ограниченных 
системах к проблеме образования зародышей новой фазы. 
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Abstract. The processes of self-organization in neural networks with homogeneous and 
heterogeneous structures of synaptic connections between neurons are studied and the dynamic 
behaviour of neural network driven by exernal noise is analyzed. The Wilson-Cowan type 
equations for the excitatory and inhibitory activities with some modifications are used. The al- 
gorithm for searching of the optimal trajectory which connects distinct basins of attraction is 
proposed. This trajectory starts at some point (its coordinates are determined by the procedure 
based on the minimum action principle), near the stable state passes near the saddle point and 
gets into the vicinity of the other stable state. The transition probability is determined by the 
flux over the saddle point. 


1. Introduction 


Several papers have appeared in recent years dealing with the problem of noise 
influence on the dynamics of neural system (Ahn and Freeman 1975; Kryukov 1984; 
Gorse and Taylor 1989; Chinarov and Degtyarenko 1990; Hendin et al. 1991; 
Chinarov et al. 1991). In contrary to gradient methods realizing associative memory 
organization (Kohonen 1984) and using the feedback mechanism for synaptic 
alterations (Freeman 1979), relatively low synaptic fluctuations without feedback may 
provide the function of mechanisms of successive recall of ihe stored information 
(Tsuda et al. 1987). The important characteristics of dynamic systems with coexisting 
attractors are probabilities of transitions between them. Usually these attractors are 
associated with memories stored by the neural network (Tsuda et al. 1987; Grossberg 
1988; Ritter and Schulten 1988). 

It was shown earlier that transitions between attractors in neural networks with 
excitatory and inhibitory interactions (both fixed point and limit cycle types of 
aitractors) driven by stochastic inputs may change substantially the level of self- 
organization of the network (Chinarov et al. 1991). In some cases low intensity noise 
may result in considerable transformations of the phase plane of the system and even 
in the appearance of stochastic oscillations. The noise may also stabilize the processes 
of pattern recognition in the presence of feedback. 

Information generated and processed by neural systems is expanded in space 
and time. The averagee number of inter-synaptic junctions needed for processing 
information is approximately ten (Szentagothai and Arbib 1975). It is clear, that due 
іо such high connectivity of neural nets, one may treat them as dynamic systems 
characterized by the ignorance of initial conditions and possessing the stochastic 
nature of generation of impulse activity. 

Studying the properties of stochastic effects with or without feedback in the 
network may have important functional significance for the understanding of 
mechanisms of the processing information by neural networks. The questions which 
arise here are as follows: what type of interaction between different neural structures 
may ensure the activity phases synchronization for providing the necessary control 
function for the system, how must corresponding neural network for generation of the 
activity pattern with preassigned parameters self-organize, what is the nature of de- 
pendence of transition probabilities between attractors in neural networks on the 
system parameters, and whether they are in accessible domains of state space of a 
dynamical system. 

In this work we study the processes of self-organization in neural networks which 
appear both in spatially homogeneous and heterogeneous structures of synaptic 
connections between neurons. Such self-organization arises due to the spatial and 
temporal character of nonlinear interaction between elemenis. The proposed scheme 
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eg es ee а По але. 
of network structural organization can provide generation of stable rhythmic activity. 
The model possesses the temporal and phase characteristics which exist in the real 
Systems such as movement activity generators and can explain phase dependent 
switching over their parameters under the influence of signals with various sensory 
modality. 

The dynamic behaviour of neural network driven by external noise is also 
analyzed. The problem is actual for many neurophysiological systems, in particular, 
for olfactory bulb of mammals (Ahn and Freeman 1975; Freeman 1979; Li and 
Hopfield 1989) and system controlling movement activity (Matsuoka 1987; Chinarov 
and Degtyarenko 1990). As a base model, we have used the Wilson-Cowan (Wilson 
and Cowan 1972) type equations for the excitatory and inhibitory activities with some 
modifications. 


2. Neural network model 
Let us consider the neural network with its dynamic state to be described by a 


point on the phase space of a system. The evolution of this point (system’s trajectory) 
is described by the following system of equations (i= 1, 2,..., п) 


dX; n 
"d = — Х; +5 È QikXk + 19) + 10 ex) , (D 
k k=1 
dU : n 
TUN =-U + Хах, (2) 


where the activities {X}; are the proportions of excitatory and inhibitory cells, 
respectively, firing per unit time, /(1) are the input signals to the excitatory and 
inhibitory populations, and the coefficients а; describe the average strength of 
excitatory (а; > 0) and inhibitory (а; < 0) synapses (Wilson and Cowan 1972), the 
sigmoid function S(x) in (1) represents the neuron response to the inputs from the 
adjacent neurons and has a form 

S(X) = TE + exp(— B(X — 9)]. (3) 
where the parameters © and В describe, respectively, the average threshold and its 
dispersion over the whole population of neurons. 

The function U has a property of a control function due io the presence of a 
feedback loop in the system. It connects different subpopulations of neurons 
summarizing their output information. The corresponding, evolution equation (2) may 
be considered as one that describes the change of an average membrane potential of 
terminal’s fibers along which the excitatory and inhibitory influences on the network 
are entered (Degtyarenko 1986; Chinarov and Degtyarenko 1990). Parameter T and 
coefficients f, characterize, respectively, the corresponding membrane time constant 
and average action of the excitatory and inhibitory outputs on the terminals. 

Such presynaptic influence is a standard mechanism for the change of efficiency 
of the synaptic transmission at different levels of nervous system (Bayev and Kostyuk 
1981). Modulation of the polarization level at various terminals by the spinal genera- 
tors has been demonstrated during the rhythmic movements of locomotion and 
scratching when strong dependence of the signal entered to the generator's structures 
upon the intensity of the scratching observed (Degtyarenko 1986). 

The system of Eqs. (1) is a modification of the Wilson-Cowan model (Wilson and 
Cowan 1972; Grossberg 1988) which allows to analyze the dynamics of disconnected 
populations of excitatory and inhibitory elements. It describes cyclically connected 
through a feedback loop heterogeneous and extended in space network consisting of 
homogeneous nets of excitatory and inhibitory neurons where each net may have the 
property of multistability. Corresponding scheme is shown in Fig. 1. 
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Fig. І. Schematic diagram of the network with interconnected excitatory and inhibitory nets controlled by the 
function U. | 


It can be easily shown using an affine transformation of the form Z = СХ, that 
system (1) may be reduced to a standard Wilson-Cowan (for the average neural 
activity) model (Wilson and Cowan 1972) or the Grossberg additive model (Grossberg 
1988) (for i= 1,2, ..., n) of the form 

dy; Ш n 
uro c би D S| Пу +4), (4) 
j=l =0 


in the limiting case when T>0. The subpopulations of excitatory and inhibitory 
neurons will be then interconnected via the elements ry. 

Equations (1)—(2) describe the dynamics of а network consisting of 
homogeneous nets of excitatory and inhibitory neurons where each net may have the 
property of multistability. An example of the phase plane of a system with average 
excitatory (E) and inhibitory (J) activities for the case of two different types of 
attractors (i = 2) is shown in Fig. 2. 


Fig. 2. Phase plane of a system with average excitatory (E) and inhibitory (/) activities. 


Simulations on the base of Eqs. (1)—(2) show that this network possesses 
complex dynamic behaviour with different coexisting attractors such as limit cycles 
and fixed points. Moreover, one can also observe a synchronization of activity 
patterns in ensembles with excitatory neurons only (а; > 0). 

The control function U may be also considered as the function governing the 
processes of recalling and storing of information by attractors of the neural network. 
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The term «information» defines here the quantity that allow to remove (or, at least, to 
reduce) the uncertainty inherent to a System at the instant when it passes from one 
stable regime of functioning to the other, from the less ordered state to the more 
ordered one, or properly, from the less complex level of its organization to the state 
with the higher level of complexity. 

The uncertainty can be implied by the possible instabilities of dynamic system, 
or it might be a direct consequence of the non-linear nature of the system dynamics 
generating various complex patterns of behaviour. The whole activity pattern may 
change them evolving into a stable pattern that corresponds to minima of the energy 
surface. These minima might correspond in turn to remembered recall of information 
stored there previously. As a rule, non-linear multistable systems. possess a high 
sensitivity to control parameters leading to multiplicity of solutions beyond instability 
points and giving rise to bifurcations, and even to chaotic behaviour via the bifurcating 
cascade (due to sensitive dependence of the System's variables on initial conditions). 
Some ground for this standpoint of view can be found from the study of 
synchronization of neuron activity in large networks (Hoppensteadt 1986) and 
convergence properties of self-organized mapping algorithms for the formation of 
topologically correct feature maps (Ritter and Schulten 1988). 


3. Noisy dynamics of neural net 


We are interested here in the case when the values of parameters a;, ЛО), cx are 
such that the system (1) has several steady states. These states may be either focuses 
(nodes) or limits cycles. The steady state may be characterized by the domains (or 
basins) of attraction on the phase plane (with the boundary in a form of unstable limit 
cycle), or may possess the separatrix passing through the saddle point (unstable 
steady state) and dividing different basins of attraction. The. information about 
pattern activity behaviour is contained entirely in these domains. 

We do not discuss here the learning rule in its usual meaning as an algorithm of 
parameter change which adjusts the coupling matrix, so that a minimum of functional, 
measuring the discrepancy between the given fixed point and the desired fixed point 
is achieved. We will only try to find here the conditions of transitions between 
attractors of the network described by Eqs. (D —(3) using the information about their 
structure and to develop an approach for calculation of the corresponding transition 
probabilities. 

Due to the possible large-scale fluctuations of activities in (1) there may occur, 
besides the relaxation to steady states, relatively large rare fluctuations causing the 
transitions between them (Chinarov et al. 1993). We will assume that characteristic 
probabilities W of such transitions are much lesser than the reciprocal relaxation 
times z, to the steady state. As a result of transitions between steady states, some 
stationary distribution over the system states is worked out and dependence of 
statistical characteristics of a system on its parameters becomes single-valued when 
the summary activity of neural network (or, rather its ensemble average) depends 
uniquely on the input signals f; We will have then the following inequality when the 
network is to be used as a memory unit 


t, Stsw~ (1/f) ТА, т< We, (5) 


where Tsw is a switching time and то is a characteristic time between the switchings. 

The size of the minimal characteristic scale of the switching time is 1,(t, < Tsw) 
for arbitrarily fast changing of control parameter f; The times И! define «the safety 
limit» for the storing of information by neural network. Thus, for the reliable 
performance of the switching element in the working range of parameters, the 
following condition must be satisfied 


И<<г 1, (6) 


Practically, this is a condition on the relative weakness of noise. It leads to the 
situation when «activation energy» (this definition is of conditional nature because in 
the non-equilibrium system transitions are not discribed by the standard Arrhenius 
law) of the escape from the steady state is much greater than the noise intensity and 
exponential factor determines, to logarithmic accuracy, the nature of dependence of 
transition probabilities W on system parameters. In this case it has sharp exponential 
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dependence on the control parameters f; However, in the vicinity of bifurcation 
points, when condition (6) is not valid, «the safety limit» for the storing of 
information sharply decreases as the working range of control parameter becomes 
nearer to the bifurcation point. On the other hand, the closer the system approaches 
the bifurcation point, the smaller is the characteristic work of the switching and 
energy expenditure in memory units. 

The existence of above-mentioned alternative specifications of parameters for 
performance of optimal ratio between reliability and memory storage capacity of the 
network makes the problem of determination of the transition probability W as a 
function of system parameters (in particular, of the control parameters /(1)) to be of 
great importance. 


4. Fluctuation transitions between network attractors 
To study the noise influence on the network dynamics, we consider first a case of 


the multiplicative noise &(t) in the right hand side of Eqs. (1) — (2). We can rewrite 
them for this case in the following form 


dX; | 
LX (7) 
dt T КАХ), 
Ч = U(X) + 0), 8) 
t 
where 
N 
K(X) = — X; її акХ + ft) + £49) +U Eon. (9) 
k k= 
N 
L(X) = — U(X) + 2 B. (10) 


Let (f be Gaussian zero mean random process. The probability density 
functional for its realization has the form (Feynman and Hibbs 1965) 
Piso] = ехр[– У f f atat By(t — EEC’) |, 
ik ар 
where function В is defined by relations 


Bul) = (B^ o) а» Во) = fat exp(iot)Ba(t), (12) 


where Ві) is a spectral representation of the correlation function of random process 
E(t) 
«ЕІ — т)» = Юд(т)Ва (т), (13) 
where D is the noise intensity. 
For times т, << t << W`}, quasistationary distribution p(X1, X2), arising in the 
phase space has the form 
> > 
p(X1, X2) = W(X, Xa)» (x = (ХІ, x2), (14) 
where 


, S PRODO (Xt - X) »[&o] 


о) = 


> X(— 
W(X, Xa) = DE DEX)P ED 


The integrals here are taken over those fluctuations which lead the system (1) to 
the point X at a moment t. Boundary conditions are such that for 
t>-—, Х(- о) = X, (X, is a point on the initially filled attractor). 

In the case of the low intensity noise limit (0 << D, the system with 
overwhelming probability approaches one or another stable state, while the transitions 
between steady states do not actually happen during such a time. At larger times 


(15) 
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t(r >> И)! in the small vicinity of one of the stable states the quasistationary 
distribution is formed. For the times т, > > under the large (and, respectively with 
small probability) «outburst» of the noise, the system goes out from the vicinity of this 
stable state and goes into the vicinity of the point of observation X(t), U(X.(f)). The 
problem we deal with consists of the obtaining such a distribution and calculation of 
corresponding escape probabilities. In general case the considered network must be 
viewed as part of a larger system — the network and its environment. This gives the 
network a memory, rendering the activity process ХГ), by itself non-Markovian. 
However, the joint process (X;, U) is Markovian one. 

For the case of white additive noise in Eqs. (7) —(8), the most simple approach 
to consider the transition probability density for the joint stochastic process (X;, U) 
consists of seeking the solution of corresponding Fokker-Planck equation. The 
quasistationary solution for this equation can be found in the form p ~ exp(—R/D), 
where R — R(X;, U) satisfies the non-linear equation (for the purpose to clarify the 
ideas of developed approach we consider below only the case of single compartment 
neural net (; — 1)) ((Chinarov et al. 1993) 


Н(К,, Ки; х,и) = 0, (16) 
where 
H = 2 вові + AB,,R,R, Р В.К + K,R, + ІК, (17) 
OR OR 
К, = 2, К, =, (18) 


The function А(х, и) can be associated with a mechanical action of а two- 
dimensional auxiliary system with the coordinaies x, u, and then Eq. (16) is a 
Hamilton-Jacobi equation. The related Hamiltonian equations for the «coordinates» x, 
u, and «momenta» R,, R, describe the extreme paths of the system. The probability 
of the escape from a stable state / is given, to logarithmic accuracy, by the expression 


ИХ, И) = Wo exp(—R/D), К = R(x, uj), (R (xp, u) = 0), (19) 


where Xss uy, are the values of the activity X and control function И of the initial 
system in the stable state f and at the saddle point s. Function R satisfies the 
conditions: R(x, и) = 0, дК/д хі, = дК/ди|; = 0; the point x» и; (on the phase 
plane) is the global minimum of R. 

A quite general approach to obtain equations for the optimal path of the escape of 
a system driven by weak Gaussian noise is based on the Feynman's idea of the 
interrelation between the paths o the system and the driving noise (Feynman and 
Hibbs 1965). The physical picture is as following: the system fluctuates mainly in the 
small vicinity of a stable state f. The transition probability is defined by such a 
fluctuational trajectory which moves from the vicinity of f to the vicinity of the saddle 
point s. This trajectory represents the most probable suitable path which corresponds 
to the optimal fluctuation of the noise. The equations for this trajectory are 
Hamiltonian equations (16). We shall call it the fluctuational trajectory. 

In the case of low intensity noise, the probabilities of defferent noise realizations 
strongly differ among themselves and main contribution into the continual integral 
gives the small vicinity of the noise optimal trajectory to which a certain trajectory of 
the dynamic system (1) corresponds [5]. The equation for this trajectory may be 
obtained using variational principle for P[£(f) ] 


2, f di&a [а/а = 0. 60) 


Solving (7) with respect to 51) we obtain from (20) (for U = 0, i= 2) 


2 : 2 
à» [S IX; X) — Х,а, > 7 n + X) – Х,а -f| =0. qn 
ik k= = 


34 


Stochastic Modelling of Transitions Between Neural Network Attractors 


This equation in common with appropriately chosen boundary conditions defines 
the optimal trajectory X""'(r) in the phase plane (X,, X5) which starts up on the initial 
attractor and goes on to the saddle point vicinity, when the limit of domains of 
attraction is the separatrix passing through the saddle point. For the case of the white 
noise, the transition probability is determined to logarithmic accuracy by the flux over 
the saddle point vicinity in the phase plane and has the following form 


W; = const : exp E 55] dt Y, eros) ; (22) 
ik 
where 
2 
ЕРИ) = 57! [axlat + xw] - Y aX? — ЛФ. (23) 
k=1 


In order to calculate the fluctuation of the membrane potential, we may consider 
the simple case when the dendritic tree configuration enclosing neurons on the scale 
of order of branching length J, (distance between bifurcations), has а spherical 
symmetry Cy,, where Мо is characteristic constant of a given dendritic greed. This 
symmetry decreases with increasing of its radius and is proportional on the scale of 
diffusion length l; to the ratio /;/</,>. Fluctuations of synaptic activity on the surface 
of such configuration determine then in a self-consistent manner the fluctuation of the 
membrane potentials (i. e., fluctuations of dynamic variables X,, X5). Their statistical 
properties may be calculated knowing the nature of damping of input inside the 
enclosing dendritic configuration and the intensity of fluctuations on the synapses 
placed at the distances of order l4, where they are statistically independent. The ratio 
М = 1,/<i,> isa parameter one is to use for the averaging of membrane potential 


«V,» = «WayX; fi» = (1/N) У И. (24) 
j i 


An approach to the calculation of activation energies R of fluctuational transitions 
is based on numerical analysis of the Hamiltonian equations for an auxiliary system 
(16) —(18) or the corresponding equations for the «coordinates»! x, и and momenta 
R,, R, of this system (Chinarov et al. 1993; Chinarov and Smelanski 1993) 


X —0H/OR,, и = 0H/0R,, В, = — 0H/óx, В, = — дН/ди. (25) 


To find R it is necessary to solve the boundary-value problem for Eqs. (25) with 
X Xy u- ugs at f£ —o and xx, u-u,, at f —o (x, и; and Xs, и, are the 
values of x, и at the metastable state and the saddle point). 

The extreme paths starting from a fixed point f make a single-parameter set. 
This can be easily verified by noting that in the immediate vicinity of a point f 
Eqs. (25) can be linearized in xy = x — хь x32u— ир and the solution of (25) is of 
the form | 


R(x, и) =- i 5 АХ 1x12! > 0. (26) 
2 


i=l, 


The coefficients of the matrix | | Ду | | are given by the equations 
8H Age wd SH. 
mp 2 OX;0.X; : (27) 
f f 


ФН 

АА» + = 

whe mi 3E dn 
where the subscript f means that the derivatives are calculated Юг 
X; x2 = Ry, = Rẹ = 0. The initial conditions for Eqs. (25) can be chosen at finite 
х1, х2 and finite t= fy. For {> —o a path parameterized in this way approaches 
X; =X) = 0 (i.e, хә х, u> up, and all points (ху, x2) already passed can be equally 
chosen as the initial ones for the same path. Hence, the problem of calculating an 
activation energy comes to finding the value of the parameter for which the extreme 

path arrives to a saddle point to a given accuracy. 
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Fig. 3. Typical example of fluctuational trajectory (optimal path) which moves from the vicinity of stable state 
f (focus) to the saddle point S. This is a case of two-compartment network with average excitatory (E) and 
inhibitory (Г) activities. 


In Fig. 3 a typical example of fluctuational trajectory which moves from the 
vicinity of f (focus here) to the vicinity of the saddle point s is shown for the case of 
two-compartment network with average excitatory (E) and inhibitory (J) activities. 
This trajectory represents the most probable suitable path which corresponds to the 
optimal fluctuation of the noise. 


5. Summary and Conciusion 


We have studied here the dynamical behaviour of neural networks with spatially 
homogeneous and heterogeneous structures of synaptic connections between neurons 
arising due to the spatial and temporal character of non-linear interaction between 
elements. As a base model we have used the Wilson-Cowan (Wilson and Cowan 1972) 
type equations for the excitatory and inhibitory activities with some modifications that 
describe the change of the membrane potential of the terminal's fibers along which 
the excitatory influences on the neural network are entered. 

The proposed scheme of structural organization for the neural network may pro- 
vide generation of the stable rhythmic activity. Such network organization possesses 
the temporal and phase characteristics which exist in the real movement activity gene- 
rators and can explain phase dependent switching over their parameters under the 
influence of signals with various sensory modality (Degtyarenko 1986; Chinarov and 
Degtyarenko 1990). 

The important characteristics of non-linear dynamic systems with coexisting 
attractors are probabilities of transitions between them. The algorithm for searching 
of the optimal trajectory іп the phase. space of the network which crosses the saddle 
point and connects two distinct basins of attraction, is proposed here. The optimal 
trajectory starts near the focus or the node at some point (its coordinates are 

.determined by the searching procedure based on the minimum action principle), 
passes in a small vicinity of ihe saddle point and immediately gets into the vicinity of 
the second focus or node. The transition probability is determined by the flux over 
the saddle point vicinity. 

We assume here that fluctuations of dynamic variables X; are determined by 
information sourse &(1) due to fluctuations of the synaptic activity of dendritic tree 
enclosing neurons and its topology. The problem of influence of the dendritic tree 
topology on the network activity dynamics we consider in the framework of the 
approach taking into account the fluctuation of the membrane potential averaged over 
the close volume of neurons. Fluctuations of averaged non-equilibrium potential are 
calculated in the following approximation. Characteristic dendritic scale of the 
«dressed» neuron which define the damping length of fluctuations (diffusion length), 
is much lesser than the mean inter-neuron distance and, therefore, the fluctuations of 
the membrane potential for each neuron are independent. 

The developed approach may be useful for providing the optimal control! in the 
neural networks and may be used in various cases concerning the problems of 
processing information by networks with different structural organization. 
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СТОХАСТИЧНА МОДЕЛЬ ПЕРЕХОДІВ МІЖ АТРАКТОРАМИ НЕЙРОННИХ ЛАНЦЮГІВ 
В. ЧИНАРОВ, Т. ГЕРГЕЛІ, С. СКУРСЬКИЙ 


Досліджуються процеси самоорганізації в нейронних ланцюгах з гомогенною та гетерогенною структу- 
рами синаптичних зв'язків між нейронами, а також динамічна поведінка нейронних ланцюгів, що ке- 
руються зовнішніми шумами, Рівняння типу Вільсона-Кована з деякими модифікаціями застосовується 
до прискорюючих та уповільнюючих функцій. Пропонується алгоритм для дослідження оптимальної 
траєкторії, що пов'язує різні царини тяжіння. Ця траєкторія починається в деякій точці біля 
стаціонарного стану (ii координата визначається методом, що виходить 3 принципу мінімальної дії) 
проходить біля сідлової точки та входить в окіл іншого стабільного стану. Ймовірність переходу визна- 
чається завдяки потоку над точкою сідловини. 


СТОХАСТИЧЕСКАЯ МОДЕЛЬ ПЕРЕХОДОВ МЕЖДУ АТТРАКТОРАМИ НЕЙРОННЬІХ ЦЕПЕЙ 
В. ЧИНАРОВ, Т. ГЕРГЕЛИ, С. СКУРСКИЙ 


Исследуются процессы самоорганизации в нейронных цепях с гомогенной и гетерогенной структурами 
синаптических связей между нейронами, а также динамическое поведение нейронных цепей, которые 
управляются внешними шумами. Уравнение типа Вильсона-Кована с некоторыми модификациями 
применяется к ускоряющей и замедляющей функциям. Предлагается алгоритм для исследования 
оптимальной траектории, которая связывает разные области притяжения. Эта траектория начинается 
в некоторой точке возле стационарного состояния (ее координата определяется методом, который 
основывается на принципе минимального действия), проходит возле точки седловины и входит в 
окрестности другого стационарного состояния. Вероятность перехода определяется потоком над точкой 
седловины. 
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Abstract. The mechanism of living objects supersensivity to low intensity factors is discussed. 
Numerical assessments are made of sensitivity and selectivity to microwaves. Threshold effect of 
low intensity factors action is investigated. 

1. There аге many cases of extremely high sensitivity to low-level factors. in 
some cases the low-level factors may be beneficial. For example homeopathic pills, or 
manual therapy. In some cases ihey may be harmful, such as technogenic environ- 
mental factors. Other cases represent very sensitive communication. For example, 
communication by odorants between some butterflies, or cases of so called extra- 
sensory perception. 

2. It's interesting that high sensitivity is mostly observed for the living objects. 
The living objects are opened systems. That means they exchange substances with the 
external world. 

All living systems, even a single cell, possess a sort of free will. Speaking in 
physical language, they are muitistable. The multistability means ihat svstem can 
have several different internal states while the external conditions remain the same. 
In other words, the external conditions cannot determine uniquely the internal state 
of a living object. It is namely this fact that is treated as a sort of free will. 

3. The effect of microwaves on living objects has been studied experimentally for 
more then 20 years 11,21. During the last 10 years microwaves are used for healing 
[3,4]. The microwaves traditionally belong to physical factors and their effect have to 
be reduced to purely physical mechanisms. Till now such mechanisms has not been 
found. The problem is that microwaves are applied to in a very low doses. At the level 
of a single organic molecule this low-intensive microwaves produce the effect which is 
well jammed by thermal oscillations of the molecules. The common physics usually 
deals with gases, or liquids, or solids, and it has nothing to say in our situation. Our 
opinion is that the answer may be found in the internal complexity of the obiect the 
microwaves are acting upon. 

4. This paper is aimed to explain how the weak microwave radiation may produce 
a pronounced effect in a coherent system comprising a large number of molecules. 
The following model chemical system is con 


(1) 
B+X2C* (2) 
СРС (3) 


Here Х, А, В, С denote molecules of different species C* and C denote active anc 
non-active states of the same molecule C. The A,:B, and C concentrations 
maintained constant. Concentration of X is expected tò be established 
via the chemical reactions not nere. 

This system is open. Indeed, as A, B, С can be produced or consu : 
chemicals have to be removed or added to the system in order to Kc constant 
concentrations. 

This system is multistable. 


It is established mathematically that under fixed 
conditions two different concentrations of X, are possible as its stable internal state: 
[Х] =ж, ог [X]=x,, x, < x, Both x, and x; are self-sustaining due to positive 
and negative cooperativity presence in reaction. (1). Switching between x, and x; 
states may happen due to external influences, or internal fluctuation processes. 
5. If the a big volume, two stable concentrations may coe 
a reactor, occupying different parts of the Serag We consider the sm: 
reactor іп y single stable concentration is possible. This is 


ж Робота ініційована співробітництвом з міжнародною програмою Vidhuk. 
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reactor. The coherent reactor responses to the external influences as a single unit, 
namely, it switches from one stable state to another one in the whole volume at once. 
The number of C-type molecules in the coherent volume is essential for sensitivity: 
the more molecules, the more sensitive the system is. This number is estimated for 
biologically realistic parameters as М ~ 10? + 1012, 

6. It is known that microwaves are able to shift equilibrium in some chemical 
rcactions, but the shift may be very small. This is because molecular vibrations in this 
frequency band is damped by viscous friction. If the mw radiation power surface 
density is I = 1mW/cm? we may expect the concentration increment to be 

x 
AE ant 4) 

ІС") 
As we know, the system (1) — (3) is able to switch from its stable state to another one 
due to fluctuations. This process is characterized by a mean waiting time for the first 
switching, Т1. If an external influence is added, the mean waiting time will 
change. How much it changes depends on the system sensitivity. Mathematical 
treatment under some simplifying assumptions (for details see [5,6 D gives the figure: 


T(without field) ~ exp(10*) 
T(with | mW/cm? field) Ps ^ (5) 


Comparing (4) and (5), we may conclude that singie receptor molecules are very little 
influenced by microwaves. But if a large number of such receptors are incorporated within 
a living-like multistable coherent system, the whole system may react very sharply. 

7. In connection with this conclusion there are two interesting observations: 

The first one is the reductionism problem [7]. If we have here such a sensitive 
system, is it possible to reduce its properties to the properties of its parts. Theoretical 
basis for such reduction has been offered in this paper. But what may be seen in 
experiments? The C-type molecules realizing the primary reception of microwaves 
show extremely small effect, whereas the whole system manifests a very high effect. 
It seems incredible that the effects are in a casual relation. 

The second observation is that the model system we have considered, has at least 
two common features with usual sensory systems, such as olfaction, taste and others 
[8]. Namely, in analogy with sensory systems, this system is hierarchical. The first 
hierarhichal level here is the C-type molecules, and the second one is the level on the 
whole system, including the cooperatively coupled X-molecules. Also, in analogy with 
known sensory systems, we have here a threshold-type behavior: a certain threshold 
must be overcome for switching. 


ОЦІНКА СЕЛЕКТИВНОСТІ TA СЕНСИТИВНОСТІ КООПЕРАТИВНИХ СИСТЕМ ПО 
ВІДНОШЕННЮ ДО НИЗЬКОІНТЕНСИВНОГО МІКРОХВИЛЬОВОГО ВИПРОМІНЮВАННЯ 


О. ВІДИБІДА 


Обговорюється механізм надчутливості живих об'єктів до факторів з низькою інтенсивністю. Наводять- 
ся числові оцінки чутливості та селективності до мікрохвиль. Досліджуються порогові ефекти дії фак- 
торів з низькою інтенсивністю. 


ОЦЕНКА СЕЛЕКТИВНОСТИ И СЕНСИТИВНОСТИ КООПЕРАТИВНЬІХ СИСТЕМ ПО 
ОТНОШЕНИЮ К НИЗКОИНТЕНСИВНОМУ МИКРОВОЛНОВОМУ ИЗЛУЧЕНИЮ 


А. ВИДЫБИДА 


Обсуждается физический механизм сверхчувствительности живых объектов к низкоинтенсивным факто- 
рам. Приводятся численные оценки возможной чувствительности и избирательности биологических систем 
по отношению к микроволнам. Исследуется пороговый эффект воздействия низкоинтенсивных факторов. 
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Abstract. Mathematical model of the immune response to the introduction of the antigen is 
costructed. The model includes both cell-mediated and humoral responses. Use of biochemical 
kinetics makes it possible to include the regulations of conjugated processes and to exploit the 
enzymatic nonlinearities reflecting the cooperative character of biosynthesis and the cell 
growth. 

(Received August 20, 1995) 


The development of the immune response to the virus infection of the organism 
goes through the stages of specifying information about the virus itself in the direction 
of biosynthesis of the cell-killers and antibodies [1—3]. The identification of the 
regimes which show simultaneously all regulatory mechanisms for survival leads to 
constructing the kinetics of conjugated and parallel processes in regulation of the 
immune response. To this end we use the methods developed earlier for the 
biochemical processes [4, 5]. The mathematical model obtained continues the 
Marchuk-Petrov: model to include regulations of the conjugate processes and 
biosynthesis on the basis of the kinetics of biochemical processes. Thus, one succeeds 
to avoid introducing time lags in the equations of the model. The time lags do exist 
and are represented by nonlinearities in the form of rational functions which make 
the dynamic system smooth. The regimes of various level stability as well as the 
unstable ones that show chaotic oscillations are obtained. High sensitivity of the 
kinetics of biosystem depending on the choice of its initial state is shown also. 

The similar approach was developed by Lumb J. R. [6], Kaufman M. and 
Thomas R. [7], De Boer R. J. [8]. They used the nonlinearity in the form of Hill- 
function for the study of the steady-state multiplicity and also investigated the 
kinetics of an immune response at the initial stage. We have adopted some ideas from 
the theory of biochemical processes so as to study a great deal of the regulations 
including conjugated processes and biosynthesis of the immune response components. 
Such an approach allows us to construct regimes of the best survival or homeostasis 
which acts as the adaptor of the biosystem with a complex structure [6, 10—12]. The 
simulation performed in this work demonstrates a high sensitivity of the biosystem to 
both the nature of the initial state and the open conditions of external medium. The 
phase-state space of the biosystem departs into adoptation domains when each of 
them approaches its adaptor. 


1. THE MODEL 


The main dynamic variables of the model will be given by concentrations of cells, 
virus particles, antibodies and macrophages that reflect changes in their number in a 
biosystem and their functional characteristics in the immune system [1—4, 11, 12]. 
All participating cells and particles can have different localization in separate organs 
and this will be partially reflected in the model disregarding the homogeneous de- 
pendence of all dynamic variables on the space coordinate. Such an interpretation of 
real biosystems is based on high specificity of the recognition of the target cells by the 
immune cells and cooperative character of catalytic processes in the living matter. 
Although we do not distinguish enzymes of the immune pathways, their participation 
in biosynthesis of antibodies, in functioning of receptors, in reproduction and 
elimination of cells makes it possible to apply the kinetics of the biochemical processes 
[4, 5] to the discription of the immune homeostasis at the level of a chosen biosystem. 
As a result, the immune system is represented by the two types of the immune 
response-humoral one and cell-mediated one. The immune system adjusts tunes to 
killing the virus particles circulating in blood and the virus infected cells. The virus 
particles ¢ and the infected cells C are thus separated in space. Viral particles infect 
the cells and the infected cells reproduce the virus particles in the blood. This way the 
virus infection of the organism gains its positive regulation for self-survival under the 
open conditions of the biosystem. Since a separate organ is meant we shall assume 
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that concentrations of the virus infected cells C and target cells c are in the depend- 
ence preserving the integral of motion 


C(t) + c(t) = const = Ly. (D 


The constant L; has the concentration dimension and reflects the total concentration 
of cells that undergo infection. 

The macrofages M prove the beginning of antigen recognition. Their total 
concentration L, in the blood of the organism is assumed to be permanent. They аге 
the first cells recognizing the antigens and they are thus believed to be in the 
biosystem as free cells M and also as complexes m with virus particles 


M(t) + т(0 = const = Ly. (2) 


These two conservation laws allow one to decrease the dimension of the system of 
model equations up to 10. The macrophages participate both in humoral and in cell 
immune responses. Together with virus particles they make up those two components 
of the model through which there proceeds the interaction of the two subsystems — 
humoral one and cell one [2, 3]. Here we do not compare the model to experiment 
[11, 12]. : 


2. THE KINETIC MODEL. THE SYSTEM OF IMMUNE RESPONSE 


The system of immune response of the organism to virus infection includes both 
the mecanisms of virus recognition and specific reaction to the destruction. We shall 
describe the kinetics of virus represented by particles freely circulating in blood. 
Localization of the virus particles will be in the form of complexes of the virus 
particles and immune-competent cells and also the virus infected cells. The 
conditionality of the complexes formation will be given by dynamic variable $(1) 
determining their concentration or their number, and also by concentrations of the 
other participating particles M(t), C(t), E(t), F(t). Under the complex 
Мф or C-ó, Е-ф, Еф we shall mean a generalized act of the interaction, 
integrating stages of the biochemical processes to recognize the infected particles. 
Application of the theory of biochemical processes and indirect participation of 
enzymes at different stages result in the appearance of nonlinearities which determine 
the terms of concentrations of the probabilities to create the corresponding complexes. 
These are: 


и = M) =P —, Vn с) 0, WOVE, Vt E 
1. ф + Ls’ 2 ф+ 13 ’ , $-L 


here we used the notation V(Z) = Z/(1 + Z). Restricting ourselves by 10 dynamic 
variables we try to make up the model apprehension of regulations and processes 
through the equations and nonlinearities in terms of these variables only. This 
assertion acquires more complete meaning when we consider the equations for 
changing the rate dv(1)/dt using the kinetic laws 


4 = nl W(CW(E) - nii ИЕ) sta — m Ys - С) = e raa 


g? Ms (3) 
— nay3V(L, — М) ZIL +m $* M; + паВ.С — аф. 


where, for the sake of generality, there is a term (parameter m) reflecting permanent 
introduction of virus in the organizm as well as its permanent outflux from the blood 
system (parameter a). The numbers т, n, пі, М2, пз, пд, а represent а part of the 
contributions coming from various sources and expenditure of virus particles in the 
immune response of the organism. The virus particles can reproduce themselves in the 
infected cells forming concentrations of virus particles in the blood (parameters n and 
nj. The parameter l, is responsible for the destruction of infected cells C by 
T-limphocytes-effectors E. The second term (parameter f) determines the appearan- 
ce of virus particles in the cells lysis process. The growth of concentration of virus 
perticles is suppressed in two ways. The $ expenditure for cells infection L, — C 
(parameter n; does not mean the beginning of the immune response yet. The 
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beginning proves to be the formation of complexes М - ф (parameter пз) with rate of 
the adsorbtion probability on the macrophage surface L4 — M of the virus ¢ in the 
form | a Й 


TICE — M) E 
y3V(Li М) Db. 

The second degree ф' reflects cooperativity and, as its result, the retardation in 
generating the primary signal about the appearance of virus particles in the organism. 
The degree and the parameter Ls give an integral characteristic of this retardation at 
the initial stage of the immune response..In the given case n, means the. number of 
virus particles necessary to generate the primary signal as a result of the interaction 
with macrophage. Such an interpretation operates also with certain notions and has 
the model meaning only. In this sense ihe parameter n reflects the growth of ihe 
number of free virus particles due to the outtlux from the infected cells where their 
infection occurred. Some virus particles are spent for penetration into healthy cells 
L, — C (parameter n). A humoral way is exhibited by antibodies F participation in 
ihe neutralization of free virus particles (parameter n,). This balance of reproduction 
and output of virus particles from the organism should also be supplemented with 
initial concentration of virus which penetrated the organism ¢(t=0)=¢". The 
mechanism used to stimuiate macrophages to transmit the information has enzymatic 
character when the signal from virus is received by acceptors and is realized by the 
enzymatic system. Thus the number of stimulated macrophages under conservation of 
the general number of macrophages increase according to the law (2). Stimulated 
macrophages for double recognition, i. e., for the formation of complex with 
T-lymphocyte ог B-lvmphocyte, are described by the kinetic equation 


‚2 


ЎА, 4,2 
aM. = У; V(L, — M) EE ii азМ, (4) 
ae d й і Q^ + Ls 
supplemented with the initial data taking into account the conservation law (2), 
M(t 2 0) = М, M? x Lı. The removal of stimulated macrophages from the organism 
leads to the corresponding reproduction of free macrophages according to the law (2). 
Receiving signal by T-lymphocytes for reproduction H, and B-lymphocytes for 
production A, has also cooperative character which leads to expand time lag. Our 
model interpretation is represented by the following kinetic equations 


LS C EVODVOM с на QS REP А 
dE 2 "9 А Н+ М, Н+ Ку. 3 (УМН, — ан, (5) 

dH; — M^; H1 > 

—— = kV XN VM чор V(X) VM) VAD) – оН), " 


where the regulating growth of Hy and H, is taken into account by the levei of © 
stimulated macrophages, and hence by the level of virus particles conceniration in the 
organism. This is represented by the parameters M,, Kı and М», К. 

The rate of the growth Н, is defined by the initial concentration of X? and cell 
division by the law У(Х) with the rate coefficient kı. All other factors in this term are 
responsible for the regulation of growth by stimulated macrophages М and their 
interaction with cells Z;. The retardation of the cell division H, is regulated пої 
only by macrophages М but also by their interaction determined by the factor 
Hi/(Hi + K,). Together with another factor M,/(H{ + Му) they define the level of 
the minimum concentration being the starting point for the H, growth and the 
maximum level of H, being the starting point for the inhibition of the cell growth Hj. 
Together with dissipations- (parameter c) all this results in stabilization of the level 
of Н, cells. Other fall of rate growth H, is due to their (parameter Аз) participation in 
the division of E cells, as defined by the equation | 

a = lk; V(X)V(M)V(H1) - із V(C)V(E) - «Е, (7) 
the E cell growth is:stimulated by the double signal from the macrophages M and 
from the H;-helpers. The additional factors regulating growth are not included here 
except for their expenditure for elimination of the infected cells C (parameter /,) and 
their permanent output from the blood system. The delay in the formation of the cell 
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immune response can be revealed only in the numerical experiment for the fixed 
values of the parameters. The cell homeostasis acquires a specific meaning in the 
blood biosystem represented in the model by regulations of the conjugate processes. 

The kinetics of the virus infected cells is regulated by their level (1) in the 
organism and the evolution of the desease (parameter 8,) with the following immune 
response (parameter /3) as shown іп the kinetic equation 


Te = Виа ~ O BF ВИСО) - фі + aC. (8) 


The infected cells are removed due to lysis (parameter В.) and also due to the 
output from the organism (parameter a4). The constant L, regulates both the cell 
infection by virus and it elimination by T-lymphocytes and this results in partial 
infection of the functioning organ. The infected cells of defferent origin are 
conventionally taken here into account 


The growth of this biomass is also an additional signal for generation of the 
immune response. 

Formation of a humoral immune response starts with proliferation of T-lympho- 
cytes Н, according to Eq. (6). This signal together with macrophages and B-cells 
stimulates the P-cells into the synthesis of antibodies P as reflected in the kinetic 
equations 


GF = Bok VOD VO) VUE) VO?) — ав. (10) 
А З і і F + М 7E, 
ЧЕ E. E 
deo!) pee Э етер леа (12) 


where the second term (parameter y,) means that free virus particles ф are eliminated 
by the antibodies F. 

The system of kinetic equations should be supplemented with the initial data 
taking into account the integrals of motion (1) and (2) 


p(t =0)=¢°; M(t=0) = М; H,(t=0) є НЗ; H(t = 0) = НЗ; E(t = 0) = E? 
P(t = б) = Р; Bt = 0) = В; F(t = 0) =F; C(t = 0) = 73 X(t = 0) = AG (13) 


where C? < L,, M? < L, and all introduced parameters ф0, МО, ..., X? are nonne- 
gative and finite. 


3. COMPUTER SIMULATION 


The model of the arising immune response to introducing virus in the organism 
and its futher evolution is fixed by the two groups of parameters: the internal 
parameters characterizing the autoregulation of the conjugate processes in the 
biosystem ; 

Li = 3; L = 1.3; Із = 0.7; Ly = 0.5; L; = 0.5; М, = 0.04; 
М) = 0.3; М, = 0.05; Mg = 0.05; К, = 1.2; К, «0.3; І- 1; 4 = 0.01; 
паб; п = 1; п. «0.1; пз. 20.1; и = 1; 
В; = 0.0031; В, = 0.013; у; = 0.006; уз = 0.01; 
К, = 1.5; ky = 0.2; К, = 0.01; К, = 0.07; (14) 
Po = 0.26; Е = 0.1; By = 0.027; 


and the external ones 
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т = 0.001; а = 0.00012; ол = 0.001; a, = 0.0001; аз = 0.002; 
24 = 0.0001; аз = 0.0002; а, = 0.0002; 
a5 = 0,0001: аз = 0.0001; ао = 0.005; 
із of the Diosysiem. These parameter 
al state of the biosystem 


BEN 


that determine the regime of homeost 
necessarily be suppiemented with t 
t= 0 when the latter is being virus infected 
фі = 1.56; МО = 2.05; Н = 2.16; НЗ = 2.11; FO = 2.28; 
В? = 0.58; P? = 0.11; Р = 0.26; С? = 0,51; X? = 1.61. (16) 


This terminology as well as the initial state of the biosystem have not appeared 
unexpectedly, their meaning has arisen on the basis of the regime selec in the 
appropriate way on the phase portrait of dynamic system determined above by. the 
fixed parameter values. Among various regimes of the biosystem of immune response 
to the evolution of virus infection we single out the conditions of the external medium 
io which the organism is to be adapted under conditions of the continuous ingress into 
organism of virus particles in limited concentrations (the parameter m) to the regime 
of the best survivai with creating homeostasis of conjugated processes in the biosystem 
revealing the immune response. The existence of such regimes underlies choosing the 
appropriate system of the parameter values and the initial state. 

The availability of the double immune response generates high sensitivity to the 
initial infection by virus. Рог ihe values $? < 0.135 and ф0 > 0.14 ihe immune 
response is distinguished as: for ФО < 0.135 the cell immune response (see Figs Та) is 
determined by high ievel of tymphocyte-effectors E, by their interaction with infected 
cells C. In spite of the fact that no virus cells ф growth is observed (see Fig. 1b) their 
stationary level testifies 1 ibsence of humoral immune response. This is true 
during fixage of the remaining components of the biosystem in the form of (16), 


0.80 5 


^ — £00090.00 


stem on the initial concentration of virus сей: The depend- 
2; 0.125; 0.13; 3; 10, corresponding here and in the following 


endence of the kinetics of virus cells in free form. 


ence of the kine! 
figures to num 
ф = 0.13; 0.135; 0.14; 0.15 


р" > 0.14 one observes the induction of antibodies biosynthesis 
(see Figs. la, b). Simultaneous appearance of both systems of the immune response 
makes it possible for ihe organism to overcome high initial concentrations of virus 
particles preserving homeostasis of the conjugate processes in the biosystem, The 
continuous influx of virus parti tasis biosystem of 
ihe immune responsc. 


For the values $" 


cles (parameter m) results in homeos 


4. PARTICIPATION OF MACROPHAGES IN THE FORMATION OF SIGNAL OF 
VIRUS INFECTION ` 


The dependence of ihe kinetics of macrophages participation in organizing the 
immune response is determined by many factors among which one distinguishes the 
regulations influencing the level of their general concentration Z4 in the organism. 
The numerical experiment reveals high degree of the infection of organism at 
comperatively low concentrations L, (see Fig. За), Іл < 2,7. Along with this, at 
increased concentrations, L, > 3.65, the same picture is observed. The reason for the 
weak immune response (curve 1) is the weak signal from macrophages whose 
concentration is lowered. In the second case under the increased concentration of Ly 
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the immune response is determined by the fact that macrophages occupy a new 
stationary level. The mechanisms of the appearance of stationaries at І = 2.6, 2.65 
and Г. = 3.7 with an intermediate region of nonstationary regimes among which 
there are also the nonstable ones can be identified in Fig. 2b. At L, = 3.7 the cell 
immune respon Se is guar ranted ee ise in iE 2 RE е и. ds шова! 


in @ " var’ 
kineti i із closely 
stilts given in this figure 
reflects simultaneous 
15 in our conventional 


e ed in M 


poe em, "he app peara 
asso d with the fall of humoral imm 
show that homeostasis formation 
func tioning of both the cell and hun 
frame ої reference (14) ) 
levels and in the follo 


є response 
the immune 
i 1e тезро 
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5. 2. The influence of macropnag ges on the kinetics of biosystem: a) Arising of stationary regimes (curves |, 
2 2 and 5) in @ at L1 = 2.6; 2.65; 2.7; 3.65; 3.7; b) General picture of the biosysiem at L1 = 3.7. 


5. PARTICIPATION OF INFECTED CELLS IN ORGANIZING THE IMMUNE 
RESPONSE 


cyte iubes ihesis Re in pra ce outpi ut of virus Cte from ris organism bur in such 
a way that the following law is satisfied (1): 


n of th 
? means the 


Biosynthesis is inducted both by C ceil elimination and reproducti 
concentration up to the level І». The C cell eliminati on by the killers 
reproduction of their population in the form of C according 1o the kinetic e uation (8). 
The cell concentration (1) promotes both the lowering of virus cells and the 
reproduction of the level ф in a stationary regime. Here for £5; = 1.297 and L; = 4.0 
the stationary regimes are observed. The fig thie: bi ication bin 
L, € [1.297; 1.298 |, that shows how the oscil 


ady 
state one, is called the Hove найом з ‚ Неге ‘the 
parameter Г. (see Fig. 3b) а : second 
steady state regime for L, = F Te the | ter rmines ‘the period of the 
evolution. retardation (see Fig. 3b). T high dependence і eryed ` for 
І, Є [3.5; 3.66]. The point L,.& [3.62; 3.64] determines the second point "Hopf 


bifurcation creation of oscillations from steady state regime when i, changes from 
Ір = 3.64 in the direction of the point ZL, = 3.62, The region of complex oscillational 
regimes including the unstable ones, is situated in the range between these values. 
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Fig. 3. The kinetics of virus infected cells: а) Participation of virus particles іп the formation of the immune 
response for Lz = 1.297; 1.298; 1.3; 3.0; 4.0; b) Formation of steady state regimes at increased 
concentrations: [2 = 3.5; 3.6; 3.62; 3.64; 3.66; c) The role of B-cells in the formation of kinetics of H2 at 
12 = 1.25; 1.3; 1.7; 2.5; 3.5; d) Dependence of the kinetics of plasma cells in the formation of the immune 
response for Lz = 1.297; 1.298; 1.3; 1.5; 1.7; e) The cell immune response of E and gradual onset of the 
humoral response effect on the whole process with the increase in L2 = 1.25; 1.3; 1.7; 2.5; 3.5; f) Appearance 
of antibodies Е coincides with bifurcation according to the parameter L2 = 1.25; 1.3; 1.7; 2.5; 3.5. 


The levels of general concentration of the cells L; can influence the formation of 
the regimes in different components. The figure Зс shows how the low level L; = 1.25 
results in the appearance of steady state regime at high level H;. As a result this 
regime leads to the decrease in cells P and simultaneously to the decrease in 
antibodies (see Fig. 3d, curve 1). The Hopf bifurcation on the interval 
L, Є |1.297; 1.298] is also shown here. Various regimes according to their asymptotic 
stability are represented in Figs 3c and 3d. 

These complex regimes are reflected both in the cell (see Fig. 4e) and in humoral 
responses of the biosystem (see Fig. 4f). The absence of oscillations on curves 1 in 
both figures evidences that fall of antibodies F results in a stationary with high 
level E. 

It is shown in curves 2 that the appearance of antibodies F leads to the 
interaction between the cell and humoral subsystems characterized by complicated 
oscillations in the antiphase. The oscillation region has stable vibrations (curve 4) and 
the asymptotic output onto the limiting point (curves 3). Curve 5 leads again to 
complex oscillations. The structure of the oscillation will be studied in the next paper. 
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6. DEPENDENCE OF THE IMMUNE RESPONSE ON VIRUS INFECTION OF 
ORGANISM 


In the case when the virus infection is constantly supported by the medium the 
external regulation of the immune response takes place. The infection of organism 
depends not only on one-fold introduction of virus cells but also on their constant 
supply to the organism. This leads to adaption of the biosystem to the regime of best 
survival under the open conditions of the input (parameter m) and output (parameter 
a) of virus cells. A continuous background of virus infection with participation of the 
medium is maintained by these two parameters. 


Fig. 4. Influence of the immune background of the 
medium (parameter т) on the kinetics of biosystem; 
m = 0.00095; 0.00096; 0.0015; 0.0023; 0.0030; B-cells 
determine the discontinuity of biosynthesis of 


antibodies in connection with the presence of rigid T te ~200060.00 


bifurcation of a stationary where F = 0; t 


The external parameters can be chosen according to our model interpretation of 
virus infection, and the numerical experiment can predict the behaviour of biosystem 
in response to the levels of virus cell concentration in the organism. Various immune 
system responses to the external behaviour of virus cells are shown in Fig. 4. The 
kinctics of the biosysiem has two stationary regimes (curves 1 and 5). There are two 
Hopf bifurcations on the intervals т Є (0.00095; 0.000965 1 апа те (0.0023; 0.003 | 
at с = 0.00012, resulting in stable vibrational regimes (curve 2 and 4). The reason for 
arising of the steady state regime at m = 0.00095 is so low level of B-lymphocytes 
that there occurs the fall of humoral response, complete absence of antibodies F. The 
kinetics of self-organization of conjugate processes for the formation of humoral 
response is shown in Fig. 4 as a sequence of organizing the humoral immune 
response. When m - 0.00095 the system is eliminated. This kinetic dependence of 
F-lymphocytes is CORREIA with the kinetics of macrophages M, B-lymphocytes are 
returded to such an extent that they are unable to reproduce their population. 

At the same time for m = 0.00096 the biosystem organizes the immune response 
the humoral one including (see Fig. 4) that allows a comparatively low level of 
concentration of the virus cells. 


CONCLUSION. 


performed has revealed the existence of 
the phase space regions where the cell- 
mediated and humoral ways of the 
immune response are functioning 
simultaneously. At the same time the 
fall of B-cells below certain level can 
lead to regimes with the fall of humoral 
response. Possibilities of the model are | 2798. 
presented in Fig. 5 where possible Fig. 5. 

regimes with B-cells (parameter 

ав = 0.00005; 0.0001; 0.00015; 0.0002; 0.0003) are illustrated. The results obtained 


allow one to carry out further construction of the immune response model. 
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МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ КЛІТИННОГО ТА ГУМОРАЛЬНОГО ПРОТИВІРУСНОГО 
ІМУННОГО ВІДГУКУ. 1. РЕЖИМИ, ЩО САМООРГАНІЗУЮТЬСЯ 


В. П. ГАЧОК, М. П. ЧЕРНЯК 


В роботі побудовано математичну модель імунної відповіді на ураження вірусом. До моделі включено 
як клітинний, так і гуморальний відгук. Завдяки використанню біохімічної кінетики виникає мож- 
ливість застосувати регуляції спряженими процесами та використати нелінійності, що відображають як 
кооперативний характер біосинтезу так і ріст клітин. 


МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ КЛЕТОЧНОГО И ГУМОРАЛЬНОГО 
АНТИВИРУСНОГО ИММУННОГО ОТВЕТА. САМООРГАНИЗУЮЩИЕ РЕЖИМЫ 


В. П. ГАЧОК, Н. П. ЧЕРНЯК 


В работе построена математическая модель иммунного ответа на вирусное поражение. Модель включа- 
ет как клеточный, так и гуморальный ответ. Применение биохимической кинетики дает возможность 
включить регуляции сопряженными процессами и использовать нелинейности, отражающие коопера- 
тивный характер биосинтеза и клеточного роста. 
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Сравнительньй качественньй анализ эволюции системы сталкивающихся электрона и протона 


СРАВНИТЕЛЬНЫЙ КАЧЕСТВЕННЫЙ АНАЛИЗ ЭВОЛЮЦИИ 
СИСТЕМЫ СТАЛКИВАЮЩИХСЯ ЭЛЕКТРОНА И ПРОТОНА 
С УЧАСТИЕМ ЭЛЕКТРОМАГНИТНЫХ ИЗЛУЧЕНИЙ 
В КЛАССИЧЕСКОМ И КВАНТОВОМ ПОДХОДАХ 


В. С. ОЛЬХОВСКИЙ, С. П. СИТЬКО 


Й 


COMPARATIVE QUALITATIVE ANALYSIS OF EVOLUTION OF A SYSTEM OF COLLIDING 
ELECTRON AND PROTON WITH PARTICIPATION OF ELECTROMAGNETIC RADIATIONS IN 
THE CLASSICAL AND QUANTUM APPROACHES 


V. S. OLKHOVSKY, S. P. SIT'KO 


Summary: Evolution of a system of the colliding electron and proton is analyzed qualitatively 

and conceptually taking into account the bremstrahlung and radiation (PHOTON) transitions 

in framework of the classical and quantum approaches. We have found theoretically, using the 

quantum approach, processes of multiple formating excited and high-excited (Rydberg) states 

of a hydrogen atom with gradually increasing excitation energies and dimensions with time. 

Such processes are a direct consequence of processes of internal excitations and ionizations of 

the formated atom by a photon emitted in preceding radiation captures of the electron by the 

proton to underlying atomic levels. Chains of all these processes compete with direct radiation - 
decays of all initially excited levels and bring to retarding of the final formation of the ground 

state of a hydrogen atom. ` 


ОТ РЕДАКЦИИ 


Проблема образования атома водорода (а тем более других атомов) с набором 
стабильных энергетических состояний до сих пор не обсуждалась в научной лите- 
ратуре. Решение этой проблемы открывает путь для установления и понимания 
тех процессов, которые следует отнести к первоначалам вещей и которые не 
гипотетически, а научно обоснованно определяют закономерности возникновения 
химических элементов, лежащих в основе образования всех имеющихся в приро- 
де веществ. Естественно, начинать решение указанной проблемы следовало с 
атома водорода, причем самого простого, протия. Это и было реализовано автора- 
ми рецензируемой статьи, которая, на мой взгляд, представляет несомненный 
научный интерес. 

Академик НАН Украины М. П. Лисица. 


1. Введение. Хорошо известно, что классическая физика не в состоянии опи- 
сать устойчивый атом. В квантовой же физике давно разработаны методы расчета 
вероятностей радиационного захвата протоном движущегося электрона в любые 
состояния атома водорода, из которых только основное состояние является устой- 
чивым, а любое возбужденное состояние является относительно метастабильным 
и испытывает радиационный распад на нижележащие (по энергии) состояния. 
Однако в полном объеме эволюция столкновения электрона и протона с формиро- 
ванием суперпозиции всех возможных состояний атома водорода и всеми воз- 
можными процессами испускания и поглощения фотонов до сих пор явно не изу- 
чалась даже в принципе. Для изучения синергетических процессов и явлений 
самоорганизации в физике живого временной ход (эволюция) образования устой- 
чивых и метастабильных состояний ядер и электронов, начиная уже с простейшей 
системы одного протона и одного электрона, имеет принципиальное значение. 

Настоящая работа посвящена качественному и концептуальному анализу за- 
дачи поведения во времени движущихся (в системе центра инерции) электрона и 
протона в рамках классического и квантового подходов. Для простоты и выявле- 
ния существа дела ограничимся: простейшими приближениями, оговариваемыми 
по мере необходимости ниже. 

2. Классический подход. В рамках классической теории, если не учитывает- 
ся тормозное излучение электромагнитных волн, изучение движения электрона в 
поле притяжения протона сводится в системе центра инерции к изучению движе- 
7 m. : m, 
m, + m, 


ния одной частиць c приведенной массой А = ; практически совпадающей 


49 


B. С. ОЛЬХОВСКИЙ, C. П. СИТЬКО 


с массой злектрона т, (m, — Macca протона), в кулоновском поле, описываемом 
формулой 


е? 
U(r) = — (1) 
r 

в точечном приближении и без учета магнитного взаимодействия (e — заряд 
электрона, r — расстояние между протоном и электроном). Решение уравнения 
движения UT = —er/r? дает траекторию, лежащую в одной плоскости и описывае- 

мую геометрическим уравнением [1] 
r=M?/u- е?(1 + е : cosp), (2) 


где ги рф — полярные координаты частицы, € — эксцентриситет орбиты, равный 
УІ + 2ЕМ?/и е, E = Ex, + О(г) — полная энергия частицы (Egg = и · v?/2, 
у — скорость относительного движения протона и электрона), М — момент коли- 
чества движения частицы относительно центра поля. При движении в поле (1) 
сохраняются Е и М. 

При М - 0 (лобовое столкновение) г -> 0, т. е. происходит падение электрона 
на протон. При М > 0 и Е > 0 движение инфинитно (по гиперболе при Е > 0 и по 
параболе при Е = 0, т. е. когда частица начинает свое движение из состояния по- 
коя на бесконечности). При М > 0 и Е < 0 частица движется по эллипсу с осями 
а=е?/2 · IEI и bz М/У2и · IEI. Период движения T по эллиптической орбите 
равен Т = л · e?VM/2IEI?. При наименьшем допустимом значении энергии TAKO- 
го движения, равном Big = —uet/2 М?, эллипс превращается в окружность 
(є = 0). Учет тормозного излучения электромагнитных волн делает невозможным 
устойчивое движение по эллипсу; при этом происходит падение электрона на 
протон по спирали. Время падения по спирали при начальной скорости у, направ- 
ленной перпендикулярно начальному расстоянию а между электроном и прото- 
ном, определяется формулой 
wea? 

4e4 ' (3) 
Для иллюстрации и в дальнейшем использовании по принципу соответствия 


в квантовом подходе приведем примеры численных значений времени падения 
электрона на протон при М = 0 (без учета тормозного излучения) по формуле [1] 


ta = 


UA 


RA , (4) 
BS io. 
где а — начальное расстояние, с которого падает электрон с энергией 


Е = Eua — е/а?. При Ехинна- = 0,1 эв и а = 1 см получаем t; = 0,6`. 1077 сек; 
при Екиннач = 0,128 и а = 100 А получаем to = 0,6 · 10-13 сек. В случае движения 
по спирали время падения электрона, определяемое по формуле (3) са = 100 А, 
оказывается равным t, = 2,7 · 1077 сек, т. e. электрон будет падать сравнительно 
очень медленно. 

3. Квантовый подход. В квантовой теории (см., например, [2, 3p электрон и 
протон, рассматриваемые в системе их центра инерции в вакууме изолированно 
от всех других частиц и взаимодействующие друг с другом только посредством 
кулоновского поля (1), могут находиться либо в состояниях, отвечающих диск- 
ретному спектру отрицательных значений полной энергии, образуя атом водоро- 
да, либо в состояниях, отвечающих непрерывному спектру положительных значе- 
ний полной энергии при инфинитном движении их друг относительно друга. При 
более детальном рассмотрении учитывается метастабильность всех возбужденных 
состояний атома водорода с конечными временами жизни 


в = В/Г,, 220 (п=1,2,....) (5) 
(где Г, — радиационная ширина n-ro уровня энергии E,, определяемая полной 
вероятностью радиационного распада на нижележащие уровни с испусканием фо- 


тонов) и излучение фотонов тормозного электромагнитного излучения при дви- 
жении электронов в непрерывном спектре энергий, т. е. вводится в рассмотрение 
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влияние электромагнитного вакуума, обусловливающее рождение фотонов тор- 
мозного излучения и в радиационных распадах уровней. 

Отметим выполнение принципа соответствия классической и квантовой тео- 
рии в описании свойств системы протона и электрона. При E = 0 выполняется 
обычный принцип Эренфеста (см., например, [2, 3D и относительное движение 
электрона и протона инфинитно и в классической и в квантовой теории. При 
Е < 0 как в классической, так и в квантовой теории электрон и протон располага- 
ются и движутся друг относительно друга в конечной области вокруг г = 0. И 
только при Е > 07 эта область расширяется до бесконечности (движение по пара- 
боле или прямой в классической физике и увеличение среднеквадратичных ради- 
усов возбужденных связанных состояний. атома водорода до бесконечности в 
хвантовой теории). 

Если начальное состояние описывается волновой функцией (волновым паке- 
том, нормированным на 1) WVe(r; 19)! 0,7, где R — средний волновой вектор OTHO- 
сительного движения электрона и протона, tp — начальный момент времени 
1,10,> — волновая функция вакуумного состояния злектромагнитного излучения 
(в представлении вторичного квантования), то с течением времени возникают не- 
нулевые вероятности рассеяния электрона протоном (с сохранением кинетической 
энергии Ey, = 472/2 в асимптотически удаленной области при пренебрежении MO- 
терями на тормозное излучение) и радиационного захвата электрона протоном с 
образованием суперпозиции различных состояний атома водорода и разлетаю- 
щихся фотонов i 


APO BUS UB SS (6) 


где C,(t) — амплитуда вероятности перехода системы из Фет to) 1O,> в cocros- 
ние а! Ру», v.(r) — волновая функция п-го связанного состояния атома водо- 
рода, | P,> — волновая функция испущенного фотона (для простоть ограничимся 
обычным приближением однофотонного излучения в каждом акте излучения) CO 
средней энергией P, + c= E,— E, Величины C,(D определяются матричными 
элементами 


«раб P, а) 15 t), Оу», к (7) 
удовлетворяют условию 21C,(t)! 2 « | и при достаточно больших t (однако мень- 


ших, чем т,) приближаются к предельным значениям С? · e i(EatPye)/b c некото- 
рими константами С°, Поскольку величины матричных элементов (7) существен- 
но определяются перекрытиём волновых функций war, thu w(t), а среднеквад- 
ратичный радиус состояния v.) пропорционален (п + 1)? [2, 4], то в Х в (6) 


вносят основной вклад члены, отвечающие значениям п = 0, 1,2,... и T. д. вплоть 
до значений, отвечающих тем состояниям, среднеквадратичный радиус которых 
порядка a. Различные методы расчета амплитуд вероятностей С.Й) изложены, 
например, в [2, 3, 5]. 

До сих пор излагалась в принципе известная элементарная схема описания 
столкновения электрона и протона (с участием электоомагнитного излучения) 
При этом обычно ограничивались фрагментарным и BO всяком случае’ неполным 
рассмотрением эволюции даже такой простейшей системы в пределах одного-двух 
и изредка нескольких переходов OT радиационного захвата до радиационного сня- 
тия возбуждения образовавшегося атома. Залавшись целью выяснения временного 
хода всех возможных процессов, происходящих с атомом водорода и злектромаг- 
нитным излучением, образовавшимися в результате столкновения электрона C 
протоном, рассмотрим, какие типы процессов, какой длительности и в какой по- 
следовательности происходят.с таким атомом вплоть до ето конечного состояния, 

Пусть вначале электрон (электронный волновой пакет» движется навстречу 
протону (протонному волновому пакету), начиная со среднего макроскопического 
расстояния (скажем; 1 см) co средней кинетической энергией относительного дви- 
жения, близкой к 0 (скажем, 0,1 эв). Для качественной оценки средней длитель- 
ности движения на пути 1 см вплоть до совмещения центров обоих волновых 
пакетов можно в соответствии с растиренным принципом Эренфеста [6 | восполь- 
зоваться классической оценкой (4), согласно которой, если сближение происходит 
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с расстояния l см, to =0,6- 1077 сек, а если с расстояния 100 A, то 
to = 0,6 - 107 сек. Известно, uro в процессе сближения электрона с протоном 
имеет место рассеяние под разными углами, сопровождающееся тормозным излу- 
чением малой интенсивности. Для нас наиболее важным в эволюции является то 
обстоятельство, что в процессе сближения электрона и протона происходят отсут- 
ствующие в классической физике процессы радиационного захвата электрона в 
атомные состояния фа. Процесс радиационного захвата в каждое атомное со- 
стояние, вероятность которого определяется соответствующим матричным эле- 
ментом (7), сам по себе происходит довольно быстро, в основном при прохожде- 
нии электроном пространственной области, занимаемой тем же атомным элект- 
ронным состоянием, т. е. длительность захвата можно качественно оценить как 


; = . 10-8 
= м = О сек = (0,5 — 5) - 1075 сек, 


где Ar — линейные размеры (толщина) области, занимаемой атомным электрон- 
ным состоянием (~1 — 10 A). 

Наряду со всем этим, идут также процессы радиационного распада возбуж- 
денных состояний атома водорода. Радиационные ширины уровней можно грубо 
оценить с помощью соотношения [5] 


Тзахв 


(8) 


3 ие 
ree (2) | E,, Е, = рі n2,3,.. (9) 


T. е. Г, < 02-107 /n? эв и т, = В/Г, > 0,3:107!! n? сек. При n > 10 (когда радиус 
такого возбужденного состояния > 100 A) Г, < 0,2-10 5535 и т, > 0,3- 10^? сек. 
При п > 10°, когда радиус такого высоковозбужденного (ридберговского) атома 
> 10? см, Г, < 0,2 - 1079 эвит, > 0,3 · 10-5 сек; при n > 10^, когда радиус pun- 
берговского атома > 1 см, Г, < 0,2 - 10 "зви т, > 0,3 - 107? сек. Такие оценки 
согласуются C экспериментальными данными по ридберговским атомам (CM., на- 
пример, [4]. 

Таким образом, если начало отсчета времени (t= 0) выбрать с момента 
«приготовления» нашего начального состояния, то при движении электрона к 
протону на пути в | см с положительной энергией (т. e. в непрерывном спектре) 
в интервале времени от 0 до ~ 0,5 - 1077 сек будет происходить вследствие после- 
довательных радиационных захватов электрона в связанные состояния формиро- 
вание возбужденного атома водорода с вкладом всех ридберговских состояний 
вплоть до п = 107? и, естественно, всех нижележащих возбужденных и основного 
связанного состояния. С течением времени все возбужденные состояния, которые 
в реальности являются метастабильными, будут постепенно распадаться. 

Теперь обратим внимание на возможность принципиально новых, ранее не 
обсуждавшихся процессов. При вылете фотона, излучаемого при радиационном 
захвате электрона на какой-либо из уровней атома (скажем, п-й), имеется нену- 
левая вероятность фотовыбивания электрона (внутренней ионизации атома) из 
вышележащих уровней (n + v, v = 1, 2,...), входящих в общее соотношение (6), c 
переходом его (электрона) в непрерывный спектр энергий, причем с большим зна- 
чением средней кинетической энергии, чем вначале (Et) = E, + (Е +» — En), 
v= 1,2, ...). Вьбитьй электрон будет вначале удаляться от протона до бесконеч- 
ности, постепенно замедляясь (строго говоря, следует говорить о движении волно- 
вого пакета, описывающего вероятность движения электрона). При этом по пути 
движения имеется ненулевая вероятность радиационного захвата в возбужденные 
и далее ридберговские состояния увеличивающегося радиуса. Вылетающие при 
таких захватах фотоны будут в свою очередь выбивать электрон из еще более 
вышележащих уровней. И этот процесс эволюции вероятности образования воз- 
бужденного атома водорода будет продолжаться до бесконечности с возрастанием 
числа формируемых связанных состояний за счет радиационного захвата в про- 
цессе движения электрона, последовательно выбиваемого из все более высоких 
возбужденных состояний последовательно вылетающими фотонами. Следователь- 
но, и размеры возбужденного атома водорода будут расти до бесконечно больших 
значений, однако, с сильным затуханием вероятности формирования такого ато- 
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ма (в силу сравнительно небольшой вероятности радиационного захвата и фото- 
выбивания). 

Разумеется, наряду с внутренней ионизацией должны иметь место и про- 
цессы внутреннего возбуждения атома (тоже с увеличением его размеров), когда 
фотон, излучаемый при радиационном захвате электрона на тот или иной уро- 
вень атома, возбуждает вышележащие уровни в еще более высоколежащие и 
далее вылетает из атома с потерей энергии, испытывая таким образом неупругое 
рассеяние. 

Во всех этих процессах, как в принципе нетрудно показать с помощью эле- 
ментарных подсчетов, закон сохранения энергии выполняется, а рост энергии вы- 
биваемых электронов и излучаемых при их последующем захвате фотонов ком- 
пенсируется увеличением относительного вклада фотонов с меньшей энергией. 

При t > © вероятность формирования ридберговского атома практически бес- 
конечных размеров затухает до нуля, а вероятность найти атом в основном состо- 
янии (т. е. минимальных размеров) достигнет некоторого конечного максималь- 
ного значения, определяемого амплитудами вероятности первого процесса захва- 
та, а также радиационными ширинами всех уровней, возбуждаемых при первом 
захвате. 

В принципе наряду с описанными процессами имеется ненулевая вероятность 
еще одной цепочки многократно повторяющихся процессов внутренней иониза- 
ции (или внутреннего возбуждения) — фотовыбивания (или фотовозбуждения) 
электрона с промежуточных уровней излучением, испускаемым при девозбужде- 
нии (радиационном распаде) более высоких уровней на самые низкие уровни, что 
тоже дает вклад в вероятность формирования ридберговских атомов увеличиваю- 
щихся размеров. 

В случае если начальная скорость электрона направлена под ненулевым уг- 
лом к начальному расстоянию между электроном и протоном, движение, которое 
в классическом случае будет проходить в течение гораздо большего времени (в 
~ 4 · 10? раз дольше при угле 90°), и в квантово-механическом случае будет идти 
также гораздо медленнее в соответствии с теоремой Эренфеста. К тому же и 
пересечения состояний y,(r) будут происходить не при движении по радиусу к 
центру, а при движении по спирали, что дополнительно увеличит время движе- 
ния между областями локализации соседних состояний. Поэтому в течение на- 
чального этапа движения (по спирали) в первично образующемся атоме будет со- 
существовать до радиационного распада меньшее число ближайших друг к другу 
уровней атома, чем в случае движения электрона и протона по прямой навстречу 
друг другу. В остальном же принцийиальных отличий эволюционной ка ртины не 
будет. 

Полноты ради отметим, что даже после первой внутренней ионизации выби- 
тый электрон может удалиться на очень большое расстояние и при возвращении 
после остановки он может захватиться в состояния ридберговского атома очень 
больших размеров. Соответственно и в результате первого внутреннего возбужде- 
ния может образоваться ридберговский атом очень больших размеров. Более того, 
энергия выбитого электрона (при первой или последующей внутренней иониза- 
ции) может оказаться такой большой, что электрон улетит с положительной ки- 
нетической энергией на бесконечность и не возвратится. 

Учет тормозного излучения, в силу его малой вероятности даже по сравне- 
нию с радиационным захватом, фотовозбуждением и фотовыбиванием, качест- 
венно не меняет общей эволюции столкновения и происходящих при этом про- 
цессов, количественно слегка ослабляя рост размеров последовательно формируе- 
мого атома и замедляя эволюцию системы, 

4. Заключение. Проведенное в настоящей работе концептуальное и качест- 
венное рассмотрение процесса столкновения электрона и протона позволило (на 
фоне подтверждения давно известных фактов невозможности формирования ус- 
тойчивого и даже метастабильного атома водорода в классическом подходе и не- 
нулевой вероятности формирования суперпозиции устойчивого основного и отно- 
сительно метастабильных возбужденных состояний атома водорода) впервые тео- 
ретически установить (предсказать) многократные процессы формирования с 
ненулевой (хотя малой и уменьшающейся с течением времени) вероятностью воз- 
бужденных и высоковозбужденных (ридберговских) состояний атома водорода с 
постепенно увеличивающимися энергиями возбуждения и размерами. Эти про- 
цессы, в свою очередь, теоретически устанавливают (предсказывают) эффект за- 
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медления и даже уменьшения вероятности окончательного формирования основ- 
ного состояния атома водорода за счет внутреннего возбуждения и внутренней 
ионизации, конкурирующих с процессами прямого радиационного распада уров- 
ней. Подчеркнем, что если отбросить влияние злектромагнитного вакуума (или, 
другими словами, без внешнего воздействия, обусловливающего радиационные 
захваты), начальные электрон и протон не смогут образовать ни устойчивого, ни 
метастабильного атома водорода даже в квантовой теории. 

В заключение считаем своим приятным долгом выразить благодарность 
Е. А. Андрееву за участие в постановке задачи, а также Л. С. Брижик и 
А. К. Зайченко за помощь и участие в работе. 


ПОРІВНЯЛЬНИЙ ЯКІСНИЙ АНАЛІЗ ЕВОЛЮЦІЇ СИСТЕМИ ЕЛЕКТРОНУ ТА ПРОТОНУ, ЩО 
СТИКАЮТЬСЯ, ЗА УЧАСТЮ ЕЛЕКТРОМАГНІТНИХ ВИПРОМІНЮВАНЬ ЗА КЛАСИЧНИМ 
ТА КВАНТОВИМ ПІДХОДАМИ 


В. С. ОЛЬХОВСЬКИЙ, С. П. СІТЬКО 


На якісному та концептуальному рівні проаналізовано еволюцію системи електрону та протону, що 
стикаються, за участю гальмуючого електромагнітного випромінювання та радіаційних (фотонних) 
переходів у класичному та квантовому підходах. В межах квантового підходу теоретично встановлено 
наявність процесів багатократного формування збуджених та вельмизбуджених (рідбергеровських) 
станів атому водню з поступово зростаючими з плином часу енергіями збудження та розмірами. Про- 
цеси подібного гатунку 3 прямими наслідками процесів внутрішнього збудження та внутрішньої 
іонізації утвореного атому фотоном, що випромінюється внаслідок попередніх актів захвату електрону 
протоном на рівнях атому, що лежать нижче. Ланцюги всіх подібних процесів, що конкурують з про- 
цесами прямого радіаційного розпаду усіх первісно збуджених рівнів, сприяють процесу гальмування 
остаточного формування основного стану атому водню. 


СРАВНИТЕЛЬНЫЙ КАЧЕСТВЕННЫЙ АНАЛИЗ ЭВОЛЮЦИИ СИСТЕМЫ СТАЛКИВАЮЩИХСЯ 
ЭЛЕКТРОНА И ПРОТОНА С УЧАСТИЕМ ЭЛЕКТРОМАГНИТНОГО ИЗЛУЧЕНИЯ 
В КЛАССИЧЕСКОМ И КВАНТОВОМ ПОДХОДАХ 


В. С. ОЛЬХОВСКИЙ, С. П. СИТЬКО 


Качественно и концептуально проанализирована эволюция системы сталкивающихся электрона и про- 
тона с участием тормозного электромагнитного излучения и радиационных (фотонных) переходов в 
классическом и квантовом подходах. В квантовом подходе теоретически установлено наличие процес- 
сов многократного формирования возбужденных и высоковозбужденных (ридберговских) состояний 
атома водорода с постепенно увеличивающимися энергиями возбуждения и размерами с течением вре- 
мени. Такие процессы являются прямым следствием процессов внутреннего возбуждения и внутренней 
ионизации образовавшегося атома фотоном, испускаемым в результате предшествующих актов радиа- 
ционного захвата электрона протоном на нижележащие уровни атома. Цепочки всех таких процессов, 
конкурирующие с процессами прямого радиационного распада всех первоначально возбужденных 
уровней, приводят к замедлению окончательного формирования основного состояния атома водорода. 
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Based on experimental results, the mathematical model of pacemaker bursting activity has been 
formulated. The model includes slow-wave generating mechanism, spike-generating mechanism, 
inward calcium current, accumulation Ca-ions in cytoplasma and their absorption by specific 
depots. Stationary [Ca? * Ji,- inhibited calcium current in this model allows us to obtain a more 
complete model of bursting activity. The model demonstrates the majority of known phenomena 
and allows us to investigate the role of separate components in generation of bursting activity in 
nerve cells. Routes to chaos and mode transitions depending оп. control parameters and stimuli 
have been observed. 


1. INTRODUCTION 


A nerve cell is the structural and functional unit of the nervous system, and its 
characteristics determine to a large extent the properties of neuronal assemblies which 
are responsible for perception, information processing and storage. Model neuron in 
neuronal networks is often described as a threshold element with several inputs, and 
its responses are represented as binary, linear, exponential or sigmoidal functions. 
Indeed, real neuronal responses to input signals may be rather complicated and non- 
predictable due to non-linear characteristics of a nerve cell. Single nerve cells in CNS 
of different animals can spontaneously generate slow membrane potential oscillations 
which are the cause of spikes bursts generation [1—6 |. This pattern is called bursting 
activity. At present, physiological role and membrane mechanisms of bursting activity 
of nerve cells are not quite clear. 


2. MECHANISMS OF PACEMAKER ACTIVITY GENERATION 


Rhythmic processes in a living organism has always attracted attention of re- 
searchers. The goal of our work is to study membrane mechanisms of generation of 
oscillatory bursting activity in CNS nerve cells of snail, Helix pomatia with sub- 
sequent mathematical modeling these processes and computer simulation. 

Fig. 1 demonstrates examples of bursting activity registration in nerve cells RPal 
(right parietal ganglion, cell 1 [6р in two different preparations. It is seen that 
membrane potential has wavy changes and on reaching the threshold value for action 
potential (AP) generation a cell generates burst of spikes (Fig. 1A). 

Amplitude of membrane potential (MP) waves in RPal neuron in different 
preparations can vary from 1 to 30 mV, durations of both AP burst and interburst 
intervals range from | to 200 s. Although parameters of MP waves remain constant in 
a majority of preparations, nevertheless one may sometimes observe and modulation 
of bursting activity (Fig. 1B). Analysis of physical models of MP waves generation are 
given in [5] and their mathematical versions in [7]. 


50mV | 
20s 


Fig. 1. Bursting activity of the identified bursting pacemaker neuron RPal in different preparations of snail 
Helix pomatia. 


55 


А. О. KOMENDANTOV, М. I. KONONENKO 


It is still Jess known about electrical bursting activity in integrative information 
processing in brain. This type of activity is believed to provide a more effective release 
of neurotransmitter by terminal of bursting neuron [8, 9]. This neuron is suggested 
to play the role of a very sensitive device [10], and very small influences on cell 
membrane as both depolarization and hyperpolarization lead to a considerable change 
of the activity pattern. ; 

Investigations of bursting activity in neurons of different vertebrates and 
especially invertebrates show that depolarizing phase of MP wave is due to specific 
stationary potential — dependent sodium conductance gu, (V). Therefore, the mem- 
brane of a bursting neuron has a constant tendency for depolarization and its current- 
voltage relationship (CVR) is always below zero. CVR has a so called negative 
resistance region (NRR) [11 ]. 

As for the mechanism of generation of MP wave hyperpolarizied phase, there is 
no general point of view in this case. In our studies it has been established that 
potential- and time-dependent outward (Ej, = —58) current, Ip, activated by hyper- 
polarization exists in the bursting snail neuron [12]. Having achieved the maximal 
value, this current decayed with a time constant of about 1 sec. At present, the ionic 
nature of this current is not understood. Thus, there are reasons to believe that 
mechanism of slow-wave generation is based on the presence of specific stationary 
potential-dependent Ма’ conductance, gy,(V), taken together with stationary 
potential-independent Ма" and K* conductances (gna + gx), and hyperpolarization- 
activated conductance, gg. These conductances in intact bursting neuron are regulated 
by neuropeptide concentration near the cell receptors. 


3. MATHEMATICAL MODEL OF PACEMAKER ACTIVITY IN A BURSTING 
NEURON 


Basing on the experimental results and hypothesis about slow-wave generation 
mechanisms in RPal Helix neuron mathematical model of pacemaker activity was 
formulated [13]. The model includes: a) slow-wave generation mechanism (к, Ina, 
Ina(V), Ip); spike-generating mechanism (the Hodgkin-Huxley sodium TT X-sensitive 
(INacrrx)) and potassium TEA-sensitive (Ixcrea)) currents); intracellular Ca ions, 
[Ca?* ],,, and their buffering; calcium currents (transient potential-dependent (Ica) 
and stationary [Ca?* J,-inhibited (Ic, c.) ]. Equations and parameters are presented 
in Appendix. The model has been studied by computer simulation. 

The model successfully reproduces different patterns of electrical activity of 
bursting neurons under physiological conditions — monomodal beating activity, 
bursting activity with different parameters of oscillations and silent state without 
activity. It also describes such known experimental phenomena as adaptation of 
electrical activity to polarizing current [14], stationary ionic currents which appears 
under voltage clamp conditions [11], hysteresis of current-voltage relations of mem- 
brane [12]. 


ое Fig. 2. Bursting pattern оѓ 
electrical activity in the model 
(Сази) [13]. Тіте course of membrane 
potential (top graph) and 
[Ca?*]i, (bottom graph) 
0.00 during the first 50 s. 
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Fig. 2 shows. the time courses of MP апа [Ca?* |, in the model. In presented 
model, inward Ca ions under AP generation play a key role for adaptation of a cell to 
polarizing current and stabilization of electrical activity under the changes of some 
system parameters. These processes are mediated by influence of Ca ions on 
properties of cellular membrane. 


4. CHAOS IN BURSTING NEURON 


The electrical activity of RPal cell is determined by a set of currents running 
through a number of different types of ion channels across the membrane. Presence of 
the processes with different time scales is a possible cause of origin of the irregular 
regimes. Indeed these regimes arise under certain conditions and parameters of the 
model [15]. Parameters for a chaotic mode are presented in Appendix. The ones of 
their number are chemosensitive conductances (ема, Єв and Бу (У)) of transmembrane 
ionic currents forming slow-waves. These conductances are chosen as bifurcations 
parameters. The model demonstrates regular beating, bursting and chaotic activity 
depending on values of gna, gu, (V) and ев. Interspike intervals (ISI) which аге a basis 
for information processing in neural nets were used for analysis of different types of 
neuronal activity. 
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Fig. 3. Bifurcation diagram for the model neuron. 
Bifurcation parameter (abscissa) is chemosensitive 0 


conductance, gNa( V). The ordinate is ISI (s). 


Fig. 3 demonstrates the bifurcation diagram obtained for вка (У) as a bifurcation 
parameter. One can see on the diagram the transformation of beating activity into the 
bursting one taking place through a cascade of period-doubling bifurcation and region 
with chaos. Time series analysis as well as quantitative and qualitative diagnoses of 
this pattern show that routes to chaos were generated. It should be noted that the 
bifurcations demonstrate some sudden and intermittent transitions from periodic to 
aperiodic behavior. These patterns of pacemaker activity in the model are very 
sensitive to initial conditions and that is a feature of deterministic chaos. Therefore, 
parameter-independent mode transitions from periodicity to chaos (or quaziperiodi- 
city) and vice versa evoked by short-lasting polarizing current or shift of con- 
centration of intracellular Ca ions have been observed (Fig. 4). These influences can 
be interpreted as activation of synaptic inputs or as functioning of intracellular Ca 
depots. 


Fig. 4. The mode transition under short depo- 
larization [15]. Time course of membrane potential. 
Abscissa: time (s); ordinate: V (mV). Horizontal thick 
bar: current pulse (-0.5 nA, 2s) is applied. 
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5. CONCLUSION 


The experiments conducted evidence that oscillations of membrane potential are 
connected with existence in membrane specific channels for inward sodium currents 
and potential-dependent channels for outward current. The results obtained let us un- 
derstand clearer the mechanisms of oscillatory electrical activity in nerve cells and the 
role of separate membrane and intracellular processes in these events. Our simulation 
studies show that deterministic model of pacemaker activity in bursting neuron under 
certain values of parameters exhibit periodic and chaotic behavior. A good agreement 
with experimental results has been observed. Some intact neurons in CNS exhibit 
bursting, beating or silent activity depending on their state, stimulation of peptidergic 
interneurons. It corresponds to changing of gna, gw,(V) and gg in physiologically 
significant regions in the model. Our simulation studies allow us to hypothesize that 
chaotic patterns of electrical activity in intact neurons could exist due to the condition 
of a cell membrane as well as properties of intracellular calcium depots. 

Chaos in a nerve cell could not be the consequence of convergention of excitatory 
and inhibitory postsynaptic potentials of this cell working in neuronal assemblies. A 
single neuron with deterministic characteristics can function as a chaotic generator 
and be a source of «random processes». It demonstrates that a neuron functions in the 
unity of deterministic and probabilistic mechanisms. Ability for chaotic behavior is 
one of the fundamental properties of excitatory systems with complex type of 
electrical activity which have slow and fast components. That corresponds to variables 
for membrane conductances on two time scale in models [16—20]. In our model, 
under certain values of their parameters, the small variations in initial conditions or 
model parameters can produce essential changes in behavior. Similar results are 
obtained in other models of electrical activity of excitable cells [21, 22]. Our studies 
demonstrate the essential role of intracellular calcium ions in mode transitions. 


APPENDIX 
Mathematical model of pacemaker activity in bursting neurons 
Basic equation governing membrane behaviour is 
— Cyd V/dt = Іма(ттх) + Їк(твА) + Ik + Ina + Їма(му + Ig + Ica + Lc; са 


Here, C,, is membrane capacity, V denotes membrane potential. 
a) spike-generating mechanism (the Hodgkin-Huxley sodium TTX-sensitive and 
potassium TEA-sensitive currents): 


Іма(ттх) = BNa(TTX) m?h (У — Vna); [kerea) = EK(TEA) n*(V - Ук); 
dm/dt = (1/(1 + exp(— 0.4(V + 31))) — m)/0.0005; 

dh/dt = (1/(1 + exp(0.25(V + 45))) — h)/0.01; 

dn/dt = (1/(1 + exp(— 0.18(V + 25))) — n)/0.015; 

b) slow-wave generation mechanism: 

In (V) = gRA(V)(1/(1 + exp(— 0.2(V + 45))))(У — Ука); 

Ik = g&(V — Ук); ha = gw(V — Ум); Ip = gomghs(V — Ув); 
dmp/dt = (1/(1 + exp(0.4(V + 34))) — mg)/0.05; 

dhy/dt = (1/(1 + exp(— 0.55(V + 43))) — hy/1.5; 


с) calcium currents (transient potential-dependent (Ic, and stationary 
[Ca?* lin — inhibited (Ica-ca): 


Іса = gc4mc,2(V — Уса); dmc,/dt = (1/(1 + exp(— 0.2V)) — me,)/0.01; 


"T ЗЛ З MESI S 1 
[се-се = 8-611 E exp(— 0.06(V + 45)) 1 + скр (Са -Ay V T Veo) 


d) intracellular Ca ions and their buffering: 
d[Ca ]/dt = — Ic,/2Fv - ks[Ca]; у = 4782/3; 
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Here, [Ca] is [Ca?* 1, Е is Faraday number, v is the volume of the cell, ks is 
the rate constant of intracellular Ca-uptake by a buffering systems. 


Parameters for Eqs. 
Vua = 40 mV; Ук = - 70 mV; Ve = — 58 mV; Vea = 150 mV; 
Си = 0.02 «Е; В 21004; К; = 0.05 !/s; kg = 301/mM; В = 0.02 mM. 


Chaotic mode: 

gk = 0.25 1/MQ; рх, = 0.0231 1/MQ; gn.(V) = 0.11 1/MQ; 
єв = 0.1372 1/MQ; ZNa(rrx) = 400 1/MQ; gk(rea) = 10 1/MQ; 
56. = 1.51/MQ; р, c, = 0.02 1/MQ. 


Initial conditions: V = —55 mV; [Ca] = 0 mM. 
НЕЛІНІЙНІСТЬ, ОСЦИЛЛЯЦІЇ TA ХАОС У ПЕЙСМЕКЕРНИХ НЕЙРОНАХ 
А. O. КОМЕНДАНТОВ, H. І. КОНОНЕНКО 


На базі експериментальних результатів було сформульовано математичну модель вибухової активності 
пейсмекеру. Ця модель об'єднує механізм генерації повільних хвиль, механізм генерації спайків, 
внутрішній потік кальцію, акумуляцію іонів Са в цитоплазмі та їхню адсорбцію у спеціальних депо. 
Стаціонарний кальцісвий потік, що Нибований са? hin, в даній моделі дозволяє одержати більш 
цільну модель пачечної активності. Ця модель демонструє більшість з відомих феноменів і дозволяє 
дослідити роль окремих компонентів у генеруванні пачечної активності у нервових клітинах. Було 
досліджено спенарії переходу до хаосу та модових переходів у залежності від значень керуючих пара- 
метрів і різних стимулів. 


НЕЛИНЕЙНОСТЬ, ОСЦИЛЛЯЦИИ И ХАОС В ПЕЙСМЕЙКЕРНЫХ НЕЙРОНАХ 
А. О. КОМЕНДАНТОВ, H. И. КОНОНЕНКО 


На основе экспериментальных результатов была сформулирована модель пачечной активности пейсме- 
кера. Данная модель включает механизм генерации медленных волн, механизм генерации спайков, 
входящий кальциевый TOK, аккумуляцию ионов Са в цитоплазме и их адсорбцию в особых депо. Ста- 
ционарный кальциевый TOK, ингибированиый [Са] в данной модели позволяст нам получить более 
цельную модель взрывной активности. Данная модель демонстрирует большинство из известных фено- 
менов и позволяет исследовать роль отдельных компонентов в генерировании пачечной активности в 


нервных клетках. Были исследованы сценарии переходов к хаосу и модовых переходов в зависимости 
от значений управляющих параметров и различных стимулов. 
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PROPAGATION AND DISTRIBUTION OF NEURAL SIGNALS: 
A MODELING STUDY OF AXONAL TRANSPORT 


FRANK RATTAY 


Technical University Veinna, A-1040 Vienna, Austria 


ABSTRACT. The functional properties of a neuron can be simulated by a system of electrical 
circuits. This technique is of use for natural situations as well as for electrical stimulation. In 
particular, a combined model for the mammalian neuron is proposed, which takes care about 
the different structures and properties of the neural membrane of the axon. The highly energy 
consuming nonmyelinated terminal region can be stimulated by the output of a myelinated axon 
if the last internodes have shorter distances. 


INTRODUCTION 


In contrast to other cells neurons are remarkable because of their shapes showing 
a very complicated fine structure and because of the electrical properties of the 
membrane which separates the ions of the intracellular and the extracellular fluids. 
Neural information is carried by the propagation of membrane voltage changes along 
the neuron. The electrical properties of the neural membrane are responsible for a 
graded signal transfer within short distances and the "ай or nothing’ law for 
information propagation by neural impulses which reaches distant locations. 

The main elements of a neuron are the soma (cell body), the dendrites and the 
axon (Fig. 1). Dendrites and soma are usually covered with synapses from other 
neurons (input zone), whereas the task of the nerve fiber (axon) is the transport and 
the distribution of neural signals. The mammalian nerve fibers with diameters from 
1—22 um are usually myelinated, i. e. only within the nodes of Ranvier’, which have 
a length of about 1 ит, the membrane is active. The internode is covered with many 
layers of membranes of Schwann cells (in the periphery) or glia cells (in the brain). 
Thus the inner diameter d of an axon represents about 70 % of the outer diameter D. 
The internode has a length of about 100 times D. The insulating properties of the 
internode allow a 20 ит mammalian fibre to send an action potential (AP) with 120 
m/s, which is essentially faster compared even to the unmyelinated giant axons of 
squid which reach diameters up 1o 1 mm. 

The individual behavior of a neuron at different sections is a result of the 
electrical properties of its membrane which depend essentially on the number and 
types of embedded ion channels and the degree of insulation by sheets of myelin. 

The functional elements of a neuron have been subjects of detailed investigations 
and there exist a plenty of modeling studies with analytical and numerical results. 
Our knowledge about the influence of synaptic activity, the passive properties of the 
dendritic tree and the famous modeling work of Hodgkin and Huxley inspired 
scientists from many fields [1—5]. New disciplines like neural network techniques in 
computer science have been generated and several hypotheses on brain functions use 
some of the empirical facts and speculate on their interactions. Most of the brain 
models use the Hodgkin — Huxley membrane dynamics (or even simpler ones) in 
order io explain e. g. regular or chaotic system behavior in spiking patterns and data 
processing. 


/ 222 -4 Fig. І. Functional scheme of a neuron. Excitation of the 
< -neuron is a consequence of the stimulation by activated 


es 


a 5 WODE OF RANVIER `\\ Synapses in the input region at the dendritic tree and the 
soma. When excitation exceeds the threshold value, a single 

/ д nerve impulse or even a train of such impulses is generated. 
— | A neural signal can propagate with high velocity along the 


myelinated axon. The nonmyelinated branching part reaches 

the target region, where the neural impulses are multiplied 
CONDUCTILE острот and they distribute the information via synapses to other 
REGION REGION REGION Cells. Besides the main axon sometimes there exist collateral 
branches. 


60 


Propagation and Distribution of Neural Signals: A Modeling Study of Axonal Transport 


In contrast to such models which loose: many of the relevant biological features 
this paper presents a general method to simulate neural reactions by using electric 
networks. An illustrative example is the situation in the terminal section. of 
mammalian axon. In the terminal section. the axon changes. from the myelinated 
structure to a branching unmyelinated fiber. We will show. that the simple 
combination of existing “minimodels’ for the myelinated and the unmyelinated axon 
fails to describe signal transport in this region. However, reducing the internodal 
length in the last sections allows the unmyelinated fiber to overcome this problem. 


MEMBRANE MODELS 


The most important element in rapid biological signal processing is the neural 
membrane. A patch of the membrane of a neuron can be described as an electrical 
circuit consisting of a voltage source, which accounts for the resting voltage, a 
capacitor (1—3 нЕ/ ст?) and a resistor. Those parts of the membrane which contain 
ion channels can be active in order to amplify the input signal. The conductance of 
the active membrane changes within a large range depending on the open/closed 
status of the ion channels. 

Four membrane models are available to simulate the mammalian axon [5—11 |. 
In this article we will concentrate оп the two of them which are similar in their 
equations: the Hodgkin-Huxley model (HH) and the CRRSS model (Table 1). Note 
that the role of potassium current in the repolarization process is of high importance 
in the HH model, whereas it is neglected im the CRRSS model [12, 13]. 


| HODGKIN-HUXLEY MODEL 


У = [-gnamPA(¥ — Ума) — grn (V — Ук) | SWEENEY et al (1987) transformed the origi- 

-gi(V — VE) +] (HH-1) паї data of Снід, RITCHIE, Rocert & STAGG 

т = [- (ањ + вт) m+ даті: k (HH-2) || (1979) from experimental temperature Т = 14°C 

n= [- (а, +8.) поволі k (HH-3) || to T = 37°С. The following equations are called 
h=[-(on+ Br) h+ an] k HH-4) | after the investigators the CRRSS model. 


(HH-4) 
|with the coefficient. k for temperature Т (in ^C). V = | -дкат' КУ - Ума) — У — и) ХЕ с 
k = 36.17-59.63 (HH-5) CRRSS-1 
ЗА m= (оњ + Es): m + Om CRRSS-2 


h = -(оһ + №) -В+аь (CRRSS-3) 


oe 
CRRSS MODEL 


and 
with 
ат = (2.5 - 0.1V)/(ezp(2.5 — 0.1V) - 1) 
Вл. 1- exp(-V/18) 
an = (1 — О EV)/(18.ezp(1 - 0.1V) - 1) ‚ Om = (97 + 0.363V)/(erp((31 — V)/5.3) + 1) 
An, = 0.125 -exp(—V/80) Ви = ат [lexpe(V — 23.8)/4.17)) 
an = 0.07 -exp(—V/20) en = Baf(ezp((V — 5.5)/5)) 
Bi, = M (ezp(3 — 0.1V) + 1) i Ba = 15.6/(1 + с2р((24 - V)/10)) 
Vrest = —70[mV] Мея = —80[mV] 
Ума = 115[mV], Vk = —12[mV], Vr = 10.6[mV] || Ума = 115{mV], V; = 0.01[mV] 
gna = 120607 сто ?), ок = 36[607 стт 2] | ома = 1445[kQ7 em], gp = 128[60- erm?) 
gn, = ЗО стт 2), с = Ци сті) с = 2.5[nF/cm?] 
m(0) = 0.05, n(0) = 0.32, A(0) = 0.6 m(0) = 0.003, À(0) = 0.75 


Table 1. Comparison of equations for membrane dynamics of unmyclinated and myelinated axons. 
V denotes the reduced membrane voltage in mV (У = Vinembrane — Vrest); Ума, Ук, VL voltages 
across the membrane, caused by different (sodium, potassium, unspecific) ionic concentrations in- 
side and ouiside the aron; gNo, ge, 9L Maximum conductance of sodium, potassium and leakage 
per cm? of membrane; т, п, h probabilities for ionic membrane gating processes; с, В opening and 
closing rates for ion channels; c the capacity of membrane per ern? ; m(0), n(0), A(0) values at the 
resting condtion where V(U) —0; i, stimulus current. 

Noic the high conductances and the missed potassium currentis in the myelinated mammalian node 
as described by the CRRSS model. 
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Original measurements on squid axons have been done at a low temperature 
(6,3 °C) to keep the axon working for several hours. In order to use the HH model for 
higher temperatures it is necessary to find a proper value for the coefficient k, 
according to the faster reactions of the gating processes in the ion channels [14]. Qio 
(which is the acceleration factor of a process at an increase of 10 °C) can be assumed 
to be 3 for all the gating variables only for lower temperatures. With k = 12 the HH 
model can describe the membrane of the mammalian nonmyelinated axon [9]. 


LUMPED CIRCUITS FOR SIMULATING A FIRING AXON 


In order to simulate a myelinated nerve fiber and its branches (conductile and 
output region in Fig. 1) these structures are divided into cylindrical segments. 
According to Fig. 2, every segment is represented by an electric circuit consisting of 
membrane conductor and capacitance, a voltage source taking into account the effect 
of different ionic concentrations of the inner- and extracellular fluids and the 
axoplasm conductances to the neighbored segments. The membrane conductance 
essentially depends on the gating mechanisms of the ion channels. 

For simplicity we assume the internode to be a perfect insulator and the 
extracellular potential to be at constant level У,, = 0 for every segment. The reduced 
voltage V, across the membrane of the n-th segment becomes 


V, = Vin > Vrest- (1) 


UNMYELINATED 
FIBER 


MYELINATED 
FIBER 


OUTSIDE 


MEMBRANE 


INSIDE 


BRANCHING 
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4 
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1 
, 
1 
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Fig. 2. Axons without and with insulating sheets of myelin are modeled by the same electrical network. Both 
types of nerve fibers are segmented into cylinders of length Ax and diameter d. For myelinated fibers the 
segment length Ax is the internodal distance but the active membrane exists only at the shaded area in the 
node of Ranvier (nodal length L is about 1 шт) and there the ion currents can cross the membrane. The 
membrane of every segment is simulated by an electric circuit consisting of capacity Cm, voltage source, and 
nonlinear resistor Gm Ve, and Vin are the external and the internal potentials at the n-th segment, Ga denotes 
the conductance of axoplasm between two segments. The begining of a branching part of a nerve fiber is 
shown in the electric network by broken lines. 


62 


Propagation and Distribution of Neural Signals: A Modeling Study of Axonal Transport 


The currents of the n-th segment are caused by membrane capacitance (C,), 
membrane conductor (G,) and the axonal conductors (G,,—1, Gan) to the neighbored 
segments: 
dV, 
dt 


C, > + С, з (У, = Vi-1) л Сап 1 Ж Vin Vint) а Gan T (...) = 0 (2) 
where С, · У, is the ion current through the n-th segment and for (...) the current to 
an eventually existing branching fiber has to be substituted. 

Because the equations for the ion currents in Table | are formulated as current 
densities, it is convenient to transform (2) by setting 


ла? 
4р Ах 


G, = апа С = слаг 

where р; is the specific resistance of axoplasm, which is about 0.1 КО. cm, с is the 
capacity per ст? of free membrane, and L is the nodal gap width. We can treat the 
case of the myelinated and the unmyelinated fiber in parallel, if we set L = Ax for 
unmyelinated fibers. If d and Ax are constant along the axon, (2) together with units 
used in membrane models becomes 


Ev налета 


ms 
p ст | 
4p; | КОст | - Ax [em] ст] 


- (V,—1 = ЗУ, + Va+1) [nV] + (3) 
І, ША] 


ad lem Lien | / ст ЕЙ ст?) 


where /,, is the stimulating current, which is often applied in modeling studies in the 
first segment only in order to generate a propagating action potential there. 

By analysis of eqn. 3, it can be seen, that from assuming 1. the internode to be 
an ideal insulator 2. constant L and 3. internode length proportional to d it follows 
that propagation velocity of myelinated fibers becomes proportional to d, whereas it is 
proportional to Vd for unmyelinated fibers. Furthermore, it is seen that the shape 
and duration of an action potential is independent of fiber diameter for both the 
myelinated and the unmyelinated fiber. For details see Ref. 9, p. 110—115. 

The nodes of Ranvier have a length L of about 1 ит, and they are very short 
compared with the internode, which obeys the empirical rule Ax = 100 D giving a 
length up to 2 mm. As a consequence, propagation velocity is rapid and energy 
consumption low compared with the situation of the non-myelinated fiber for fibers 
with d > 0.25 um. However, for the theoretical value d= 0.25 um nonmyelinated 
axons become faster than myelinated ones because of their nonlinear diameter 
velocity relation [9]. This may be the reason why the thinny fibers in the terminal 
branching regions are unmyelinated (comp. Fig. 1). 

Quick and coworkers have reported that the last few internodes before reaching 
the unmyelinated terminal are shorter [15]. By computer simulation it is 
demonstrated that this shortening of the internode is necessary to allow an action 
potential to propagate into the unmyelinated part of an axon. Fig. 3 and Fig. 4 show 
a computed action potential which is generated in the beginning of a 12 node 
myelinated axon. This signal can pass the border into the unmyelinated part because 
the last internodes are gradually shortened until half of the standard length is 
reached. The lowest part of Fig. 4 shows the geometric situation. Simulation was done 
with the models described in Table 1, however the unmyelinated axon was computed 
with k=12 and doubled channel density resulting in. gy, = 240 КО" ст" 2, 


вк = 72 КО" ст"? and g = 0.6 KO! cm. 
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Although the ion current dinsities in the firing membrane at the nodes of Ranvier 
are about ten times that of the unmyelinated axon (comp. the lower curves in Fig. 3 
and Ref. 11), the unmyelinated membrane needs essentially higher membrane 
currents per cm (Fig. 4). The reason is that the first 7 internodes in this model have 
a length of Ах = 857 ит and for the last internode Ах = 428 ит, but with 
L= 1.5 ит it turns out that only 0.175 to 0.35 У, of the myelinated axon have an 
active membrane. 


DISCUSSION 


The technique of equivalent electric circuits proposed in this paper is a general 
tool which is of use to simulate arbitrary elements of a neuron. Hot spots, i. e. active 
membranes in the dendrites can be simulated as well as the membrane at the axon 
and at the soma, if the density and types of ion channels are known. Furthermore, 
many results have been obtained by taking the extracellular voltage У,, in Fig. 2 
according to an applied electric field in order to simulate artificially evoked neural 
reactions [9, 16]. 

Combining a simple model for the myelinated fiber with the double-channel- 
density Hodgkin-Huxley membrane demonstrates that shortening the internodes 
allows an action potential to travel into final regions of an axon. Additional 
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investigations are necessary to show the influence of myelin [17], which is no a 
perfect insulator and the temporal reactions of the combined fiber to firing sequencies, 
because current leakage through myelin as well as differences in refractory behavior 
can lead to loss of action potentials in the terminal section. 
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Fig. 4. Voltage and current distribution across the axon membrane as a function of fiber's length coordinate. 
The same situation as in Fig. 3, i. e. a mixed type of 1 cm fiber with gradually shortening of the last 4 
internodes. At the bottom of the figure the position of the 12 nodes is marked by circles and the filled part 
shows the unmyelinated axon. The curves show snapshots in intervals of 0.1 ms of membrane voltage and total 
membrane current per cm of axon (lion + icapaciry in pA/ cm). Note the low propagation velocity and the high 
energy consuming in the unmyelinated part of the axon. 


ПОШИРЕННЯ TA РОЗПОДІЛ НЕЙРОННИХ СИГНАЛІВ: МОДЕЛЮВАННЯ AKCOHHOTO 
ТРАНСПОРТУ 


ФРЕНК РЕТТЕЙ 


Функціональні властивості нейрону можуть модулюватися за допомогою системи електричних ланцю- 
гів. Ця техніка придатна як для природних процесів, так і процесів електричного стимулювання. 
Зокрема, запропоновано об'єднану модель нейрону ссавця, яка бере до уваги різні структури та вла- 
стивості нейронної мембрани аксону. Немієлізована термінальна ділянка, яка потребує більшої енергії 
може бути стимульована сигналом мієлізованого аксону при умові, що останні міжвузлові відстані 
менші попередніх. 
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РАСПРОСТРАНЕНИЕ И РАСПРЕДЕЛЕНИЕ НЕЙРОННЫХ СИГНАЛОВ: МОДЕЛИРОВАНИЕ 
АКСОННОГО ТРАНСПОРТА 


ФРЭНК РЭТТЕЙ 


Функциональные свойства нейрона могут моделироваться с помощью системы электрических цепей. 
Эта техника пригодна как для естественных процессов, так и для процессов электрического стимули- 
рования. В частности, предложена объединенная модель нейрона млекопитающего, которая учитывает 
различные структуры и особенности нейронной мембраны аксона. Потребляющий большую энергию 
немиелинизированный терминальный участок может быть стимулирован выходным сигналом миели- 
низированного аксона при условии, что последние межузельные расстояния меньше предыдущих. 
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This investigation treats the problem of denaturative effect of temperature on 50 kDa frag- 
ment of myosin 51 in the presence of various nucleosidetriphosphates (NTP) and 
nucleosidediphosphates (NDP). Denaturation changes were evaluated by trypsinolysis extent. 
It is shown that by protecting action NTPs and NDPs may be arranged as follows: 
ATP > CTP > UTP > GTP (1), ADP > GDP > CTP > UDP (2). It is known that by the 
fluorescence increase caused mainly by tryptophanyls of the 50 kDa fragment various NTPs 
may be arranged just as in (1); the same is the sequence of lifetime decrease of longliving 
intermediate complexes M**NDP - Pi. These results suggest that the decrease in the lifetime of 
M**NDP - Pj complex and fluorescence increase are due to the strength of the NTP binding in 
the active centre, resulting in the increase in the 51 rigidity. Especial importance has the 
interaction of the NTP nitrous base and the «hydrophobic pocket» of the ATPase centre. The 
decisive factor of this interaction is probably the NH2-group (proton donor) in position 6 of the 
purine ring of ATP and position 4 of the pyrimidine ring of the CTP, and =N-group (proton 
acceptor) in position 1 and 3, respectively; this stipulates the leading position of ATP and CTP 
in the series (1). In the case of NDP interaction with ATPase site these groups lose their 
importance and hydrophobic forces become more valid. Therefore the series (2) is headed by 
ADP and GDP. Summing up these results we suggest the existence of at least two 
conformational states of the ATPase site «hydrophobic pocket». One of them is responsible for 
decelerating the ATPase reaction on the stage of M**ADP - Рі, Another state takes place at 
actin binding when the structure of «hydrophobic pocket» undergoes reverse changes. 


INTRODUCTION 


Elementary contraction act occurs through the interaction between myosin and 
actin molecules. Fundamental significance in this interaction belongs to myosin head 
(subfragment-1) where ATPase and actin-binding sites are located. Subtragment-1 
(51) consists of one heavy chain and two light ones. The latter seem to play ancillary 
role as SI devoid of them retains the ATPase activity and the ability to actin 
crosslinking [5]. Limited cleavage of the S1 heavy chain yields three fragments (27, 
50 and 20 kDa) [3], arranged along the primary sequence from N- to C-terminal 
respectively. Under non-denaturating conditions all the three fragments are 
interrelated in the structure of such «fragmented» (tryptic) S1, ATPase properties of 
which are practically unaltered while its affinity to actin decreases [4]. More and 
more data are proving that proteolytic fragments of the heavy chain are separate 
structural domains [7]; 27 and 50 kDa domains participate in formation of the 
ATPase site, 50 and 20 kDa domains take part in formation of actin-binding centre. 
The ATPase centre includes the site of polyphosphate binding and.the «hydrophobic 
pocket», responsible for the nitrous base binding [1]. 

51 proteolytic fragments have different sensitivity to high temperatures and 
chemical agents [6]. Thus, at trypsinolysis of 51 preincubated at 35—40° C, 
degradation of the 50 kDa fragment is much more rapid than that of 27 and 20 kDa 
fragments. Degradation of the 50 kDa fragment correlates with the decrease of SI 
ATPase activity. Denaturation of the 50 kDa fragment by denaturating agents (high 
temperatures, methanol) is markedly reduced in the presence of actin, ATP and ADP. 
Alteration of the 51 structural lability in the presence of:specifically interacting agents 
(actin, ATP and ADP) may be of functional significance. Р 

То throw light on the features of myosin ATPase site organization, we have 
studied the process of limited trypsinolysis of SI preincubated at high temperatures in 
presence of various NTPs and NDPs. The results indicate that the structure of the 
ATPase site «hydrophobic pocket» is markedly changed under the influence of 
substrate hydrolysis. 


MATERIALS AND METHODS 


Myosin was isolated from rabbit skeletal muscle [12 ]. Myosin subfragment-1 was 
obtained as described previously [10] with some modifications. 

The purity of protein preparations was checked by SDS-gel electrophoresis, 
spectrophotometry and spectrofluorimetry. The concentration of proteins (mg/ml) was 
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determined spectrophotometrically using relations (Aso — 1.5 А0) /0.533 and 
А зво! 0.75 for myosin and 51 respectively, where А,зо and A329 are absorbtions at 280 
and 320 nm [9]. 

S1 preparations (2 mg/ml) were preincubated for 30 min at 39? C in buffer 
containing (mmol/D: Tris-HCL — 200, pH 7.8; NaN; — 0.5; DTT — 0.5; NTP or 
NDP — 5. After heat treatment S1 preparations were cooled and subjected to limited 
tryptic cleavage at 25? C and 51 to trypsin ratio 100:1 (wt/wt). Trypsinolysis was 
stopped by addition of soy-bean trypsin inhibitor (trypsin/inhibitor=1:2). For 
electrophoretic investigations to aliquots of solutions equal volumes of denaturating 
buffer were added containing (mol/l): НВО; — 0.77; Tris-HCL — 0.2, рН 7.0; 
SDS — 5 75, B-mercaptoethanol — 10 %, sucrose — 40 %. Electrophoresis in 6— 
30 % polyacrylamide gradient gel slabs in presence of SDS was performed as 
described in [13]. Gels were stained with Coomassie Brilliant Blue R-250. After 
washing off the unbound dye the stained gel zones corresponding to the 50 kDa 
fragment were cut out and placed for destaining into solution containing (95): 
isobuthanol — 40, trichloroacetic acid — 10. Optical density of corresponding gel 
bands was estimated by absorbtion Aso, at maximum of absorbtion spectrum of the 
extracted dye with «Specord M-40» (Germany). a-chymotrypsin, TPSK-trypsin, soy- 
bean trypsin inhibitor, NTPs and NDPs were purchased from «Serva». Other reagents 
used were of Soviet production no less than analytical grade. 


RESULTS AND DISCUSSION 


As mentioned above, native structure of the 50 kDa fragment is characterised by 
increased lability and its denaturation changes are already initiated at 35° C and 
higher. Due to the changes this fragment becomes more sensitive to trypsinolysis the 
extent of which allows to determine the degree of its denaturation. The level of 
trypsinolysis may be estimated by the density of the 50 kDa band in 51 
electrophoregram: the.less is the band density the higher are the severity of 
trypsinolysis of the 50 kDa fragment and denaturation changes after heat treatment. ' 
Heat treatment (39* C) denaturation of the S1 in the presence of various NTPs was 
studied (Fig. 1). 
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Fig. I. a) Electrophoregrams of 51 preparations, preincubated in the presence of various NTP at 39° C and 
after cooling subjected to limit tryptic cleavage at 25° C. Experimental conditions are described in the text. 

b) Optical absorption at 590 nm of the dye extracted from gel bands corresponding to 50 kDa frag- 
ment of 51 on the electrophoregrams. On abscissa the following additions are shown: 1-АТР, 2-GTP, 3-UTP, 
4-СТР, 5-РРРІ, 6-151 control (in the absence of NTP, NDP, PPPi and PPi), 7-2nd control (the same as the 
first control, but the sample wasn’t preincubated at 39° C and was stored at 25° C). 
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All NTPs exert protective effect, i. є. at their addition into reaction mixture 
denaturative changes of S1 markedly decrease. Preceeding from the density of the 50 
kDa fragment electrophoretic bands and the intensity of absorbtion of the dye 
extracted from these bands (Fig. 1b), NTPs might be arranged by their protective 
effect as follows: 


ATP > CTP > UTP > GTP. (1) 


It is known that nucleosidetriphosphate hydrolysis in the presence of Ме?" 
includes several ый [2]: 


3 4 5 
М+МТР= М. ТР М". NTP > М", МОР. Ріе М" МОР : РІ = 
lw NDP4PleM- NDP 2 M + NDP (2) 


where M* and M** are conformational states of myosin or its active subfragments 
with increased levels of their intrinsic fluorescence. Stage 4 limits the rate of the NTP 
hydrolysis. At the excess of substrate the long-living complex M** - МОР - РІ 
accumulates in the solution. It is shown that by the degree of increase in intrinsic 
fluorescence, due mainly to the 50 kDa tryptophanyls, various NTPs can be arranged 
just as in (1). In the same sequence there decreases the life-time of long-living 
intermediate complex M** - NDP - P1 which varies inversely to the rate of various 
NTPs hydrolysis by myosin and its active fragments [15]. 

All these results may be explained proceeding from the fact that at NTP binding 
the S1 structure becomes denser and more rigid. Responsible for intrinsic fluorescence 
tryptophanyls on the condition of such 51 structure should have higher quantum yield 
manifesting itself by an increase in 51 fluorescence after NTP addition. Such 
structure should also be more resistant to denaturating factors including high 
temperatures; this observation is consistent with both our and other worker's results 
[6]. These structural transitions are probably connected with the most labile SI 
region, namely the 50 kDa fragment. 
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Fig. 2. a) Electrophoregrams of 51 preparations preincubated in the presence of various NDP at 39" С апа 
after cooling subjected to limited tryptic cleavage at 25° C. 
b) Optical absorbtions at 590 nm of the dye extracted from gel bands corresponding to 50 kDa frag- 
ment of 51 on the electrophoregrams. 
On abscissa the following additions are shown: 1-ADP, 2-GDP, 3-UDP, 4-CDP, 5-PPi, 6-ATP, 
7-control (additives are absent). 
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In presence of NDP in incubation mixture denaturation changes of S1 are also 
considerably decreased (Fig. 2). But in contrast to NTP, nucleosidephosphates may 
be arranged, according to their protective effect, as follows: 


ADP » GDP » CDP » UDP. (3) 


Figs. 1, 2 show the results of the trypsinolysis of SI preincubated at 39° C in 
presence of triphosphate or pyrophosphate and in absence of these agents or any 
other nucleotides. Fig. la, b also shows the results of the control experiment where 51 
was subjected to trypsinolysis without heat pretreatment in absence of nucleotides. 
Triphosphate and pyrophosphate also protect the native structure of SI from heat 
denaturation. In absence of NTP, NDP, PPPi or PPi in the solution heat 
denaturation of 51 is much rapider, so the 50 kDa band on the electrophoregram 
disappears after trypsinolysis. 

From the results presented above, it is suggested that the increases in life-time of 
the long-living complex M** - NDP - РІ and its fluorescence intensity result from the 
stability of NTP binding to the active centre which leads to the enhancement of 
density and rigidity of 51 structure. The interaction between the nitrous base of NTP 
and the ATPase «hydrophobic pocket» is of prime importance. The key role here 
might be played by the NH-group (proton donor) in position 6 of the ATP purine ring 
and corresponding position 4 of the CTP pyrimidine ring, and =N-group (proton ac- 
ceptor in position 1 and 3, respectively (Fig. 3). This provides the leading position of 
ATP and CTP in series (1). It must be noted that the above mentioned groups are 
responsible for formajion of hydrogenous bonds between complementary bases in the 
double-stranded DNA [14]. 


RH. 
2 ИН 
| i A R ! ) 
276 TN ae 
li У Е 
4th Е | di 
3 9 ly £i 
"SN Р 5/ в ae 
их ale 
| АДЕНИН Н 1 цитозин к x 
| 0 
| 0 
| М. Nv ; Д | 
| i ] 
E | | | 
№: XN М п 
З H.1 у чит < 
ГУАНИН NS УРАЦИЛ ( H i Fig. 3. Structural formula of purine and 


pyrimidine bases. Circles correspond to 
hydrogene atoms released at binding to 
shugar residue. 


- 


In cases where NDPs interact. with the ATPase site forming M* - NDP complex 
these groups become invalid and hydrophobic interactions gain more importance. 
Therefore the series (3) is opened with ADP and GDP containing nitrous bases with 
solid purine ring (Fig. 3), possessing pronounced hydrophobic features. The evidence 
to this is in the fact that series (3) is exactly the same as the sequence of retention of 
adenosine, guanosine, cytidine and uridine on ion-exchange chromatographic column 
filled with substance possessing pronounced hydrophobic properties [11]. Due to 
hydrophobic interactions with the column bed adenosine and guanosine are eluted 
much later than cytidine and uridine. 

Goodno et al. [8] reported that at interaction of ADP with myosin in the 
presence of VO}~-ions the complex M** - ADP - УО; is formed similar to complex 
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M** - ADP - Р,. This might be the case in the presence of other NDPs in the reaction 
mixture. Therefore one may expect that in the presence of VO}~ protective effect of 
NDP would be analogous to that of NTP. Indeed, as we have found (see Fig. 4), in 
this case NDPs may be arranged by their protective action as follows: 


ADP > CDP > UDP > GDP. (4) 


This series is identical to series (1). 
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Fig. 4. a) Electrophoregrams of 51 preparations preincubated in the presence of various NDP and 5 mM 
УОҘ at 39° C and after cooling subjected to limited tryptic cleavage at 25° C. 

b) Optical absorption at 590 nm of the dye extracted from gel bands corresponding to 50 kDa frag- 
ment of S1 on the electrophoregrams. 

On abscissa the following additions are shown: 1-ADP, 2-GDP, 3-UDP, 4-CDP, 5-control (additives 
are absent). 


Thus, binding of the nucleotide base to the «hydrophobic pocket» of the myosin 
ATPase site results in additional stabilization of the 50 kDa fragment and diminishes 
the effect of denaturating factors. Maximal stability is characteristic for the structure 
of intermediate complex M** - ADP - Рі; at its formation essential changes may occur 
both in the «hydrophobic pocket» and the 50 kDa fragment on the whole. During ATP 
hydrolysis ATPase site of myosin provides some conformations differing by the 
pattern of interaction with substrate and by reaction products. Myosin conformations 
in intermediate complexes М" - АТР, М" - ADP - P, and М" - ADP, forming on 
stages 2, 4 and 5 (scheme 2), are possibly identical. On binding of the whole molecule 
of substrate or product to the ATPase site polyphosphate and rybose groups must also 
have importance (as may be seen in Fig. 1, 2, PPP; and PP, also have protective 
effect). But the effect of these groups must be the same for all NTPs (NDPs). 

These data suggest that at interaction of ADP with the ATPase site the nitrous 
base binds to the «hydrophobic pocket» primarily due to weak hydrophobic forces. 
Probably it may be true for ATP up to stage 2 (scheme 2). Such interaction is 
accompanied by a small increase in fluorescence. At the moment of macroergic bond 
splitting (stage 3) a conformational transition of the «hydrophobic pocket» occurs 
resulting in its accommodation to the spatial structure of ATP ring and development 
of polar interaction between -NH- and =N-groups of the nitrous ‘Базе. and 
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complementary groups of the «pocket». As a result the base forms rather strong bond 
with the «hydrophobic pocket» and gives rise to the long-living complex 
M** - ADP - Р; characteristic by maximal increase of its own fluorescence. The 
strength of the cytidine base binding in the «hydrophobic pocket» at CTP hydrolysis 
must be the same because spatial arrangement of -NH- and =N-groups in cytidine is 
analogous to that in adenine. As to the UTP and GTP, binding of their bases in the 
«pocket» occurs due to hydrophobic forces despite conformational transitions of the 
«pocket» at hydrolysis of these nucleosidetriphosphates because in positions of their 
nitrous bases analogous to positions 6 and 1 of adenine the groups witn the opposite 
polarity are located (Fig. 3). Therefore complementary interactions with polar groups 
of the «hydrophobic pocket» are impossible for them. 

lt is known that at actin binding to the myosin head the ATPase reaction is 
promoted more than by an order of magnitude. The kernel of regulatory effect of 
actin on the ATPase site must consist in reverse changes of the «hydrophobic pocket» 
at binding of SI. As a result polar interactions of complementary groups with the 
«pocket» break down, the bond of the nitrous base with the «pocket» weakens 
abruptly and the macroergic energy accumulated in 51 structure is utilized for 
clementary contraction act. 


ФУНКЦІОНАЛЬНО РІЗНІ СТАНИ «ГІДРОФОБНОГО KAPMAHY» АТФазного ЦЕНТРУ 
МІОЗИНУ | У 


A. M. ФІЛЕНКО, B. М. ДАНИЛОВА, Е. Б. БАБІЙЧУК 


Досліджено денатураційний вплив підвищеної температури (39° С) на фрагмент 50 кДа субфрагменту 
І міозину в присутності різних нуклеозидтрифосфатів (НТФ) та нуклеозиддифосфатів (НДФ). Ступінь 
денатураційних змін цього фрагменту оцінювали за глибиною його трипсинолізу. Показано, що за 
протектуючою дісю НТФ та НДФ розташовуються таким чином: АТФ > ЦТФ > УТФ > ГТФ та 
АДФ > ГДФ > ЦДФ > УДФ, відповідно. Такий порядок розташування різних НТФ та НДФ обумов- 
лений ступенем зв'язування їх основ з «гідрофобним карманом» АТФазного центру міозину. Зв'язуван- 
ня основи приводить до додаткової стабілізації структури фрагменту 50 кДа та меншого впливу на неї 
денатуруючих факторів. Одержані результати разом з літературними даними свідчать про те, що є два 
стани «гідрофобного карману» АТФазного центру, один з яких відповідає за уповільнення АТФазної 
реакції на стадії довготривалого комплексу M** - АДФ - Фу. Сутність регуляторного впливу актину на 
АТФазний центр є, певно, в тому, що при зв'язуванні актину міозином «гідрофобний карман» актив- 
ного центру зазнає конформаційних змін, які приводять до послаблення зв'язування основи 3 цим цен- 
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тром i, відповідно, до розпаду комплексу М” · АДФ · Фу. 


ФУНКЦИОНАЛЬНО РАЗЛИЧНЫЕ СОСТОЯНИЯ «ГИДРОФОБНОГО КАРМАНА» АТФазного 
ЦЕНТРА МИОЗИНА 


А. М. ФИЛЕНКО, В. М. ДАНИЛОВА, Э. Б. БАБИЙЧУК 


Исследовали денатурационное влияние повышенной температуры (39° С) на фрагмент 50 кДа суб- 
фрагмента 1 миозина в присутствии различных нуклеозидтрифосфатов (НТФ) и нуклеозиддифосфа- 
тов (НДФ). Степень денатурационных изменений этого фрагмента оценивали по глубине его трипси- 
нолиза. Показано, что по протектирующему воздействию НТФ и НДФ располагаются следующим об- 
разом: АТФ > ЦТФ > УТФ > ГТФ и АДФ > ГДФ > ЦДФ > УДФ, соответственно. Такой порядок 
расположения различных НТФ и НДФ обусловлен степенью связывания их оснований с «гидрофоб- 
ным карманом» АТФазного центра миозина. Связывание основания приводит к дополнительной стаби- 
лизации структуры фрагмента 50 кДа и меньшей подверженности ее влиянию денатурирующих фак- 
торов. Полученные результаты совместно с литературными данными свидетельствуют о том, что име- 
ется по крайней мере два состояния «гидрофобного кармана» АТФазного центра, одно из которых 
отвечает за замедление АТФазной реакции на стадии долгоживущего комплекса м** . АДФ + Фи. 
Сущность регуляторного воздействия актина на АТФазный центр состоит по-видимому в том, что при 
связывании миозином актина «гидрофобный карман» активного центра претерпевает конформацион- 
ные изменения, приводящие к ослаблению связывания основания с этим центром и, соответственно, к 
распаду комплекса M** · АДФ · Фи. | 
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ВПЛИВ pH І БІЛКОВОГО ОТОЧЕННЯ НА ФЛУОРЕСЦЕНТНІ 
ПАРАМЕТРИ 1НДО-1 
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THE ACTION OF THE pH AND PROTEIN SURROUNDING ON THE FLUORESCENT 
PARAMETERS OF INDO-1 


V. L. ZYMA, O. M. DYACHOR, A. P. VETROV 


It was established that in non calcium medium (2mM EGTA) with pH < 8.0 the integral 
fluorescent spectrum of indo-1 consists of three defferent spectra of individual forms: calcium 
free unprotonated L-form (тах = 477 nm), partially protonated LH1- (max = 446 nm) and 
protonated LH2-form (max = 412 nm). Acidosis leads to change in percentage relation of these 
forms. At the pH5.68 the LH2-form is dominated (63 %) in the integral fluorescent spectrum 
by other forms. This predominance initiates a marked increase in calibrated parameter Rmin 
and corresponding errors in intracellular calcium [Са2+]; registration. In calcium saturated 
mediums ({Ca2*] = 2 mM) with pH from pH9.0 to pH5.2 the whole spectrum of indo-1 and 
calibrated parameter Rmax is altered insignificantly. It is shown that in medium with serum 
albumin (pH > 8.0) besides spectra of L-form, the LP-form fluorescent spectrum of indo-1 
with Amax = 445—450 nm is detected. This LP-form is due to binding indo-1 with serum 
albumin. The red shift of Amax fluorescent spectrum of indo-1-protein complex upon the 
exitation at the long wavelength absorbtion edge is established. The red shift testifies that the 
shift of indo-1 fluorescent Amax toward short wavelenghts is caused not only by hydrophobic 
microsurrounding but also by decreasing in protein dynamics in the indo-1 binding site. 


ВСТУП 


Флуоресцентний зонд індо-1 широко використовується для визначення KOH- 
центрації вільного кальцію [Ca** | в клітинах, але кількісні розрахунки [Ca?* ]; 
обмежені певними ефектами внутрішньоклітинного середовища при зв'язуванні 
іонів кальцію індо-1 [1—7 ]. В запропонованому методі [1] [Ca?* | знаходиться за 
умови, що 1ндо-1 в розчині перебуває в двох формах: в кальцій- вільній /-формі 
i кальцій-зв'язаній 1М-формі. Чутливість методу обумовлена важливою вла- 
стивістю індо-1 зміщувати максимум спектру флуоресценції від 477—485 нм (у 
відсутності іонів кальцію, 2 MM ЕГТА) до 415 нм (при зв'язуванні іонів кальцію 
1ндо-1) [1, 3, 8, 9]. Реєструється двохвильовий параметр А=[.10//лоз, де [410 — 
інтенсивність флуоресценції ГМ-форми при 410 нм та Z495 — інтенсивність флуо- 
ресценції Г-форми індо-1 при 495 нм. 

Було відмічено, що спектр флуоресценції індо-1 в клітинах значно 
відрізняється від спектру флуоресценції індо-1 в буферних розчинах з малою 
концентрацією кальцію [10]. Аналіз показав, що зміщення в короткохвильову 
ділянку спектру флуоресценції індо-1 в клітинах обумовлено не тільки за paxy- 
нок зв'язування внутрішньоклітинного кальцію, але певний зсув спектру флуо- 
ресценції може спричиняти зв'язування індо-1 з цитоплазматичними білками [8, 
11]. Відомо також, що спектр флуоресценції індо-1 і константа дисоціаціїї Ka 
змінюються при закисленні середовища [12]. В дослідах in vitro (моделювалась 
зміна рН і додавались різні білки) виявлено існування 4-х індивідуальних форм 
індо-1: кальцій-вільна депротонована Г-форма при pH > 8.0 (Amas = 477 нм [9], 
Алах = 485 нм [8 D, кальцій-зв'язана |М-форма Amas = 415 нм), протонована LH- 
форма Umax = 455 нм) при pH < 8.0 і ІР-форма (Amas = 438 нм), яка утворюється 
за рахунок зв'язування 1ндо-1 з деякими білками [8]. Крім того було встановле- 
HO, що спектр флуоресценції LP-cbopmu індо-1 представляє собою комбінацію двох 
спектральних форм: ІР-форми з А.х = 445 нм, яка обумовлена зв'язуванням 
індо-1 в гідрофобних «кишенях» білків (виявлено в розчинах сивороткових аль- 
бумінів) i LP*-cbopma з Amas = 454 нм, яка обумовлена електростатичною взаємо- 
дією індо-1 з позитивно зарядженими залишками лізінів і аргінінів в білках (ви- 
явлено в розчинах гістонів та полі-І-лізіну) [11]. Використовуючи значення 
квантових виходів і спектри флуоресценції індивідуальних форм індо-1, розробле- 
но метод комп'ютерного розкладання складних спектрів флуоресценції індо-1 в 


74 


Вплив pH і білкового оточення на флуоресцентні параметри індо-1 


інтактних клітинах, що дозволяє оцінювати внесок різних форм в інтегральний 
спектр і більш точно визначати концентрації цитозольного кальцію [Ca?* |; [8]. 

Незважаючи на певні успіхи, ще не зовсім з'ясовано зміни спектральних ха- 
рактеристик індо-1 і механізм зв'язування 1ндо-[ з білками. Інтерес викликає 
вплив закислення середовища на флуоресценцію індо-1, оскільки відомо, що ряд 
фізіологічних процесів супроводжуються зміною pH внутрі клітини, яке 
реєструється рН-чутливими флуоресцентними зондами [5]. З метою більш точно- 
го врахування внеску некальцієвих форм 1ндо-1 в інтегральні слектри флуорес- 
ценції індо-1, які реєструються в клітинах, B цій роботі досліджено флуоресцентні 
характеристики індо-1 при різних рН розчину і вивчено зв'язування індо-1 з си- 
воротковими альбумінами. 


МАТЕРІАЛИ І МЕТОДИ 


Спектри флуоресценції реєстрували в | см кварцевій кюветі на спектро- 
флуориметрі СДЛ-2 («ЛОМО», Росія). Флуоресценцію 1ндо-1 збуджували на 
А = 350 нм. Всі вимірювання проводились при кімнатній температурі. 
Вимірювання pH розчинів проводилось на рН-метрі (pH-150). 

В роботі використовувались індо-1 («Sigma», США), сиворотковий альбумін 
бика («Fluca Chemical», Швейцарія) і сиворотковий альбумін людини («Reanal», 
Угорщина). Всі розчини індо-І та білків готувались на бідистильованій воді 3 пи- 
томим опором 2 МОм · см. 


РЕЗУЛЬТАТИ ТА ОБГОВОРЕННЯ 


Внутрішньоклітинна концентрація вільного кальцію визначається за форму- 
лою (при умові, що в цитоплазмі є тільки L- та LM-cbopma): 


(Car l = Ка 7 В і (R Ha Ras / (Rmax T R), (D 


ne А=1,10/ 145 — відношення інтенсивностей флуоресценції 1ндо-1 на 410 i 495 нм; 
В = 14о5/ Гьдоѕ — відношення інтенсивностей флуоресценції І-форми до .М-форми 
на 495 нм; Ку — константа дисоціації комплексу Са-індо-1; Rmin — ДВОХВИЛЬОВИЙ 
флуоресцентний параметр 1ндо-! у відсутності кальцію (B розчині 2—5 мМ 
ETA): Rmax — двохвильовий параметр при наявності насичуючих концентрацій 
кальцію, 2 мМ И]. Правильне визначення [Ca?* |, вимагає точної реєстрації 
калібровочних параметрів А, Ry, В i Ка. Але виявляється, що їх значення 
суттєво залежать від зміни температури, іонної сили, pH, білкового оточення. 

Ми дослідили вплив закислення середовища на флуоресценцію 1ндо-1. B без- 
кальцієвому розчині (2 MM ETTA) при зміні pH від рНІ0.0 до pH8.0 реєструється 
депротонована кальцій-вільна L-dopMa (Amas = 477 нм) індо-1, спектр флуорес- 
ценції якої представлено на рис. 1 (спектр 1). Зменшення рН < 8.0 викликало 
зміщення максимуму спектру флуоресценції індо-і в короткохвильову сторону по 
відношенню до спектру індо-1 при pH10.0 i для розчинів індо-1 з pH7.0 та pH6.53 
Алах були 469 нм і 460 нм, відповідно (рис. l, спектри 2 i 3). Значні зміни в 
спектрі флуоресценції індо-1 спостерігались при закисленні середовища до pH5.5. 
Як видно з puc. 2, y безкальцієвому розчині з рН5.68 спектр флуоресценції індо- І 
значно зміщений (ух = 420 нм) і розширений за рахунок значного зростання 
короткохвильового краю спектру. Зміна інтегрального снектру флуоресценції 
індо-1 при закисленні розчину від pH8.0 до рН5.2 пояснювалась [8] перероз- 
поділом спектрів флуоресценції двох форм індо-1 — кальцій-вільної І форми 
max = 485 нм) і протонованої ІН-форми Й пах = 455 нм), рівновага яких характе- 
ризується pK, = 6.18. Але встановлений нами характер змін спектру флуорес- 
ценції індо-1 при рН < 8.0 відмінний від даних роботи [8 |. Це навело нас на дум- 
ку, шо інтегральний спектр флуоресценції компонується більш складно i3 
індивідуальних форм 1ндо-1. 
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Рис. І. Комп'ютерне розкладання інтегральних спектрів 1ндо-! (5 MKM) в розчинах з рНІ0.0 (1), 
pH7.0 (2), pll6.53 (3) на спектри Е-форми Qa, За), ІЛІї-форми (26, 36) i 1Н2-форми (26, 36). 


736 = 350 нм; температура 20° С. Ip -- інтенсивність флуоресценції в відносних одиницях. 


Для перевірки висловлених припущень ми провели комп'ютерне розкладан- 
ня інтегральних спектрів флуоресценції індо-1 на спектри індивідуальних форм. 
Прийнявши спектр флуоресценції індо-1 при pH10.0 як індивідуальний (конт- 
рольний) спектр Ё-форми, ми провели нормування на довгохвильовому крилі 
цього спектру із інтегральним спектром 1ндо- при pH7.0, де вони добре 
співпадають (рис. 1, спектр 2a). Після вирахування спектру 2a із інтегрального 
спектру 2 ми отримали спектр флуоресценції, який не співпадав 3 спектром 
| Н-форми в роботі [8] i який мав більш складну форму. Труднощі з його розкла- 
данням на складові спектри було знято при аналізі інтегрального спектру флуо- 
ресценції індо-І при pH5.68 (рис. 2, спектр 2). Виявилось, що він добре розкла- 
дається на два контрольні спектри — спектр І-форми (рис. 2, спектр 2a) і новий 
спектр з Ag = 412 нм (рис. 2, спектр 2в), який ми віднесли до протонованої 
ІН2-форми індо-1. Це дозволило нам розкласти інтегральні спектри флуорес- 
ценції індо-1 в розчинах 3 РН < 8.0 на три контрольні спектри 1-форми 
(Дах = 477 нм), частково-протонованої ІНІ-форми (ах = 446 нм) і протонованої 
ІН2-форми (Amas = 412 нм) (рис. 1 і 2). Ці три форми дають різний процентний 
внесок в інтегральний спектр флуоресценції індо-1 при різних pH: рН10.0 — 
100 % І-форми; pH7.0 — 67 У, (D, 21 % ЯНЬ, 12 У, (LH2); pH6.53 — 54 % 
(L), 21 % ЫН), 25 % (LH2); pH6.2 — 44 % (D, 12 % ЯНЬ, 44 % (LH2); 
pH5.68 — 36 % (Ш, 1 % (LH1), 63 % (LH2). 
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Рис. 2. Комп'ютерне розкладання інтегральних спектрів 1ндо-! (5 мкМ) в розчинах з pH6.2 (1), 
рН5.68 (2) на спектри і-форми (la, 2a), МИ-форми (16, 26) i 192-форми (Ів, 28). 456 = 350 нм; 
температура 20° С. 
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Спектр флуоресценції .Н2-форми Ama: = 412 нм) близький до спектру флуо- 
ресценції 1М-форми (Amas = 415 нм). Ми не виключали певного внеску в 
інтегральний спектр індо-1 при рН5.68 1М-форми, яка утворюється завдяки 
зв'язуванню молекулами індо-1 невеликих кількостей іонів кальцію в бідистиляті 
(за даними роботи [12] загальна слідова концентрація іонів кальцію біля 
3-107? моль/л). Концентрація вільного кальцію [Ca?* |, при заданій концент- 
рації ЕГТА буде залежати від рН розчину, оскільки змінюється константа ди- 
соціації (чи константа асоціації) комплексу Са-ЕГТА. Розрахунок констант 
асоціації К'єзргта при різних рН проводився за формулою [13]: 


Kear 


d PEE - - Fe РН + 
K'cwrr = 34 (Ht -K, + [H* P-K K; + [Ht PK; Ka K; + [H* J- Ki: Ka: Ks K4 
i Kea2 2 
UIH*FK)*1*IH FK; + НР. К, K; + [HFP KK; K)? O 

де Ki, К, Кз i Ky — константи асоціації для окремих протонів; Кса i Ke; — 


індивідуальні константи асоціації для іонів кальцію. При температурі 20° C ta 
іонній силі 0.1 М використані такі значення констант K, = 9.47; К, = 8.85; 
Кз = 2.66; К. = 2.0; Kea = 10.97; Ке = 5.3. Розрахувавши за формулою (2) 
К'єавгта, МИ Знайшли як змінюється концентрація вільного кальцію [Ca?* |, при 
наявності в розчині слідової концентрації кальцію 3: 1079 моль/л i ETTA 2 ммоль/л: 
при рНІ0.0 — 1.7 : 10-15 моль/л; pH7.0 — 3.57 110719; pH6.53 -- 3.1:107°; 
pH6.2 — 1.4- 1078; pH5.68 — 1.4- 1077 моль/л. Отримавши [Ca?* | = 1.4: 1077 моль/л 
при pH5.68, можна оцінити за формулою (1) величину R і порівняти її з експери- 
ментальним значенням В, отриманим із інтегрального спектру флуоресценції 
1ндо-1 при рН5.68. В формулу (1) підставляли такі калібровочні параметри, отри- 
мані в дослідах in vitro: Rmin = 0.1; Rmax = 2.0; В = 1.18. Константу дисоціації 
К, = 3160 нМ при рН5.2 взяли із роботи [12]. Розрахунки за формулою (1) дали 
значення R - 0.1, в той час, як експериментально виміряна величина дає К = 1.62 
(рис. 2, спектр 2). Це свідчить, що в розчині з 2 MM ETTA участь LM-dopmu не 
суттєва і ссновний внесок в інтегральний спектр флуоресценції при pH5.68 дає 


Lio 


спектр 7) 2-форма 1ндо-1. 
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Рис. 3. Залежність двохвильового флуо- 

ресцентного параметру індо-1 A (1, 2) та 45 
параметру Py (співвідношення LH2- до 

ІНІ-форми) (3) від рН розчину. 1 i 3 — 0 
у відсутності [Са2 7] (2 MM ЕГТА); 2 — кон- 


центрація [Са2+] = 2 мМ. 456 = 350 нм; 5 6 7 8 9 10 
температура 20° С. 
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В безкальцієвому розчині (2 мМ ЕГТА) двохвильовий параметр А, який є 
Rmin У формулі (1) [1], при зміні рН від pH10.0 до рН8.0 залишається постійним 
i дорівнює Rmin = 0.1. При зменшенні рН < 8.0 Rmin швидко зростає i при pH5.54 
він досягає значення 1.77 (рис. 3, крива /). Ми порівняли тенденцію зміни napa- 
метру Rmin, отриманого із інтегральних спектрів 1ндо-1, з поведінкою параметру 
Pa = 1412/1446, який характеризує співвідношення LH2-cdopmu до ІНІ-форми 
індо-1. При зниженні pH від pH8.0 до рН7.0 параметр Р, мало змінюється, 
оскільки ще невеликий процентний внесок LH2-cbopmu, i в цих межах рН незнач- 
но зростає Rmin Подальше закислення розчину (рН < 7.0) параметр Р, різко зро- 
стає завдяки переважному утворенню 1Н2-форми 1ндо-1. З цим добре корелює 
значне зростання Rmin Тому використання Rmin із інтегрального спектру індо-1 в 
розчинах з рН < 8.0, коли не враховується процентний внесок LH2- i ІНІ-форм, 
приведе до значних похибок у визначенні [Ca?* |; за формулою (D. 

Кількісне визначення [Ca?' |; обмежується певними ефектами внутрішньо- 
клітинного середовища при зв'язуванні іонів Са?" з індо-1. Установлено, що 
максимум спектру флуоресценції індо-І в клітинах зсунутий в голубу ділянку 
спектру порівняно з спектром флуоресценції індо-1 в бозбілкових розчинах при 
нейтральних pH [10]. Такий зсув спектру флуоресценції індо-і зумовлюють 
деякі білки, зокрема сивороткові альбуміни, гістони |8, 9], які зв'язують 1ндо-1, 
утворюючи флуорссцюючу [Р-форму 1ндо-1. Нами виявлено, що вирішальну 
роль в зв'язуванні індо-1 в неполярних «кишенях» сивороткових альбумінів 
відіграють гідрофобні взасмодй i неполярне мікрооточення обумовлюс коротко- 
хвильове положення A4, = 445—450 нм. Частково в стабілізації комплексу 
відіграють електростатичні взаємодії індо-1 з позитивно заряженими лізіновими і 
аргініновими залишками амінокислот. Короткохвильове положення Ах індо-І 
при зв'язуванні з сиворотковим альбуміном також може бути обумовлено 3MCH- 
шенням рухомості білкових груп в ділянці зв'язування індо-1, як це було виявлс- 
но для таких флуоресцентних зондів як AHC, THC, ДНС у в'язких розчинах або 
в твердих матрицях [14]. Експериментально установлено, що в системах із 
сповільненою релаксацією мікрооточення навколо флуоресцентного зонда спо- 
стерігається зсув А. флуоресценції в довгохвильову сторону при «червоному» 
крайовому збудженні. На рис. 4 зображена залежність максимуму спектру флуо- 
ресценції Amas індо-1 від довжини хвилі збуджуючого світла Аб в воді (крива /) i 
при зв'язуванні з сиворотковим альбуміном (крива 2). У воді ефект червоного 
краю не спостерігасться, шо означас швидку рухомість полярних молекул води 
навколо індо-І. При зв'язуванні індо-1 з сиворотковим альбуміном CNO- 
стерігасться сфект червоного крайового збудження — при зміні 4,5 від 370 нм до 
400 нм Ag, індо-Ї зсувається в довгохвильову сторону від 450 нм до 458 нм (рис. 
4, крива 2). Такий спектральний зсув А.а флуоресценції при червоно-крайовому 
збудженні вказус на те, що крім гідрофобності середовища певну роль відіграє 
зростання жорсткості мікрооточення при зв'язуванні індо-1 з білком. Висловлю- 
сться сумнів, що при зв'язуванні іонів Са?" з 1ндо-|, коли утворюється LM-dop- 
ма, ця форма HC взаємодіє з сиворотковим альбуміном [8]. Ми показали, що 
1.М-форма взаємодіє з сиворотковим альбуміном, утворюючи [МР-комплекс, npo 
що свідчить зміщення Ag, Від 415 нм (І.М-форма) до 412 нм ((МР-форма). 

Наявність флуоресцюючої ІР-форми індо-1 в розчині з pH > 8.0, коли 
відсутні LH1- i 1Н2-форми, змінює калЮровочний параметр Rmin і він зростає від 
0.1 (безбілковий розчин) до 0.58, коли додасться 50 мкМ сивороткового аль- 
буміну. Rmax і f незначно змінюються. Білкове оточення виявляє вплив на Ка B 
безбілковому розчині з pH7.3 К = 210 нМ, а при додаванні білкового екстракту із 
м'язів жаби Ky збільшується до 808 HM [15]. 

Таким чином, проведен! дослідження свідчать, що взаємодія індо-1 з прото- 
нами або білковими молекулами, зокрема сиворотковими альбумінами, змінює 
зв'язування іонів Са?" з індо-І та флуоресцентні параметри індо-1. Все це 
потрібно враховувати при визначенні калібровочних параметрів, які використову- 
ються для визначення внутрішньоклітинного вільного кальцію |Са?" |. 
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Puc. 4. Залежність максимуму спектра флуоресценції (max) індо-1 (5 MKM) від довжини хвилі збуджу- 
ючого світла 436: Í — у воді; 2 — в комплексі с сиворотковим альбуміном (50 мкМ). Температура 
20° C; pH7.0. 


ВПЛИВ pH І БІЛКОВОГО ОТОЧЕННЯ HA ФЛУОРЕСЦЕНТНІ ПАРАМЕТРИ 1НДО-1 
В. Л. ЗИМА, О. M. ДЯЧОК, О. IT. ВЕТРОВ 


В роботі установлено, що в безкальцієвих розчинах (2 мМ ЕГТА) з рН < 8.0 інтегральний спектр флу- 
ореспенції індо-1 розкладається на спектри трьох індивідуальних форм — кальцій-вільної депротоно- 
ваної [-форми Umax = 477 нм), частково протонованої LH1- (max = 446 нм) і протонованої LH2-dopmu 
(max = 412 нм). Закислення середовища змінює процентне співвідношення цих форм і при pH5.68 в 
інтегральному спектрі переважає LH2-cdopma (63 %), що спричинює значне зростання калібровочного 
двохвильового параметра Rmin і відповідні похибки в визначенні концентрацій внутрішньоклітинного 
вільного кальцію [Ca?* ]i. B кальцій-насичуючих розчинах ({Ca? 7] = 2 MM) інтегральний спектр та 
калібровочний параметр Rmax майже не змінюються при зміні рН від рН9.0 до рН5.2. Виявлено, що в 
розчинах з сиворотковим альбуміном (pH > 8.0) крім спектру І форми реєструється спектр флуорес- 
ценції LP-cbopmu 1ндо-1 з Amax = 445—450 нм, яка обумовлена зв'язуванням індо-1 з сиворотковим аль- 
буміном. Виявлений ефект червоного крайового зсуву Amax спектру флуоресценції комплексу 1ндо-1- 
білок свідчить про те, що не тільки гідрофобність мікрооточення обумовлює короткохвильовий зсув 
Amax індо-1, але певну роль відіграє зменшення рухомості білкових груп в місці зв'язування індо-1. 


ВЛИЯНИЕ pH И БЕЛКОВОГО ОКРУЖЕНИЯ HA ФЛУОРЕСЦЕНТНЫЕ ПАРАМЕТРЫ ИНДО-1 
В. Л. ЗИМА, О. М. ДЯЧОК, A. П. ВЕТРОВ 


Установлено, что в бескальциевых растворах (2 MM ЗГТА) с pH < 8.0 интегральный спектр флуорес- 
ценции индо-1 состоит из спектров трех индивидуальных форм — кальций-свободной депротониро- 
ванной 1-Формьт (тах = 477 нм), частично протонированной LH 1 - (max = 446 нм) и протонированной 
LH2-dopMbt. (max = 412 нм). Закисление среды изменяет процентное соотношение этих форм и при 
pH5.62 в интегральном спектре флуоресценции преобладает LH2-copwa (63 %), что вызывает значи- 
тельное увеличение калибровочного параметра Amin и соответствующие ошибки в определении кон- 
центрации внутриклеточного свободного кальция [Cat]. В кальций-насыщающих растворах 
([Ca^*] = 2 мМ) интегральный спектр и калибровочный параметр Rmax изменяются незначительно 
при изменении pH от pH9.0 до рН5.2. Показано, что в растворах C сывороточным альбумином 
(pH > 8.0) кроме спектра Г-формы регистрируется спектр флуоресценции ГР-формы индо-1 c 
Атах = 445—450 нм, которая обусловлена связыванием индо-1 с сывороточным альбумином. Установ- 
ленный эффект красного краевого сдвига Amax флуоресценции комплекса индо-1-белок свидетельству- 
ет, что не только гидрофобность микроокружения обуславливает коротковолновой сдвиг Amax индо-1, 
но определенное влияние оказывает уменьшение подвижности белковых групп в области связывания 
индо-1. 
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MEDICAL AND BIOLOGICAL APPLICATIONS 


НОРМАЛІЗУЮЧИЙ ВПЛИВ СТОХАСТИЧНОГО УЛЬТРАЗВУКУ 
НА ХЕМІЛЮМІНЕСЦЕНЦІЮ ПЛАЗМИ КРОВІ ТА АТФ-азну 
АКТИВНІСТЬ ЕРИТРОЦИТІВ КРОВІ ХВОРИХ НА 
ХРОНІЧНИЙ ТОНЗИЛІТ 
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NORMALIZING INFLUENCE OF STOCHASTIC ULTRASOUND ON CHEMILUMINESCENCE 
OF BLOOD PLASMA AND ATP-ase ACTIVITY ОЕ THE BLOOD ERYTHROCYTES OF 
PATIENTS WITH CHRONIC TONZILLITE 


KARAS A. F., ANDREICHENKO S. V., NURISCHENKO N. E., STRYZHAK У. V. 
and OSTROVSKY I. V. 


Institute of otolaringology, Kiev 
National University of Kiev-Mohyla Academy, 
T. G. Shevchenko National University 


Abstract. Investigation ої the ATP-ase' activity of erythrocyte membranes and free radical 
oxidation: in blood plasma ої patients: with: chronic tonzillite has been carried out before and 
after stochastic ultrasonic treatment. 

It has been revealed that stochastic ultrasonic therapy is effective in chronic tonzillite treatment 
producing a stabilizing effect on the ATP-ase activity of erythrocyte membranes through 
increasing or decreasing its altitude to the normal state. 

, After the treatment: is administered ten:times the increased level of free radical oxidation 
observed. in cases chronic, tonzillite is. dropped. down to normal by the stochastic ultrasound in 
patients’ blood plasma. 


Ультразвук (УЗ) широко використовується в медицині, зокрема B отоларин- 
гології, як безмедикаментозний неінвазійний метод лікування (1—6). 

Традиційні моночастотні УЗ терапевтичні апарати працюють на одній із 
стандартних, дозволених для вживання, частот від 0,88 МГц до 5,28 МГц. Водно- 
час моночастотне ультразвукове поле в біооб'єкті неоднорідне, що знижує тера- 
певтичний ефект. Крім того, клітини організму мають різні розміри, знаходяться 
в різних тканинах, які мають досить різні акустичні характеристики, що також 
знижує ефективність лікування моночастотним УЗ. У доповнення до всього, 
моночастотному УЗ притаманний небажаний фактор «звикання організму» при 
послідовній серії лікувальних сеансів, що також знижує ефективність терапії 
моночастотним УЗ, 

Зроблена спроба підвищити терапевтичну дію УЗ за рахунок використання 
принципово іншого ультразвукового : сигналу, в якому весь час випадково 
змінюються частота, амплітуда і фаза. Такий УЗ називається стохастичним, а 
тип дії на біооб'єкт -- інформаційним; При цьому клітини тканин сприймають з 
загального потоку найбільш: «підходящі» для них компоненти, маючи на увазі всі 
характеристики УЗ, тобто частоту, амплітуду і фазу. Ці складові загального УЗ 
сигналу будуть різними для різних клітин і тканин, і можуть викликати деякі 
резонансні ефекти, які змінюють стан цих клітин і тканин. 
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Лікувальна дія стохастичного ультразвуку (СТУЗ) була виявлена емпірично 
і теоретично не обгрунтованою, що стримує оптимізацію режимів ультразвуко- 
вої терапії при різних патологіях. В літературі існують поодинокі експеримен- 
тальні роботи по вивченню дії стохастичного ультразвуку (7--9), в той час як 
перспективи його використання потребують ретельного вивчення механізмів 
біологічної дії. 

Метою наших досліджень було вивчення ефективності застосування стоха- 
стичної ультразвукової терапії при лікуванні хронічного тонзиліту. Як відомо, 
одним 3 провідних механізмів розвитку запалення є активація вільнорадикаль- 
ного окислення та зв'язана з цим дискоординація ферментативної активності. 
Зважуючи на це, ми вивчили вплив стохастичної ультразвукової терапії на стан 
вільнорадикального окислення в плазмі крові та АТФ-азну активність мембран 
еритроцитів крові хворих на хронічний тонзиліт. 


СТОХАСТИЧНА УЛЬТРАЗВУКОВА ТЕРАПІЯ ХВОРИХ НА ХРОНІЧНИЙ 
ТОНЗИЛІТ 


Обстежено 25 хворих на хронічний декомпенсований тонзиліт. 

Для лікування використовували СТУЗ інтенсивністю 0,14 Вт/см? з випадко- 
вим зміненням частоти від 0,88 до 2,64 мГц. Режим роботи безперервний. Оп- 
ромінення робили по зовнішній методиці на область мигдаликів, по 5 хв на кож- 
ну. Джерелом СТУЗ був спеціально розроблений генератор з широкосмуговим 
аплікатором. 

У хворих на хронічний декомпенсований тонзиліт клінична картина визнача- 
лась за слідуючими ознаками: суб'єктивний стан почуття болю, відчуття комка в 
горлі; анамнестичні дані -- по частоті захворювань ангіною; результати об'єктив- 
ного дослідження -- фарингоскопія, стан регіонарних лімфатичних вузлів. 

Під час лікування СТУЗ спостерігається фазовість в його дії на реактивні 
прояви запалення в області піднебінних мигдаликів і піднебінних дужок. На по- 
чатку лікування у деяких хворих відмічалось погіршення самопочуття, загальна 
слабкість, які повністю проходили після 6-го сеансу. 

При фарингоскопії візуально спостерігались явища набряку і гіперемії в o6- 
ласті піднебінних мигдаликів і верхніх відділів піднебінних дужок, які після 
2-ro -- 3-го сеансу повільно посилювались і розповсюджувались до 5-го — 6-го 
сеансу, після чого відбувалась деяка стабілізація реактивних явищ. Після 6-ro — 
7-го сеансу вищеописані ознаки запалення слизової оболонки піднебінних мигда- 
ликів і піднебінних дужок зменшувались. Після 2-го сеансу посилювалась дре- 
нажна функція лакун, які до 6-го — 7-го сеансу практично повністю очищались. 
Після закінчення ультразвукової терапії піднебінні мигдалики значно зменшува- 
лись в розмірах порівняно з початковими даними. 


ВПЛИВ СТОХАСТИЧНОГО УЛЬТРАЗВУКУ НА ПЕРЕКИСНУ 
ХЕМІЛЮМІНЕСЦЕНЦІЮ ПЛАЗМИ КРОВІ 


Вільнорадикальне окислення в плазмі крові досліджували хемілю- 
мінесцентним методом на приладі ПХЛ — 01. Хемілюмінесценцію (ХЛ) інду- 
кували 3 %-ним перекисом водню. Кінетичні параметри свічіння оцінювали 
відповідно рекомендаціям [10]. | 

Результати дослідження перекисної хемілюмінесценції крові хворих на 
хронічний тонзиліт, а також при лікуванні стохастичним ультразвуком представ- 
лені у таблиці 1. 
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Нормалізуючий вплив стохастичного ультразвуку на хемілюмінесценцію плазми крові 


Таблиця 1 


Параметри хемілюмінесценції плазми крові хворих на хронічний тонзиліт 


Дослідні групи 


Умовно здорові 3,44 + 0,08 1,50 + 0,02 0,67 + 0,08 
Хронічний тонзиліт 3,81 + 0,09 1,50 + 0,04 0,67 + 0,06 
Терапія широкосмуговим 

ультразвуком 

Після 3-го сеансу 7,25 + 0,08 3,30 + 0,05 0,30 + 0,04 
Після 6-го сеансу 4,69 + 0,09 10,30 + 0,15 0,10 + 0,01 
Після 10-го сеансу 2,38 + 0,03 1,90 + 0,04 0,53 + 0,02 


eee и, 


Позначення: 5 — світлосума реакції хемілюмінесценції 
І -- амплітуда спалаху хемілюмінесценції 


При хронічному тонзиліті спостерігається деяке збільшення світлосуми ре- 
акції (5) хемілюмінесценції, що свідчить про збільшення інтенсивності окислення 
ліпідів. Амплітуда спалаху люмінесценції в групі хворих хронічним тонзилітом 
(D не відрізняється від такої в групі умовно здорових, тобто антиокислювальна 
активність (1/1) плазми крові не змінюється, що може бути викликане акти- 
вацією компенсаторних механізмів і формуванням патологічного окислювального 
гомеостазу. При лікуванні СТУЗ після 3-го і 6-ro сеансів відзначається збіль- 
шення показників як світлосуми, так і амплітуди спалаху хемілюмінесценції. 
Після десяти сеансів СТУЗ антиокислювальна активність знижується практично 
до початкових значень, а величина світлосуми -- до значень нижче рівня, від- 
міченого в групі умовно здорових осіб. 

Таким чином, стохастичний ультразвук після 10-ти сеансів знижує інтенсив- 
ність перекисного окислення плазми крові хворих на хронічний тонзиліт. 


ВПЛИВ СТУЗ НА АТФ-азну АКТИВНІСТЬ МЕМБРАН ЕРИТРОЦИТІВ 


Отримані дані вказують на вплив УЗ хвиль на структуру біологічних 
мембран, а також на структуру інтегральних мембранних білків, в тому числі 
АТФ-ази. Тому в роботі проводили дослідження активності K*, Ма" — АТФ-ази 
мембран еритроцитів у хворих хронічним тонзилітом, а також при проведенні 
стохастичної ультразвукової терапії. Визначення активності K*, Ма" — АТФ-ази 
в мембранах еритроцитів проводили спектрофотометричним методом (11). Мем- 
брани еритроцитів виділяли методом диференційного центрифугування. Кількісно 
препарати характеризували за наявністю в них білку. 

Встановлено, що показники АТФ-азної активності мембран еритроцитів хво- 
рих значно відрізняються від показників умовно здорових осіб. (таблиця 2). Ак- 
тивність досліджуваного ферменту може бути як підвищеною (1 група), так і 
пониженою (2 група), що можна пояснити, певною мірою, різним ступенем ви- 
раження запального процесу. Зміна активності Ме?" — АТФ-ази i К+, Ма" — 
АТФ-ази мембран еритроцитів свідчить про порушення транспорту іонів через 
мембрану і порушення електрохімічного іонного градієнту. 


Таблиця 2 
АТФ-азна активність мембран еритроцитів при хронічному тонзиліті і СТУЗ терапії 


азна активність, мкмоль Рп/ч на мг білк 


Групи людей 


р — АТФ-аза 
Умовно здоров! 0,0560 0,0270 0,0295 
+ 0,0020 + 0,0018 + 0,0051 
Хворі на хронічний тонзиліт 
І група 0,209 0,1220 0,0870 
+ 0,056 + 0,0080 + 0,0020 
П група 0,0019 0,0001 0,0018 
+ 0,0008 + 0,00002 + 0,0009 
СТУЗ терапія 
І група 0,0480 0,0200 0,0280 
+ 0,0013 ж 0,0009 + 0,0013 
П група 0,0460 0,0280 0,0180 
+ 0,0010 + 0,0015 + 0,0015 
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Зміна внутрішньо- і зовнішньоклітинної концентрації К”, Ма", Ca?*, які 
являють собою промотори більшості біохімічних реакцій, не може не впливати на 
нормальне їх протікання і дає змогу розвитку патологічного процесу. 

При впливі СТУЗ через 10 сеансів активність К", Ма" -- АТФ-ази 
змінюється і спрямованість цих змін залежить від початкового рівня активності 
досліджуваного ферменту. У хворих хронічним тонзилітом з підвищеною (1 гру- 
па) по відношенню до групи умовно здорових осіб АТФ-азною активністю спо- 
стерігається зниження загальної, Mg?* — залежної і оубаінінгибуючої АТФ-азної 
активності, У хворих П групи зі зниженими, по відношенню до норми, досліджу- 
ваними показниками -- їх підвищення. 

Таким чином, стохастичний УЗ нормалізує як підвищені, так і знижені при 
хронічному тонзиліті показники АТФ-азної активності мембран еритроцитів. 

Аналізуючи отримані дані, можна зробити такі висновки: 

1. Стохастична ультразвукова терапія виявляється ефективною при ліку- 
ванні хронічного тонзиліту. 

2. Стохастичний ультразвук нормалізує як підвищені, так і знижені при хро- 
нічному тонзиліті показники АТФ-азної активності мембран еритроцитів. 

3. Після десяти сеансів стохастичний ультразвук знижує підвищений при 
хронічному тонзиліті рівень вільнорадикального окислення в плазмі крові хворих. 


НОРМАЛІЗУЮЧИЙ ВПЛИВ СТОХАСТИЧНОГО УЛЬТРАЗВУКУ НА ХЕМІЛЮМІНЕСЦЕНЦІЮ 
ПЛАЗМИ КРОВІ ТА АТФ-азну АКТИВНІСТЬ ЕРИТРОЦИТІВ КРОВІ ХВОРИХ НА 
ХРОНІЧНИЙ ТОНЗИЛІТ 


КАРАСЬ A. Ф., АНДРІЙЧЕНКО С. B., НУРИЩЕНКО Н. Є, СТРИЖАК B. В., ОСТРОВСЬКИЙ L B. 


Проведено дослідження АТФ-азної активності мембран еритроцитів та вільнорадикального окислення 
в плазмі крові хворих на хронічний тонзиліт до та після лікування стохастичним ультразвуком. 


Встановлено, що стохастична ультразвукова терапія виявляється ефективною при лікуванні хронічного 
тонзиліту, нормалізуючи як підвищені, так і знижені при цьому захворюванні показники АТФ-азної 
активності мембран еритроцитів. 


Після десяти сеансів стохастичний ультразвук знижує підвищений при хронічному тонзиліті рівень 
вільнорадикального окислення в плазмі крові хворих. 


НОРМАЛИЗУЮЩЕЕ ВЛИЯНИЕ СТОХАСТИЧЕСКОГО УЛЬТРАЗВУКА НА 
ХЕМИЛЮМИНЕСЦЕНЦИЮ ПЛАЗМЫ КРОВИ И АТФ-азную АКТИВНОСТЬ 
ЭРИТРОЦИТОВ КРОВИ БОЛЬНЫХ ХРОНИЧЕСКИМ ТОНЗИЛЛИТОМ 


КАРАСЬ А. Ф., АНДРЕЙЧЕНКО С. В., НУРИЩЕНКО Н. Е., СТРИЖАК В. В., ОСТРОВСКИЙ И. В. 


Проведено исследование АТФ-азной активности мембран эритроцитов и свободнорадикального окис- 
ления в плазме крови больных хроническим тонзиллитом до и после лечения стохастическим ультра- 
звуком. 


Установлено, что стохастическая ультразвуковая терапия выявляется эффективной при лечении хро- 
нического тонзиллита, нормализуя как повышенные, так и пониженные при этом заболевании показа- 
тели АТФ-азной активности мембран эритроцитов. 


После десяти сеансов стохастический ультразвук понижает повышенный при хроническом тонзиллите 
уровень свободнорадикального окисления в плазме крови больных. 
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EVENTS 


UKRAINIAN CONGRESS OF RADIOLOGISTS HELD IN COMMEMORATION 
OF THE CENTENARY OF X-RAYS DISCOVERY 


Valeri E. Orel, Physics-Technical Laboratory, Ukrainian Research Institute of Oncology and Radiology. 


November 8, 1995 is the centenary of X-rays discovery by Wilhelm Conrad 
Róntgen (1845—1923). Ukrainian medical radiation physics and engineering have 
glorious international achievements in the use of X-rays in medicine. A physicist and 
electrical engincer J. Puluj (1845—1919) who was born in Ukraine was employed as 
Professor ai Prague German Technical! University and had invented an original 
aea ray iube, awarded a silver medal at the Paris international exhibition of 


1881. 1896 Puluj obtained an X-ray im age of a guinca- -pig, a hand and ctc., using 
an x ray Ws i11 а S. Ovoshnikov (1907—1 1994) veloped original X-ray 
equipment anc astrumenis for fil m ing human 15 EI :ngth at a reduced scale, He was 

in | 1951 of the Ukrainian 


awarded Stal 


Prise 
4 blyum was development of the 
ography as carly as 1957 [3]. These pape 
ore the well-known Nobel prise winning papers 


S 


| unian medical ra 
singled out among the besi ones during the on 
medicine in a report bi dirofessur F. Nüsslin à 
Congress in Würzburg (Germany) held on Scptember 20— 23, 1995. A1 this Сс 53 
the Association of Ukrainian medical physicists was awarded a membership in the 
European Federation of Organizations for Medica! Physics. 

emoration of “the centenary of X-rays discovery by W. С. Róntgen on 
May 1 | 5 the venue for the Ukrainian Congress of 1 
A section on the physics of medical imaging and its technology was held 
framework of this Congress. The section was attended by medical phy 
engineers and physicians from pat Russia, Germany, U SA and Hungary, as well 
as students trained ai Kiev unive sities in medical physics and engineering. Lectures 
were delivered by Professor У, E, "Orel on «The history, current role and outlook for 
medical phy in Ukraine», Dr. L. Z. Kalmykov «Dosimetric systems for medicine», 
Dr. W. Huda «Archives and communication systems in radiology», Professor 
K. S. Ternovoy «Modern high technologies of imaging using computer X-ray and 
magnetic resonance tomography». The lectures and reports evoked heated discussions. 


iation physicists and engineers were 
hundred vears of using X-rays in 
ermany) at ihe Ronigen 
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INTRODUCTION TO QUANTUM MEDICINE by S. P. Sit'ko and.L. М. Mkrtcian 
(Kiev, 1994). 


The monograph is devoted to the description of fundamental physical 
foundations of quantum medicine and the applied aspects of this principally new 
direction in the context of traditional medicine problems and achivements. 

The decade experience of application of methods of microwave resonance therapy 
elaborated by the collective of researchers and consultants of SRC «Vidhuk» 
introduction in different clinics of Ukraine, also the former USSR is summarized here. 

The results of treatment of the diseases of digestive organs, blood circulation, 
respiration, osteoarticular, nervous system, mental deviations, malignant neoplasias 
are described. In all nessery cases the references to' the perfomers of research are 
presented. 

The monograph is addressed to physicians, biologists and researches in arca of 
physics of the Alive. 


Монографію присвячено розгляду фундаментальних фізичних засад кванто- 
вої медицини та прикладним аспектам цього принципово нового напрямку в кон- 
тексті проблем і досягнень традиційної медицини. 

Узагальнено 10-річний досвід впровадження методик мікрохвильової резо- 
нансної терапії, які створено колективом співробітників і консультантів МНІЦ 
«Відгук», у різних клініках України, а також колишнього Радянського Союзу. 

Описано результати лікування органів травлення, кровообігу, дихання, руху, 
нервової системи, психічних відхилень, злоякісних новоутворень. У всіх не- 
обхідних випадках даються посилання на виконавців. 

Монографію розраховано на медиків, біологів, біофізиків, а також 
дослідників у галузі фізики живого. 


The book was published by SRC «Vidhuk» and сап be purchased directly from 
it at the following address: SRC «Vidhuk», Volodymyrska st. 61-B, Kiev 252033, 
Ukraine. 
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ІНФОРМАЦІЯ ДЛЯ АВТОРІВ 


Головним призначенням журналу є публікація статей та повідомлень, які демонстру- 
ють тенденцію переходу від традиційної біофізики до «фізики живого». Цей перехід грун- 
тується на застосуванні синергетичних та квантових принципів в дослідженні живого на 
всіх рівнях його ієрархічної організації. 


Журнал складається з трьох рубрикацій. У рубриці «Синергетика та квантова фізика 
живого» будуть друкуватися статті загального характеру з фундаментальних проблем си- 
нергетики, квантової механіки та їх можливого взаємозв'язку. Рубрика «Біофізика» при- 
свячується дослідженню живого методами традиційної біофізики, У рубриці «Медичні та 
біологічні застосування» будуть друкуватися статті прикладного характеру, які мають 
відношення до тематики журналу. 


Статті публікуються англійською, російською та українською мовами у відповідності 
з обраною в рукописі мовою. Обсяг статті має не перевищувати одного друкованого 
аркуша. 


ЗАПРОШУЄМО ДО СПІВРОБІТНИЦТВА! 
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INFORMATION FOR AUTHORS 


The journal is specially designed to publish articles and reports that show the tendency of 
transition from traditional Biophysics to Physics of the Alive. This transition is based on the use 
of synergetical and quantum principles in studying life at all levels of its hierarchic organisation. 


The journal contains three main sections. Under the «Synergetics and Quantum Physics of 
the Alive» heading we will publish general articles on problems of synergetics, quantum 
mechanics and their possible interrelation as well as their application to the Alive. The 
«Biophysics» section will be devoted to the traditional issues of this science. Under the heading 
of «Medical and biological applications» we will publish articles of applied character, relating to 
the Journals subjects. 


According to the language of the original, articles will be published in English, Ukrainian 
and Russian. The size of an article cannot exceed 1 printer's sheet. All articles are to be 
accompanied by abstracts. 


You are welcome to cooperation! 


ИНФОРМАЦИЯ ДЛЯ ABTOPOB 


Главным назначением журнала является публикация статей и сообщений, выявляю- 
щих тенденцию перехода от традиционной биофизики к «физике живого». Этот переход 
основан на применении синергетических и квантовых принципов в исследовании живого 
на всех уровнях его иерархической организации. 


Журнал содержит 3 основные рубрики. В рубрике «Синергетика и квантовая физика 
живого» будут печататься статьи общего характера по фундаментальным проблемам си- 
нергетики, квантовой механики и их возможной взаимосвязи. Рубрика «Биофизика» будет 
посвящена методам традиционной биофизики в исследовании живого. В рубрике «Меди- 
цинские и биологические приложения» будут печататься статьи прикладного характера, 
имеющие отношение к тематике журнала. 


Статьи будут публиковаться на английском, русском и украинском языках в соответ- 
ствии с языком рукописи. Объем статьи — до одного печатного листа. 


ПРИГЛАШАЕМ К СОТРУДНИЧЕСТВУ! 
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Нормалізуючий вплив стохастичного ультразвуку на 
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Ukrainian Congress of Radiologists Held in Commemoration of the 
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